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ABSTRACT 


With  the  objective  of  determining  optimum  design  criteria 
for  predetection  multichannel  filters  for  use  at  the  subarray  level  at  the 
LASA,  a  study  was  conducted  of  the  characteristics  of  the  ambient  noise 
fields  and  instrument  responses  of  five  selected  subarrays. 

Characteristics  of  the  ambient  noise  fields  which  relate  to 
design  of  multichannel  filters  and  which  were  studied  include 

•  Coherence  as  a  function  of  distance  between 
seismometers  as  well  as  a  function  of  frequency 

•  Power  density  as  a  function  of  frequency  and 
wavenumber 

•  Isotropicity  of  noise  fields 

•  Similarity  of  noise  fields  at  widely  separated 
subarrays 

•  Similarity  of  daytime  noise  fields  and  nighttime 
noise  fields 

Characteristics  of  individual  instrument  responses  studied 

include 


•  Variation  among  elements  of  a  particular  subarray 

in  amplitude  and  phase  response  at  selected  frequencies 

•  Variation  with  time  in  amplitude  and  phase  response 
of  each  instrument  at  1  cps 

•  Degree  to  which  phase  responses  can  be  equalized  by 
minimum -phase  equalization  filters 

Based  on  the  measured  characteristics  of  the  noise  fields 
and  instrument  responses,  14  multichannel  filters  were  developed  and 
compared  as  to  degree  of  improvement  in  signal -to -nois  e  ratio.  These 
filters  were  designed  in  such  a  way  that  only  one  design  criterion  differed 
for  each  pair  compared,  thus  allowing  evaluation  of  that  criterion's  effect 
on  filter  performance. 

A  multichannel  filter,  of  a  design  determined  by  the  com¬ 
parison  experiments  to  be  optimum  for  subarray  predetection  filtering, 
was  developed  for  each  of  five  subarrays.  These  filters  were  evaluated 
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as  to  frequency -wavenumber  response,  attenuation  of  ambient  noise  and 
signal  distortion.  Five  multichannel  filters  of  a  slightly  different  design 
also  were  developed.  These  filters  achieve  greater  signal -to -noise  im¬ 
provement  in  the  case  of  signals  with  horizontal  velocity  components  above 
20  km/sec. 


These  multichannel  filters,  and  certain  others  of  interest, 
were  punched  on  paper  tape  for  on-line  evaluation  with  the  Digital  Multi¬ 
channel  Filter  System  at  LASA. 

A  processing  procedure  for  detecting  signals  at  either  sub¬ 
array  or  full-array  level  is  recommended  and  illustrated  by  application 
to  selected  events  at  those  subarrays  under  consideration. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 


A.  INTRODUCTION 

The  objective  of  the  study  conducted  under  MIT  Lincoln 
Laboratory  Subcontract  325,  AF  19  (628)-5l67  Pr  ime,  was  the  specification 
of  design  criteria  for  optimum  predetection  multichannel  filters  to  be  used 
at  the  subarray  level  on  the  Large  Aperture  Seismic  Array  (LASA).  Certain 
characteristics  of  the  ambient  noise  fields  and  the  instrument  responses,  and 
the  variability  of  both  factors  with  time  were  investigated  since  they  bear 
directly  on  the  problem  of  optimum  filter  design.  Multichannel  filters  of  the 
type  determined  to  be  optimum  for  predetection  processing  were  developed 
for  five  subarrays  and  were  evaluated  through  a  comparison  of  response 
characteristics,  noise  attenuation  and  lack  of  signal  distortion. 

A  number  of  both  daytime  and  nighttime  ambient  noise  samples 
were  used  to  study  and  compare  noise  fields  and  to  develop  measured-noise 
multichannel  filters.  Routine  1  -  cps  calibrations  spanning  one  month  were 
employed  to  study  the  time  -  stability  of  instrument  responses,  and  a  special 
multifrequency  calibration  was  employed  to  study  variation  in  response 
among  the  instruments  within  selected  subarrays.  The  remainder  of  the 
data  ensemble  consisted  of  recordings  of  several  teleseismic  P-waves  which 
were  used  for  filter  evaluation.  Direct  copies  of  LASA  field  tapes  were 
used  for  the  entire  data  ensemble. 

Since  the  objective  was  the  determination  of  the  optimum  type 
of  MCF  for  predetection  processing  and  was  not  the  actual  development  of 
such  filters  for  all  21  subarrays,  only  five  of  the  subarrays  were  selected 
for  study  and  filter  development.  These  five  were  selected  to  give  the 
maximum  separation  between  subarrays  and  hence  a  good  representation 
of  the  differing  conditions  across  the  LASA.  This  number  of  subarrays  was 
considered  to  be  both  economical  and  adequate  for  the  purposes  of  the 
investigation.  Demultiplexing  to  individual  records  of  each  channel  to  be 
processed  was  necessary  and  desirable  due  to  the  multiplexed  format  in 
which  the  data  was  received  and  the  small  number  of  subarrays  to  be  pro¬ 
cessed.  The  bulk  of  data  processing  was  accomplished  on  the  IBM  7044 
computer,  with  a  small  amount  of  filter  application  and  detection  processing 
accomplished  on  the  TIAC*  (Texas  Instruments  Automatic  Computer). 
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While  the  data  ensemble  was  being  collected  and  prepared 
for  analysis,  the  wavenumber  resolution  of  a  LASA  subarray  was  measured 
by  designing  two  multichannel  filters  from  theoretical  noise  and  signal 
models.  In  one  case  the  noise  was  modeled  as  a  disk  in  the  wavenumber 
plane  and  the  signal  as  a  point  source  from  which  energy  propagated  from 
North  to  South  across  the  subarray  at  an  apparent  horizontal  velocity  of 
18  km/sec.  In  another  case  the  noise  was  modeled  as  an  annulus  in  the 
wavenumber  plane  and  the  signal  as  having  an  infinite  apparent  horizontal 
velocity.  These  two  filters  were  fully  described  and  analyzed  in  Special 
Report  No.  1  and  Special  Report  No.  2¥ .  After  an  inspection  of  the  filter 
responses,  the  conclusion  was  that  directional  signal  models  and  distributed 
signal  models  were  justified  by  the  subarray  wavenumber  resolution  below 
2. 0  cps. 

Optimum  design  criteria  for  subarray  predetection  multi¬ 
channel  filters  were  determined  by  designing  and  comparing  a  number  of 
filters  for  one  subarray.  Only  one  design  parameter  differed  for  the  pair 
of  filters  being  compared  in  each  case,  thereby  allowing  an  exact  assessment 
of  the  effect  on  filter  performance  of  that  design  parameter. 

Multichannel  filters  were  developed  for  five  subarrays,  according 
to  the  criteria  determined  by  the  comparisons  previously  made,  to  be  optimum 
for  predetection-filtering  for  12-km/sec  and  higher  velocity  signals.  These 
filters  were  evaluated  by  computing  their  transfer  functions  and  wavenumber 
responses  and  by  measuring  attenuation  of  ambient  noise  as  a  function  of  fre¬ 
quency.  Distortion  of  teleseismic  signals,  or  lack  thereof,  was  estimated 
by  applying  each  filter  to  several  signals  and  comparing  with  both  a  single 
seismometer  and  an  average  trace. 

A  detection  processing  procedure  was  illustrated  by  application 
to  three  teleseismic  signals  at  each  of  the  five  subarrays.  This  procedure 
consisted  of  squaring  the  output  of  a  Wiener  MCF  and  integrating  over  a  moving 
time  gate.  Earlier  work  on  another  project  has  indicated  this  procedure  to 
be  equally  as  good  as  the  analysis -of-variance  procedure  or  the  quadratic 
processor  procedure.  Since  the  technique  just  illustrated  is  more  easily  and 
economically  implemented  than  the  other  two  procedures,  it  is  recommended 
for  use  in  signal  detection. 


Texas  Instruments  Incorporated,  1966:  LASA  Data  Analysis  and  MCF  Support, 
Special  Rpt.  No.  1,  Feb  24. 

Texas  Instruments  Incorporated,  1966:  LASA  Data  Analysis  and  MCF  Support, 
Special  Rpt.  No.  2,  Apr. 

Texas  Instruments  Incorporated,  1966:  Array  Research  Final  Rpt.  ,  Matrix  - 
Multiply  Detection  Processing  of  Array  Data,  May  24. 
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Included  in  this  report  are  two  appendixes  containing  lists  of 
filter  weights  for  all  multichannel  filters  developed  in  the  course  of  this 
project  and  lists  of  those  average  correlation  functions  used  in  designing 
the  five  final  filters. 

The  five  optimum  filters  designed  to  pass  all  signals  of 
apparent  horizontal  velocity  greater  than  12  km/sec  were  punched  into 
paper  tape  for  input  to  the  on-line  digital  MCF  system.  Additionally,  five 
multichannel  filters  of  exactly  the  same  design  as  the  disk  signal  model 
filters,  except  for  the  use  of  infinite  velocity  signal  models  in  design,  were 
also  supplied  on  punched  tape.  These  filters  provide  better  noise  attenuation 
at  the  expense  of  some  signal  attenuation  (above  1  cps)  below  20  km/sec  in 
apparent  horizontal  velocity.  MCF  11  and  MCF  12  were  also  punched. 

B.  AMBIENT  NOISE  CHARACTERISTICS 

Although  a  complete  characte rizaton  of  LASA's  ambient  noise 
field  was  not  the  intent  of  this  project,  certain  characteristics  of  the  noise 
field  were  investigated.  Those  characteristics  investigated  were  especially 
significant  in  determining  the  optimum  type  of  predetection  processing. 

Noise  coherence  as  a  function  of  frequency  and  distance  was 
studied  first  since  this  characteristic  determines  whether  multichannel  filters 
can  achieve  greater  noise  reduction  than  a  simple  beam-steer  operation.  If 
the  noise  in  the  signal  frequency  band  were  not  coherent  between  seismo¬ 
meters,  then  beam-steering  in  combination  with  optimum  Wiener  frequency¬ 
filtering  would  provide  all  the  S/N  improvement  possible.  The  ambient  nois e  was 
found  to  be  coherent  to  a  significant  degree  in  the  signal  frequency  band 
(0.5  cps  to  3.0  cps).  Coherence  was  measured  by  developing  and  applying 
both  single-channel  and  multichannel  spatial  prediction  filters  and  measuring 
the  percentage  of  the  noise  predictable  as  a  frequency  function. 

Coherence  generally  fell  off  as  distance  increased  between 
the  seismometer  whose  output  was  predicted  and  the  seismometer  whose 
output  was  filtered.  However,  some  noise  in  the  signal  band  was  still 
predictable  at  distances  of  3.5  km,  the  subarray  radius. 

In  Section  II  the  results  of  the  noise  characterization  are 
presented  in  detail.  In  addition  to  spatial  coherence  measurements,  esti¬ 
mates  of  the  stationarity  of  the  coherent  noise  component  are  presented 
through  comparisons  of  the  percent  predictable  at  the  time  of  filter  design 
and  the  percent  predictable  six  days  later.  Frequency-wavenumber  spectra 
of  representative  daytime  and  nighttime  noise  samples  are  presented. 
Comparison  of  day  and  night  noise  has  indicated  that  multichannel  filters 
designed  from  an  average  of  daytime  and  nighttime  noise  correlations  would 
perform  satisfactorily  and  that  specific  day  and  night  filters  would  be  un- 
neces  sary. 
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Power  density  spectra  of  noise  recorded  on  the  center  seis¬ 
mometer  500  ft  deep  and  noise  recorded  on  the  seismometers  ZOO  ft  deep 
indicated  that  the  noise  was  identical  below  1.5  cps.  Above  1.5  cps  the 
center  seismometer  recorded  less  noise  than  did  the  shallower  seismo¬ 
meters;  this  difference  increased  almost  linearly  with  the  frequency  to 
7  db  at  5 .  0  cps . 

C.  VARIATIONS  IN  INSTRUMENT  RESPONSE 

Differences  in  response  characteristics  of  the  seismometer- 
amplifier  package  from  instrument  to  instrument  are  of  interest  in  the  MCF 
design  problem  because  if  uncompensated,  these  differences  can  cause  serious 
signal  attenuation.  This  signal  attenuation  is  due  to  the  fact  that  theoretical 
and,  therefore,  perfectly  equalized  signal  models  are  usually  used  in  filter 
development.  Two  solutions  are  readily  available,  however:  s  ingle  -  channel 
filters  may  be  designed  from  multifrequency  calibration  data  to  equalize  the 
response  of  each  instrument  to  some  reference,  or  the  response  variations 
as  determined  from  the  calibration  analysis  may  be  incorporated  into  the 
signal  model.  The  latter  method  is  more  easily  accomplished  and  was  used 
for  the  development  of  the  response-compensated  MCF  (Section  IV  ,  para- 
graph  C). 

To  test  the  feasibility  of  compensating  for  differences  in 
instrument  response  by  this  method  and  to  determine  the  seriousness  of 
the  differences,  a  special  multifrequency  calibration  of  each  of  five  sub¬ 
arrays  was  conducted.  Analysis  of  this  calibration  is  presented  in  Section 
III-B.  The  differences  in  both  phase  and  amplitude  response  were  found 
to  be  large  at  all  frequencies  except  near  1.  0  cps.  Compensation  for  these 
differences  would  clearly  improve  broadband  signal  response  of  the  pre¬ 
detection  multichannel  filters. 

Another  important  instrument  response  characteristic  to 
MCF  design  is  the  instrument’s  variation  with  time.  If  there  is  significant 
variation  with  time,  multichannel  filters  must  be  redesigned  at  frequent 
intervals.  Nine  calibrations  (at  1.0  cps)  of  each  of  the  125  instruments  being 
considered  were  obtained  and  analyzed  to  estimate  the  time  -  stationarity  of 
each  instrument’s  response.  These  nine  calibrations  covered  one  month 
and  during  this  period  each  instrument  was  found  to  have  a  very  stable 
response,  except  for  occasional  sudden  changes  for  individual  instruments. 
These  sudden  changes  in  response  are  largely  unexplained,  although  from 
their  nature  it  is  believed  that  they  are  due  primarily  to  adjustments,  repairs 
or  replacements  of  the  instrument  components.  Several,  but  not  all,  instru¬ 
ments  within  a  subarray  would  simultaneously  experience  these  sudden  re¬ 
sponse  changes  which  usually  better  equalized  those  instruments  whose  ampli¬ 
tude  responses  at  1.  0  cps  deviated  most  from  unity. 
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These  step-like  changes  in  instrument  response  were  very 
numerous  during  November  1965,  There  was  such  poor  correlation 
for  many  instruments  between  the  responses  at  the  end  of  November  and 
the  1-cps  response  measured  in  March  1966  when  the  multifrequency  cali¬ 
brations  were  run  that  it  was  deemed  impractical  to  include  compensation 
for  instrument  response  variations  in  the  design  of  the  five  multichannel 
filters  to  be  supplied  for  the  on-line  digital  system. 

An  evaluation  of  an  equalization  technique  which  models  the 
response  variation  among  instruments  as  the  effect  of  a  minimum -phas e 
filter  was  conducted  using  a  multifrequency  calibration  test  recorded  with 
the  Texas  Instruments  Digital  Field  System  (DFS)',C  during  May  1965.  The 
technique  consists  of  designing  for  each  channel  a  minimum -phas e  whitening 
filter  from  the  ratio  of  the  amplitude  response  of  that  channel  to  the  ampli¬ 
tude  response  of  an  arbitrary  reference  channel.  Evaluation  primarily  was 
concerned  with  the  degree  of  phase  response  equalization  obtained  with  this 
method.  Good  improvement  was  obtained  at  all  frequencies,  the  greatest 
improvement  being  near  the  1-cps  damped  resonant  frequency  of  the  seis¬ 
mometer.  Variations  of  ±20  were  typically  reduced  to  less  than  ±5°. 

D.  DETERMINATION  OF  OPTIMUM  DESIGN  PARAMETERS  FOR 
PREDETECTION  MULTICHANNEL  FILTERS 

Thirteen  multichannel  filters  were  designed  and  evaluated  in 
order  to  determine  the  effects  of  each  design  parameter  on  the  S/N  ratio 
improvement  capability  of  each  filter.  These  filters  were  designed  so  that 
only  one  design  parameter  differed  for  each  pair  to  be  compared.  The  pairs 
of  filters  were  compared  as  to  transfer  function,  frequency-wavenumber  re¬ 
sponse,  ambient  noise  attenuation,  and  signal  distortion  or  attenuation. 
Signal  distortion  was  not  considered  an  undesirable  effect  as  was  signal 
attenuation  since  the  filters  were  to  be  used  only  for  predetection -filte ring. 
The  specific  design  parameters  whose  effects  were  compared  were 

•  Ring-model  (5  and  8)  multichannel  filters  vs 
25-channel  maximum -likelihood  multichannel 
filte  rs 

•  Use  of  an  infinite  velocity  signal  model  vs 
use  of  a  disk  signal  model 

•  Simulation  of  instrument  response  differences 
in  the  signal  model  vs  an  equalized  signal  model 
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•  Variable  S/N  ratio  (permitting  filtering  on 
a  frequency  basis  as  well  as  on  a  velocity 
basis)  vs  a  constant  S/N 

•  Use  of  measured-noise  statistics  vs  use  of 
a  theoretical  noise  model 

•  1-sec  filters  vs  2-sec  filters  vs  6-sec  filters 

•  Whitened  measured  noise  vs  nonwhitened  measured 
nois  e 

•  Use  of  daytime  noise  statistics  vs  use  of  nighttime 
noise  statistics 

Results  of  these  comparisons  are  presented  and  discussed  in  detail  in 
Section  IV. 

E.  DEVELOPMENT  AND  EVALUATION  OF  PREDETECTION  MULTI¬ 
CHANNEL  FILTERS 

On  the  basis  of  the  comparisons  of  various  design  parameters, 
an  MCF  considered  to  be  optimum  for  predetection-filtering  at  the  LASA 
subarray  level  was  developed  for  each  of  the  five  subarrays. 

Five-channel  ring-model  filters  were  developed  in  each  case. 
The  noise  models  consisted  of  weighted  and  averaged  daytime  and  nighttime 
noise  correlations  which  had  been  whitened.  The  noise  was  whitened  by 
convolving  the  measured  correlation  set  with  the  autocorrelation  of  a  5-point 
whitening  filter  developed  from  the  average  of  all  autocorrelation  functions. 

Noise  correlations  were  computed  for  ±40  lags,  but  the 
whitening  operation  limited  the  available  length  to  ±36  lags.  For  this 
reason  all  filters  are  37  points,  or  1.85-sec,  long. 

The  signal  models  simulated  a  disk  in  wavenumber  space 
with  a  radius  corresponding  to  12  km/sec  and  were  band-limited  from 
1.0  cps  to  3.0  cps .  Differing  instrument  responses  were  not  simulated 
in  the  signal  models  since  comparison  of  November  1965  and  March  1966 
calibrations  at  1.  0  cps  had  indicated  significantly  different  responses  for 
many  of  the  instruments.  Also,  the  filter  comparison  discussed  in  Section 
IV,  paragraph  C.  7  showed  that  a  filter  incorporating  relative  responses 
estimated  from  March  calibrations  into  the  signal  model  yielded  a  S/N 
ratio  improvement  2  or  3  db  less  than  a  filter  designed  from  an  equalized 
signal  model  when  both  filters  were  applied  to  November  events.  A  total 
signal  energy  to  total  noise  power  ratio  of  4  was  used  for  all  filter  design. 
Since  the  noise  was  prewhitened  and  the  signal  model  was  band -limited , 
the  effective  S/N  ratio  in  the  signal  frequency  band  was  actually  nearer 
to  10  or  12. 
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These  filters  were  evaluated  through  computation  of  transfer 
functions  and  frequency-wavenumber  responses  and  through  comparison  to 
straight  summation  as  to  ambient  noise  attenuation.  Filter  response  was 
always  within  ±1.  5  db  of  unity  at  velocities  greater  than  12  km/ sec  in  the 
frequency  band  of  1.  0  to  3.  0  cps.  In  addition  to  frequency- filte ring  outside 
the  signal  band,  these  filters  can  achieve  considerable  noise  attenuation  on 
a  velocity  basis  . 

The  wavenumber  responses  of  all  filters  were  remarkably 
similar.  Since  average  noise  statistics  were  used  in  filter  design,  this 
similarity  in  response  is  an  indication  that  the  long-term  average  noise 
statistics  are  very  similar  across  the  LASA.  At  1.5  cps  and  above,  all 
responses  were  uniformly  low  outside  the  12-km/sec  disk,  as  would  be 
expected  if  the  majority  of  the  noise  at  these  frequencies  was  incoherent. 

Attenuation  of  ambient  noise  in  the  signal  band  from  1 .  0  to 
3.  0  cps  ranged  from  3  to  6  db  less  than  that  obtained  by  averaging  all  25 
channels.  This  is  not  particularly  surprising  in  view  of  the  fact  that  the 
filters  are  constrained  to  pass  a  very  large  region  in  wavenumber  space  pre¬ 
cluding  the  \Zh~  improvement  in  S/N  attained  by  averaging  the  predominantly 
random  noise  in  this  frequency  band.  For  this  reason,  a  filter  identical  to 
the  MCF  in  all  respects,  except  that  the  disk  signal  model  was  replaced  with 
an  infinite  velocity  signal  model,  was  also  developed  for  each  of  the  five 
subarrays.  Previous  comparison  of  a  disk  signal  model  filter  with  an  in¬ 
finite  velocity  signal  model  filter  showed  that  the  latter  filter  achieved  up 
to  5  db  more  reduction  in  ambient  noise  in  the  signal  band.  Thus,  the  five 
infinite  velocity  filters  may  be  expected  to  achieve  about  3  db  more  noise  re¬ 
duction  in  the  signal  band  but  will  attenuate  signals  with  apparent  horizontal 
velocities  of  20  km/sec  or  more  by  less  than  1  db  at  1.  0  cps. 

Infinite  velocity  filters  would  be  preferred,  therefore,  as  pre¬ 
detection  filters  for  high-velocity  teleseismic  events  since  they  provide  more 
improvement  in  S/N  ratio.  Additional  improvement  could  probably  be  obtain¬ 
ed  by  using  directional  point  source  signal  models  if  the  objective  were  de¬ 
signing  filters  to  monitor  particular  geographical  regions. 

The  poorer  than  straight  summation  noise  attenuation  in  the 
signal  band  attained  with  infinite  velocity  signal  model  filters  strongly  suggests 
that  detection  at  the  LASA  subarray  level  is  primarily  limited  by  the  high- 
velocity  (above  8  km/sec)  mantle  P-wave  noise.  A  study  of  this  noise  in 
regard  to  its  composition,  time- stationarity ,  spatial  -  stationarity,  coherence, 
and  frequency-wavenumber  characteristics  is  now  underway  on  the  Large 
Array  Signal  and  Noise  Analysis  project  for  AFTAC.  It  is  expected  that,  when 
the  characteristics  of  the  high-velocity  noise  field  are  better  understood,  multi¬ 
channel  filters  which  will  achieve  a  significant  reduction  in  this  noise  compo¬ 
nent  can  be  devised  for  application  to  the  large  array. 
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The  signal  detection  procedure  recommended  is  a  square-and- 
integrate  operation  following  application  of  a  Wiener  MCF.  This  procedure 
is  illustrated  in  Section  V  by  applying  it  to  three  teleseismic  P-wave  signals 
at  each  of  the  five  subarrays.  A  relatively  constant  level  output  is  obtained 
preceding  the  P-wave,  and  an  abrupt  and  proportionately  very  large  increase 
in  the  output  is  seen  at  the  time  of  P-wave  arrival.  The  events  selected  for 
this  illustration  ranged  in  reported  magnitude  from  4.  6  for  the  Mariana  event 
to  5.  0  for  the  Algerian  event.  A  clear  indication  of  the  arrival  of  coherent 
teleseismic  energy  was  obtained  at  each  event. 
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SECTION  II 


AMBIENT  NOISE  CHARACTERISTICS 

A.  NOISE  COHERENCE  STUDIES 

1.  Purposes 

The  purposes  of  these  studies  were  to  determine  the  pre¬ 
dictability  of  EASA  noise,  on  a  subarray  basis,  in  the  seismic  band  of 
0.  0  to  5.  0  cps  as  a  function  of  both  distance  and  the  number  of  seismometers 
used,  and  to  obtain  a  preliminary  estimate  of  the  time -stability  of  the  noise. 

2.  Method 


Two  noise  samples  approximately  4  min  long  were  obtained 
for  subarray  F-3.  The  first  sample  (noise  sample  A)  was  recorded  on 
21  May  1965  at  07:07  MST  and  the  second  one  (noise  sample  B)  was  recorded 
on  26  May  1965  at  22:01  MST.  The  two  samples  were  decimated  by  4  to  give 
a  Nyquist  frequency  of  5.2  cps,  and  were  approximately  whitened  by  design¬ 
ing  a  whitening  filter  for  seismometer  32  (for  each  noise  sample)  and  ap¬ 
plying  that  filter  to  all  25  channels.  Figure  II- 1  shows  part  of  the  two  re¬ 
sampled  and  whitened  noise  samples. 

Two  sets  of  filters  to  predict  the  center  seismometer  were 
designed  for  noise  sample  A: 

•  Single  -  channel  filters  for  distances 
from  0.  5  to  3.  5  km 

•  Multichannel  filters  using  either 
rings  or  arms  of  the  array 

The  filters  are  described  in  Table  II- 1. 

All  filters  were  then  applied  to  noise  sample  A  and  almost 
all  of  the  filters  were  applied  to  noise  sample  B  (Table  II- 1).  The  error 
trace  was  obtained  as  the  reference  trace  (center  seismometer)  minus  the 
reference  trace  prediction.  Bartlett-smoothed  autocorrelations  of  the 
error  and  reference  traces  were  computed  and  Fourier-transformed  to 
obtain  spectral  estimates.  The  ratio  of  the  error  to  the  reference  spectrum 
was  obtained  to  measure  the  predictability  as  a  function  of  frequency  for 
each  filter. 
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Table  II- 1 


DESCRIPTION  OF  THE  PREDICTION  FILTERS 


Filter 

Description 

Applied 
to  NS  A 

Applied 
to  NS  B 

i 

Single  -  channel,  1/2  km  from  center 
s  eismometer 

Yes 

Yes 

2 

Single  -  channel ,  1  km  from  center 
s  eismometer 

Yes 

No 

3 

Single  -  channel ,  1-1/2  km  from  center 
s  eismomete  r 

Yes 

Yes 

4 

Single-channel,  2  km  from  center 
s  eismomete  r 

Yes 

No 

5 

Single  -  channel,  2-1/2  km  from  center 
seismometer 

Yes 

Yes 

6 

Single-channel,  3  km  from  center 
s  eismomete  r 

Yes 

Yes 

7 

Single -channel,  3-1/2  km  from 
center  seismometer 

Yes 

Yes 

8 

4-channel,  ring  2,  1/2  km  from 
center  seismometer 

Yes 

Yes 

9 

5 -channel,  rings  3  and  4,  1  and 

1-1/2  km  from  center  seismometer 

Yes 

Yes 

10 

5-channel,  rings  5  and  6,  2  and 

2-1/2  km  from  center  seismometer 

Yes 

No 

11 

5 -channel,  rings  7  and  8,  3  and 

3-1/2  km  from  center  seismometer 

Yes 

No 

12 

9-channel  "arm"  case  using  seis¬ 
mometers  1/2,  1-1/2  and  2-1/2  km 

from  center  seismometer  along  arms 

1,  3  and  5 

Yes 

No 
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B.  DISCUSSION  OF  RESULTS 

1.  Application  of  Prediction  Filters  to  Noise  Sample  A 
a.  Single -Channel  Prediction 

Figure  II-2  shows  the  reference  trace  and  the  reference  trace 
prediction  (noise  sample  A)  for  the  seven  single -channel  filters.  Two  obser¬ 
vations  can  be  made: 

•  The  predicted  output  had  much  less  high  frequency 
than  was  actually  present  which  implies  that  the 
high-frequency  noise  was  mainly  incoherent. 

Considering  the  noise  field  to  consist  of  coherent 
and  incoherent  noise,  then  the  frequency  response 
of  the  prediction  filter  will  be 

$  <0 

H(f)  =  i  (f)  +  $.(f) 

C  1 

where  §c  (f)  and  (f)  are  the  coherent  and  incoherent 
noise  sources,  respectively.  As  §c(f)-»0,  then  H(f)->  0 
and  the  filter  essentially  turns  off  or  becomes  a  band¬ 
pass  filter  which  rejects  the  high  frequency.  Since  the 
high-frequency  energy  was  not  passed  by  the  filters,  it 
must  have  been  mainly  incoherent. 

•  The  predictability  of  the  low-frequency  energy 
decreases  with  increasing  distance. 

Figures  II - 3  through  II - 5  show  the  error-to-reference  spectral 
ratios.  A  highly  predictable  peak  is  seen  at  0.  25  cps  for  which  the  pre¬ 
diction  decreased  as  a  function  of  distance  (indicated  in  Figure  II-2  and 
shown  in  Figure  II-6).  The  predictability  decreased  rapidly  to  about  0.6  cps 
and  then  decreased  more  slowly  to  essentially  no  predictability  at  about 
1.4  cps  for  seismometers  spaced  0.  5  km  apart.  At  3.  5  km  these  figures 
are  0.  4  and  1.0  cps,  respectively. 

At  higher  frequencies,  predictable  energy  existed  at  2.4,  2.8 
and  3.2  cps  and  to  some  extent  at  1.75  and  3.65  cps.  The  predictability 
of  these  peaks  also  decreased  with  distance;  the  2.4-cps  and  2.8-cps  peaks 
disappeared  at  1 . 5  km  and  the  3.  2-cps  peak  disappeared  at  2.  0  km. 
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Figure  II-2.  Center  Seismometer  and  Single-Channel  Predictions  of  Same  (NSA) 
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Figure  H-3.  Single-Channel  Prediction  of  the  Center  Seismometer  Using  Seismometer  25  (NSA) 
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Figure  II-4.  Single-Channel  Prediction  of  the  Center  Seismometer  Using 
Seismometers  36,  43,  52  (NSA) 
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Figure  II- 5 .  Single-Channel  Prediction  of  the  Center  Seismometer  Using 
Seismometers  61,  74,  85  (NSA) 


II- 8 


science  services  division 


II  -  9 


science  services  division 


Figure  II-6.  Predictability  of  0.  25-cps  Peak  as  a  Function  of  Distance  (NSA) 


b.  Multichannel  Prediction 

The  reference  trace  and  the  reference  trace  prediction  for  the 
four  ring  filters  and  the  9-channel  arm  filter  are  shown  in  Figure  II - 7 . 
Figures  II-8  and  II - 9  show  the  error-to- reference  spectral  ratios  for  the 
ring  filters  and  Figure  11-10  shows  the  ratio  for  the  arm  filter. 

Regarding  the  set  of  ring  filters 

•  The  predictable  peaks  seen  for  the  single- channel 
data  were  better  defined  by  the  multichannel  filters. 

Also,  additional  peaks  at  2.  1,  4.  1  and  4.7  cps 
were  observed. 

•  The  degree  of  predictability  was  a  function  of 
distance  again. 

•  The  predictions  were  better  by  2  to  3  db  than 
were  the  corresponding  single  -  channel  predictions. 

It  can  be  concluded  that  the  coherent  component  of  the  noise  is  either  multi- 
mode  or  multisource  energy  (or  both),  because  of  the  improved  prediction 
obtained  by  using  additional  channels  and  because  of  the  decrease  in  pre¬ 
dictability  with  distance. 

For  the  9-channel  filter  (Figure  II -  10)  the  er ror -to- r efe rence 
spectral  ratio  suggests  that  1-  or  2-db  prediction  has  been  obtained  at  all 
frequencies.  However,  this  prediction  probably  was  not  real,  but  due  to  the 
total  number  of  filter  points  becoming  large,  with  respect  to  the  length  of 
the  data;  i.  e.  ,  the  prediction  was  becoming  more  deterministic  than  statis¬ 
tical  due  to  the  fewer  degrees  of  freedom  in  filter  design. 

2.  Application  of  Prediction  Filters  to  Noise  Sample  B 
a.  Single  -  Channel  Prediction 

The  reference  trace  and  the  reference  trace  prediction 
(noise  sample  B)  for  five  single-channel  filters  are  shown  in  Figure  II- 11. 
Figures  11-12  and  11-13  show  the  error-to-reference  spectral  ratios.  Very 
little  noise  in  noise  sample  B  was  predictable,  using  the  filters  designed 
from  noise  sample  A.  At  0.25  cps  there  was  about  50-percent  prediction, 
but  at  other  frequencies,  except  for  the  0.5-km  filter,  there  was  no  pre¬ 
diction.  Above  about  2.  0  cps  there  was  actually  more  power  in  the  error 
than  in  the  reference  trace,  especially  for  the  close-in  filters. 
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Figure  H-7.  Center  Seismometer  and  Multichannel  Prediction  of  Same  (NSA) 
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FREQUENCY  (CPS) 

Figure  H-8.  Four-Channel  Prediction  of  the  Center  Seismometer  Using  Ring  2  (NSA) 


Figure  H-9.  Multichannel  Prediction  of  the  Center  Seismometer  Using 

Rings  3  and  4  (5  Channels),  5  and  6  (5  Channels),  and  7  and 
8  (4  Channels),  Respectively  (NSA) 
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Figure  II- 10.  Multichannel  Prediction  of  Center  Seismometer  Using 
9-Channel  "Arm"  Configuration  (NSA) 
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igure  XX— 11.  Center  Seismometer  and  Single-Channel  Prediction  of  Same  (NSB) 


Figure  11-12,  Single -Channel  Prediction  of  Center  Seismometer  Using 
Seismometers  25,  and  43,  Respectively  (NSB) 
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Figure  11-13.  Single-Channel  Prediction  of  Seismometer  Using 

Seismometers  61,  74,  and  85,  Respectively  (NSB) 
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The  apparently  poor  prediction  at  the  0. 25-cps  peak  was 
probably  partly  due  to  differences  in  the  relative  power  of  the  reference 
trace  and  traces  used  in  the  prediction.  As  an  example,  for  noise  sample 
A  the  0.5-km  seismometer  had  about  6  db  more  power  than  the  center  seis¬ 
mometer  at  0.25  cps,  but  for  noise  sample  B,  the  two  seismometers  had 
about  equal  power  (Figure  11-14).  Thus,  about  6  db  of  the  decrease  in  pre¬ 
dictability  for  noise  sample  B  could  be  explained.  Combining  this  with  the 
prediction  obtained  brings  noise  sample  A  and  noise  sample  B  into  reasonably 
good  agreement.  Adding  6  db  to  the  0.  25-cps  prediction  for  other  noise  sam¬ 
ple  B  distances  would  also  bring  these  sample  B  values  into  fair  agreement 
with  those  of  noise  sample  A.  (This  is  not  unreasonable  because  of  the 
difference  in  depth  between  the  center  seismometer  and  the  other  seismometer.) 
Thus,  it  appears  that  at  0.25  cps  the  noise  was  reasonably  time-stationary. 

Above  2.  0  cps  the  noise  does  not  appear  to  be  time  -  stationary. 
Figure  11-15  shows  the  coherence  between  the  center  seismometer  and  the 
0.5-km  seismometer  for  both  noise  sample  A  and  noise  sample  B.  There 
were  coherent  peaks  for  both  noise  samples,  but  they  occurred  at  different 
frequencies  for  each  sample.  Thus,  the  filter  designed  from  noise  sample  A 
would  be  ineffective  in  predicting  the  coherent  peaks  in  noise  sample  B.  Note 
that  noise  sample  A  was  recorded  at  07:07  MST  and  noise  sample  B  at  10:01 
MST  so  that  the  difference  in  the  high-frequency  energy  is  reasonable. 

b.  Multichannel  Prediction 

The  reference  trace  and  the  two  ring  predictions  of  the 
reference  trace  are  shown  in  Figure  11-16.  Figure  11-17  shows  the  error- 
to-reference  spectral  ratios.  These  data  corroborate  previous  observations 
that  the  0.  25-cps  peak  was  reasonably  time-stationary,  but  the  high-frequency 
(above  2.  0  cps)  noise  was  not.  Also,  the  improvement  obtained  over  single¬ 
channel  prediction  again  indicates  a  multimode  or  multisource  (or  both)  noise 
field  at  0.  25  cps. 

C.  CONCLUSIONS 

From  this  analysis  can  be  concluded  that  for  subarray  F-3 

•  There  was  a  peak  at  0.25  cps  which  was  both 
highly  predictable  and  fairly  time -stationary 

•  There  were  predictable  peaks  at  higher  frequencies 
(notably  2.4,  2.8  and  3.2  cps)  which  were  predictable 
but  not  time-stationary 

•  Where  the  noise  was  predictable,  it  was  either 
multimode  or  multisource  (or  both)  energy 
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Figure  11-14.  Comparison  of  Power  Density  at  0.25  cps  on 

Channels  1  and  4  for  NSA  and  NSB,  respectively 
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COHERENCE  COHERENCE 


Figure  11-15.  Coherence  Between  Channels  1  and  4  for  NSA  and  NSB, 
Respectively 
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Figure  11-16.  Center  Seismometer  and  Multichannel  Prediction  of  Same  (NSB) 
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Figure  11-17.  Multichannel  Prediction  of  center  Seismometer  Using  Ring  2 

(4  Channels)  and  Rings  3  and  4  (Channels),  Respectively  (NSB) 
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D.  COMPARISON  OF  DAYTIME  AND  NIGHTTIME  NOISE  STATISTICS 


Two  noise  samples  were  selected  for  use  in  comparing  day 
and  night  noise  statistics.  Noise  sample  D  was  recorded  at  the  F-3  subarray 
on  24  November  1965  between  03:21  and  03:24  MST.  Noise  sample  F  was 
recorded  at  the  F-3  subarray  on  24  March  1966  between  12:47  and  12:50  MST. 
The  frequency-wavenumber  spectrum  of  the  nighttime  noise  sample  D  is 
shown  in  Figure  11-18  for  frequencies  in  0.25-cps  increments  (from  0.25  cps 
to  1.  00  cps).  Figure  11-19  shows  the  f-k  spectrum  at  a  frequency  of  1.  50  cps. 
The  f-k  spectrum  for  the  daytime  noise  sample  F  was  computed  by  a  different 
process  and  the  frequencies  (at  which  peaks  occured  in  the  spectrum)  shown 
in  Figure  11-20  are  0.3  cps  and  0.8  cps.  The  two  sets  of  spectra  are  similar, 
indicating  there  is  no  significant  change  in  the  noise  field  between  day  and 
night  or  from  November  to  March. 

At  0.3  cps  the  daytime  sample  has  a  high-velocity  noise  source 
situated  to  the  northeast  of  the  subarray.  All  other  noise  sources  are  down 
about  9  db  or  more  from  this  source.  The  0.25-cps  nighttime  noise  spectrum 
shows  these  same  features.  The  nighttime  noise  peak  is  broader  than  the 
daytime  peak,  but  the  nighttime  peak  at  a  higher  frequency  (0.  50  cps)  is 
broader  than  the  0.3 -cps  daytime  peak.  This  indicates  that  the  difference 
is  mainly  due  to  the  difference  in  frequencies  at  which  the  spectra  of  the 
two  samples  are  plotted.  The  difference  in  size  is  noticed  also  between  the 
0.75-cps  nighttime  spectrum  and  the  0.8-cps  daytime  spectrum.  The  pri¬ 
mary  difference  between  the  day  and  night  spectra  occurs  at  0.7  5  cps  where 
the  high  velocity  night  noise  is  essentially  isotropic  and  the  high  velocity 
day  noise  (at  0.8  cps)  is  from  the  west.  However,  the  2  km/sec  noise  sources 
to  the  northwest  and  west  are  present  in  both  samples. 

These  spectra  indicate  a  rather  stable  noise  field  between 
day  and  night  over  a  period  of  several  months. 
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Figure  11-18.  Frequency -Wavenumber  Spectra  of  Nighttime  Noise 

Sample  D;  f  =  0.  25  cps,  0.  50  cps,  0.  75  cps  and  1.  00  cps 
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Figure  11-19.  Frequency- Wavenumber 
Spectrum  of  Nighttime  Noise  Sample  D; 
f  =  1.50  cps 


Figure  11-20.  Frequency- Wavenumber 
Spectra  of  Daytime  Noise  Sample  F; 
f  =  0.  3 ,  0.8  cps 
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SECTION  III 

INSTRUMENT  RESPONSE  VARIATIONS 


A.  TIME  -STATIONARITY  AT  1  CPS 

1.  Purpose 

A  study  was  made  to  determine  if  the  LASA  seismometer 
amplitude  and  phase  responses  are  stable  as  functions  of  time.  Time- 
stability  of  seismometer  responses  is  a  necessary  requirement  for  effective 
on-line  multichannel  filtering. 

2.  Data  Used 

Nine  1  -cps  calibrations  covering  a  1 -month  period  (2  November 
1  December  1965)  were  chosen  for  subarrays  F-3,  F-4,  C-3,  E-2,  and  F-2. 
All  25  seismometers  from  each  subarray  were  demultiplexed,  and  a  few 
channels  from  each  subarray  for  each  calibration  were  plotted  to  check  the 
data.  The  signal -to-noise  (S/N)  ratios  were  large  (40  db  or  more)  on  all 
seismometers.  Table  III- la  lists  each  subarray's  calibration  length, 
which  was  constant  for  all  nine  calibrations. 


Table  Ill-la 


LENGTH  OF  CALIBRATION  FOR  EACH  SUBARRAY 

Subarray 

Length  of  Calibration 
(No.  of  cycles) 

F-3 

36 

F-4 

40 

C-3 

16 

E-2 

29 

F-2 

36 

In  addition,  a  1 -cps  calibration  performed  in  mid-March  was 
obtained  for  each  subarray.  Table  III  -lb  gives  the  date  for  each  subarray's 
calibration. 


The  seismometer  or  channel  numbers  used  in  the  text  refer  to  the 
sequence  in  which  the  channels  are  recorded  on  the  LASA  field  tape  for  a 
particular  subarray. 
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Table  Ill-lb 


l-CPS  CALIBRATION 

FOR  EACH  SUBARRAY 

Subarray 

Date 

F-3 

10  March  1966 

F-4 

14  March  1966 

C-i 

1 1  March  1966 

E-2 

1  5  March  1  966 

F-2 

9  March  1966 

3.  Method 


After  the  calibrations  were  demultiplexed,  gate  lengths  (Table 
Ill-la)  were  chosen  and  Fourier  transforms  were  computed  on  a  0.001-cps 
increment  for  all  seismometers  in  each  subarray  in  order  to  determine  the 
precise  calibration  frequency.  An  integral  number  of  cycles  were  chosen 
(within  visual  accuracy)  to  avoid  end  effects  in  the  transform.  The  amplitude 
and  phase  at  the  peak  (i.  e.  ,  calibration)  frequency  were  then  used  as  estimates 
of  the  seismometer  amplitude  and  phase  responses. 

To  remove  possible  variations  due  to  changes  in  the  signal 
input,  the  seismometers  were  referenced  to  channel  1.  (No  reference  trace 
was  available,  so  channel  1  was  arbitrarily  chosen.)  That  is,  the  amplitude 
and  phase  responses  relative  to  channel  1  were  studied: 


i  l  1 


Note  that  a  change  in  the  response  of  channel  1  (but  no  other  channel)  between 
calibrations  would  change  all  the  above  values  and  thus  give  artificial  variation, 
as  was  observed  on  two  occasions. 
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To  quantitatively  measure  the  stability  of  the  seismometer 
over  a  1  -month  period,  the  means  and  variances  of  the  amplitude  and  phase 
responses  were  computed  according  to 


n 


(where  n  =  9,  except  for  F-2,  where  n  =  8)  and  the  standard  deviation  (S)  was 
plotted  for  each  seismometer. 

Also,  the  amplitude  and  phase  responses  of  the  1  December  1965 
calibration  were  compared  with  those  of  mid-March  1966.  Percentage  changes 
in  amplitude  responses  and  absolute  changes  in  phase  responses  were  tabulated 
and  compared  with  corresponding  changes  over  the  1 -month  period. 

4.  Results 


a.  Stability  of  Calibration  Frequency 

Table  III -2  lists  the  calibration  frequencies  of  the  five  subarrays 
for  each  calibration  as  determined  from  the  Fourier  transforms.  Note  that 
F-3  has  a  significant  change  in  frequency  (.022  cps)  between  5  November  and 
16  November.  The  reason  for  the  change  is  not  known,  but  it  may  have  been 
due  to  replacing  the  F-3  oscillator  between  the  two  dates. 

All  other  calibrations  change  by  not  more  than  0.001  cps.  These 
small  changes  could  probably  be  explained  by  temperature  or  humidity  changes. 
Table  III -2  also  lists  the  calibration  frequencies  of  the  March  calibrations 
which  were  special  calibrations;  there  is  no  reason  that  these  should  agree 
with  the  others.  However,  since  the  values  will  be  relevant  in  later  discussions, 
they  are  listed  in  the  table. 

b.  Amplitude  and  Phase  Variations  Over  the  2  November  to  1  December 
Period 


In  general,  with  a  few  exceptions  the  amplitude  and  phase 
responses  showed  good  time -stability  over  the  1 -month  period.  An  attempt 


Note:  Refer  to  end  of  Section  III  for  Tables  III  -  2  through  III-27. 
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was  made  to  explain  all  anomalously  large  variations  by  relating  them  either 
to  low-quality  input  data  or  to  modifications  to  the  s eismometer  -amplifier 
system  as  reported  in  the  LASA  logs  (provided  by  Mr.  Royce  Brown),  or  by 
infering  a  possible  explanation  from  the  nature  of  the  variations. 

A  discussion  of  each  subarray  follows,  with  accompanying 
tables  and  figures  when  necessary. 

Subarray  F-3  (Tables  III  -  3  and  III -4)  was  unique  because  it  had 
a  shift  in  calibration  frequency  which  was  accompanied  by  large  variations  in 
amplitude  responses  and  small  variations  in  phase  responses.  The  amplitude 
variations  probably  were  partly  due  to  the  shift  in  calibration  frequency  and 
partly  due  to  a  lowering  of  the  channel  1  gain  on  12  November  (as  reported 
in  the  LASA  logs).  Because  channel  gains  are  not  altered  unless  they  exceed 
certain  specifications,  the  major  part  of  the  variations  probably  were  due  to 
the  gain  change  on  channel  1.  The  small  phase  variations  probably  were  due 
to  the  shift  in  calibration  frequency. 

There  are  four  anomalously  high  amplitude  ratio  values 
corresponding  to  channels  10,  14,  18,  and  25  for  the  16  November  calibration. 
The  playback  of  the  calibration  did  not  show  any  input  data  problems  on  these 
channels  and  verified  the  relative  amplitude  ratios  obtained  from  the  transforms. 
The  station  logs  did  not  have  any  reference  to  channel  adjustments,  although 
it  appears  these  four  rather  high-gain  channels  may  have  been  balanced.  Thus, 
no  definite  explanation  for  the  high  values  can  be  given. 

The  phase  difference  corresponding  to  channel  21  shows  a  large 
16°  change  between  5  November  and  16  November.  Again,  neither  the  playback 
nor  the  station  logs  provided  the  explanation  for  this  change. 

Table  III-5  gives  the  phase  means,  variances  and  standard 
deviations  for  all  nine  calibrations  and  the  amplitude  for  the  six  calibrations 
after  the  shift  in  calibration  frequency  and  gain  adjustment  to  channel  1. 

Figures  III  -  1  and  III  -  2  show  the  standard  deviations  in  graph  form.  The 
one  anomalously  high  value  on  channels  10,  14,  18  and  25  cause  their 
amplitude  standard  deviations  to  be  relatively  large.  For  other  channels, 
amplitude  response  standard  deviations  are  less  than  4  percent.  The 
phase  response  standard  deviations  are  less  than  1.  8  for  all  seismometers 
other  than  channel  21. 

Subarray  F-4  (Tables  III  -  6  and  III  -  7 )  generally  had  good  stability. 
The  channel  10  and  channel  24  amplitude  ratios  were  high  and  low,  respectively, 
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Figure  III- 1  •  Seismometer  Amplitude  Response  Variations  Over 
A  One  Month  Period.  Subarray  F -3 


Figure  HI-2.  Seismometer  Phase  Response  Variations  Over  A 
One  Month  Period.  Subarray  F -3 
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for  the  2  November  calibration.  Again,  no  explanation  was  available.  The 
phase  differences  for  channel  Z  had  relatively  low  values  on  Z  November  and 
21  November.  No  explanation  was  available  for  these  values.  The  channel  11 
phase  difference  had  much  less  stability  than  the  other  channels,  suggesting 
that  some  component  of  the  s eismometer  -amplifier  system  may  have  been 
deteriorating. 


Table  III  —  8  gives  the  means,  variances  and  standard  deviations 
for  all  24  amplitude  ratios  and  phase  differences.  Figures  III  -  3  and  III -4  show 
the  standard  deviations  in  graph  form.  The  amplitude  response  standard 
deviations  were  less  than  5  percent  on  all  channels  and  the  phase  response 
standard  deviations  were  less  than  0.4°  on  all  but  channels  2  and  11. 

Subarray  C-3  (Tables  III  -  9  and  III  -10)  had  excellent  stability  in 
amplitude  responses  except  for  a  large  abrupt  change  between  16  November  and 
20  November.  This  change  was  not  explainable  from  the  station  logs,  nor 
could  it  be  attributed  to  a  variation  in  the  channel  1  response.  (The  change 
varied  from  zero  on  channel  2  to  0.  29  on  channels  1  1  and  13.  )  Because  the 
change  varied  from  seismometer  to  seismometer  and  because  the  changes 
balanced  the  gains  at  1  cps  across  the  array,  the  most  likely  explanation  is 
that  the  s ei smometer -amplifier  systems  were  manually  adjusted. 

The  phase  differences  were  stable  except  for  the  channel  10 
response  which  changed  about  10  between  16  November  and  20  November. 

This  change  probably  was  related  to  the  amplitude  change  discussed  above. 

Also,  a  change  of  about  2  3  in  the  phase  response  of  channel  1  may  be  inferred 
between  3  November  and  5  November,  as  shown  by  the  2  change  in  all  24 
ratios . 


The  amplitude  means,  variances  and  standard  deviations  are 
not  presented  because  of  the  apparently  manually  caused  change.  The  phase 
data  (after  the  2°  phase  shift  on  channel  1  was  removed)  are  given  in  Table 
III  - 1  1  and  are  plotted  in  Figure  III -5.  Except  for  channel  10,  the  phase 
response  standard  deviations  were  less  than  0.4°. 

Subarray  E-2  (Tables  III -12  and  III  - 1  3)  had  stable  amplitude 
and  phase  responses,  except  for  eight  seismometers  which  had  abrupt  amplitude 
changes  between  23  November  and  24  November.  Again,  these  changes  appear 
to  be  due  to  manual  adjustments  of  these  comparatively  high-gain  systems. 
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Figure  III-3.  Seismometer  Amplitude  Response  Variations  Over 
A  One  Month  Period.  Subarray  F-4 


Figure  III-4.  Seismometer  Phase  Response  Variations  Over  A 
One  Month  Period.  Subarray  F-4 
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Figure  III- 5 .  Seismometer  Amplitude  Response  Variations  Over  A  One  Month 
Period.  Subarray  C-3 


Table  III  - 1 4  gives  the  means,  variances  and  standard  deviations 
for  both  the  amplitude  and  phase  responses,  except  for  the  eight  seismometers 
discussed  above,  and  Figures  III -6  and  III -7  present  the  data  in  graph  form. 
Amplitude  response  standard  deviations  were  less  than  3.  5  percent  and  the 
phase  response  deviations  were  less  than  0.3°. 

Subarray  F -2  (Tables  III  - 1 5  and  III  - 1  6)  amplitude  responses  were 
stable  except  for  the  last  calibration  on  2  December  which  was  different  for 
nearly  every  seismometer .  However,  since  the  diffe rences  were  not  regular,  it 
appears  that  they  were  natural  variations  due  to  some  physical  cause.  The 
phase  responses  were  stable. 

Table  III  - 1 7  gives  the  means,  variances  and  standard  deviations 
for  the  amplitude  and  phase  responses,  and  Figures  III -8  and  III- 9  present  in 
graph  form  the  standard  deviations  from  the  mean.  Amplitude  response 
standard  deviations  were  less  than  6  percent  and  phase  response  deviations 
were  less  than  0.6  . 

c.  Comparisons  of  1  December  and  Mid-March  Calibrations 

The  amplitude  and  phase  responses  for  the  1  December  1965 
calibration  are  compared  with  those  of  the  March  calibration  in  Tables  III  - 1 8 
through  III -22;  Table  III -Z 3  summarizes  the  variations.  About  two-thirds 
of  the  amplitude  response  variations  were  less  than  10  percent,  and  one -third 
of  the  phase  response  variations  were  less  than  2°.  For  the  1 -month  data, 
essentially  all  amplitude  and  phase  responses  were  within  these  limits 
(neglecting  the  abrupt  amplitude  response  changes  observed  on  subarrays  C -3 
and  E  -2) . 


Most  of  the  remaining  amplitude  response  variations  were 
between  10  and  20  percent,  but  ranged  up  to  a  maximum  of  65  percent  (channel 
10,  subarray  F-3).  Most  of  the  remaining  phase  response  variations  were 
between  2°  and  10°,  but  ranged  up  to  a  maximum  of  41°  (channel  6,  subarray 
F-3).  Because  of  the  substantial  number  of  large  variations,  instrument 
response  differences  were  not  included  in  the  five  multichannel  filters  designed 
for  on-line  application  (Section  V). 

As  shown  in  Table  III-2,  the  calibration  frequencies  were  about 
0.02  cps  different  for  the  two  sets  of  data.  This  difference  could  explain  a 
significant  part  of  the  phase  response  variations.  (Absolute  phase  responses, 
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AMP.  RATIOS 


Figure  III  -  6 .  Seismometer  Amplitude  Response  Variations  Over 
A  One  Month  Period.  Subarray  E-2 
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Figure  III  -  7 .  Seismometer  Phase  Response  Variations  Over  A 
One  Month  Period.  Subarray  E-2 
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AMP  RATIOS 


Figure  III-8.  Seismometer  Amplitude  Response  Variations  Over 
A  One  Month  Period.  Subarray  F -2 


10  U  14  16 

SEISMOMETER  NO. 


Figure  III- 9 •  Seismometer  Phase  Response  Variations  Over  A 
One  Month  Period.  Subarray  F-2 
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obtained  from  the  March  1966  calibrations  where  a  reference  channel  was 
included,  varied  quite  rapidly  in  the  neighborhood  of  1  cps.)  Otherwise, 
differences  in  phase  and  amplitude  response  probably  were  due  to  one  or  a 
combination  of  the  following: 

•  Natural  variations  due  to  environmental  changes 
and/or  instrument  drift 

•  Replacement  of  components 

•  Manual  adjustments 

Necessary  information  was  not  available  to  determine  which 
one  or  a  combination  of  these  three  was  the  major  cause  for  the  large  variations 
observed  in  both  amplitude  and  phase  responses.  Thus,  at  present,  no  con¬ 
clusion  can  be  made  about  the  stability  of  seismometer  responses  over  the 
longer  3-1/2  month  period. 

5.  Summary 

The  time -stability  at  1  cps  of  all  seismometers  from  each  of 
five  subarrays  was  studied  by  analyzing  their  amplitude  and  phase  responses 
relative  to  the  center  seismometers.  Nine  calibrations  recorded  over  a  1- 
month  period  (2  November  -  1  December  1965)  were  used  in  the  study.  The 
analysis  showed  that  the  seismometers  generally  had  good  time -stability  in 
both  amplitude  and  phase  responses  over  the  1 -month  period.  With  the  ex¬ 
ception  of  a  few  seismometers,  amplitude  response  standard  deviations 
(neglecting  the  abrupt  changes  discussed  below)  were  less  than  5  percent  and 
phase  response  standard  deviations  were  less  than  2°.  These  deviations  were 
small  enough  that  they  would  not  have  significantly  reduced  the  effectiveness 
of  an  on-line  system. 

In  addition  to  the  normal  variations  the  following  abrupt 
variations  were  observed: 

•  Variations  due  to  changes  in  the  seismometer- 
amplifier  system 

•  For  six  amplitude  responses  and  four  phase  responses, 
larger-than-normal  variations  which  had  no  immediate 
explanation 

•  For  subarrays  C-3  (all  seismometers)  and  E-2  (eight 
seismometers)  abrupt  changes  in  amplitude  responses 
(by  as  much  as  30  percent)  which  were  most  likely  due 
to  manual  adjustments  to  the  seismometer  system 
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The  relative  amplitude  and  phase  responses  of  the  1  December 
calibration  were  compared  with  responses  obtained  for  a  mid-March  1966 
calibration  on  all  five  subarrays.  Variations  significantly  larger  than  those 
seen  over  the  1 -month  period  were  observed.  Many  were  large  enough  to 
affect  the  performance  of  an  on-line  MCF  system.  However,  the  information 
necessary  to  determine  the  major  cause  of  the  variations  was  not  available. 
Thus,  at  present,  no  conclusion  can  be  made  about  the  stability  of  the 
seismometer  responses  over  the  3-1/2  month  period. 

B.  MULTIFREQUENCY  CALIBRATION  ANALYSIS 

1.  Purpose 

This  study  is  to  investigate  the  instrument  responses  of  125 
seismometers  as  determined  by  an  analysis  of  one  multifrequency  calibration 
of  each  seismometer.  The  results  of  this  investigation  were  used  to 
experimentally  determine  the  improvement  in  MCF  performance  to  be  obtained 
by  incorporating  these  measured  responses  into  the  MCF  design.  The  response 
variations  between  instruments  within  each  subarray  and  estimates  of  relative 
instrument  responses  are  presented. 

The  calibration  data  for  this  analysis  was  specially  recorded 
during  March  1 966  from  five  subarrays  (F-3,  F -4,  C-3,  E-2,  and  F -2)  at 
11  different  frequencies.  The  calibration  signal  length  at  each  frequency  was 
approximately  150  sec  and  was  recorded  at  a  high  S/N  ratio. 

2.  Method 


The  LASA  digital  tapes  containing  the  calibrations  were  de¬ 
multiplexed,  at  11  different  calibration  frequencies ,  to  obtain  26  traces  for 
each  of  the  five  subarrays.  The  26  traces  consisted  of  one  calibration 
reference  trace  and  the  recorded  output  of  each  25  seismometers  in  each 
subarray.  The  reference  trace  and  the  traces  of  two  selected  instruments 
were  Calcomp -plotted  for  each  subarray  to  determine  the  starting  time  and 
the  length  of  each  frequency  calibration  and  to  check  for  abnormalities  in  the 
data.  A  gate  length  of  an  integral  number  of  cycles  was  chosen  for  each 
frequency.  To  determine  the  exact  calibration  frequencies,  the  reference 
traces  were  Fourier -transformed  in  increments  of  0.002  cps  over  a  frequency 
band  centered  on  the  nominal  frequency.  Each  energy -density  spectrum  was 
plotted  and  the  frequency  at  maximum  amplitude  was  selected.  Results  of 
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this  analysis  are  shown  in  Table  III-24.  (Two  subarrays  were  not  calibrated 
at  certain  frequencies.)  The  digital  tape  containing  the  0.4-  and  0.6-cps  calibra¬ 
tions  for  F -3  was  damaged  and  the  6.0-cps  calibration  of  F-2  was  not  received. 

For  each  calibration,  after  exact  calibration  frequencies  had 
been  determined,  the  Fourier  transforms  were  calculated  for  each  channel  at 
the  corresponding  calibration  frequency.  Each  data  gate  began  at  a  zero 
crossing  and  contained  an  integral  number  of  cycles.  This  process  yielded 
the  phase  spectra  and  amplitude  spectra  for  all  of  the  125  channels  at  1 1 
frequencies /channel.  The  center  seismometer  of  each  subarray  was  chosen 
as  a  reference;  responses  were  not  referenced  to  the  calibration  reference 
channel  because  of  large  phase  differences  observed  between  the  reference 
channel  and  the  other  channels.  Since  the  cause  of  this  phase  shift  was  un¬ 
known  and  only  relative  responses  were  required,  it  was  decided  that  the  center 
seismometer  was  preferable  as  reference.  For  each  frequency,  differences 
in  the  phase  spectrum  for  the  center  seismometer  and  the  phase  spectra 
for  the  other  24  seismometers  were  formed.  Relative  amplitude  response 
was  measured  by  dividing  the  amplitude  spectra  of  the  remaining  24 
seismometers  by  the  amplitude  spectrum  of  the  center  seismometer.  By 
this  means,  relative  phase  and  amplitude  responses  were  obtained  for  all 
125  channels  at  1  1  frequencies /channel. 


Since  the  calibration  frequencies  did  not  contain  all  frequencies 
of  interest,  the  relative  responses  were  estimated  by  fitting  a  polynomial  to  the 
calibration  results  of  each  seismometer  and  evaluating  this  polynomial  at 
increments  of  0.25  cps.  The  method  of  least  squares  was  employed  to  obtain 
the  coefficients  of  a  polynomial  of  the  form 

f(X)  =  A  +  A,  X  +  AX2  +  .  .  .  +  A  Xn 
o  1  2  n 

An  investigation  was  made  to  determine  the  degree  polynomial 
best  suited  for  this  estimation.  The  set  of  measured  relative  responses  for 
one  channel  was  used  to  fit  polynomials  of  degree  3  through  9.  It  was  found 
that  a  polynomial  of  degree  5  fit  both  the  amplitude  and  phase  responses  very 
accurately  (less  than  2 -percent  error)  and  still  yielded  values  which  were  well- 
behaved  between  the  measured  points.  Then,  when  the  responses  for  a  com¬ 
plete  subarray  were  estimated  by  being  fitted  to  a  degree  5  polynomial,  the 
results  were  likewise  satisfactory.  Tables  III-25  and  III-26  contain  comparisons 
between  the  measured  relative  amplitude  and  phase  responses  and  the  poly¬ 
nomial-estimated  relative  amplitude  and  phase  responses  for  selected  channels. 
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The  relative  phase  and  amplitude  responses  of  channel  22,  subarray  F-3, 
represent  the  best  results.  The  relative  phase  response  of  channel  12, 
subarray  C-3,  and  the  relative  amplitude  response  of  both  channel  17, 
subarray  F-4,  and  channel  12,  subarray  C-3,  are  typical  results.  Poorest 
results  for  the  125  channels  are  the  relative  phase  response  of  seismometer  6, 
subarray  E-2,  and  seismometer  17,  subarray  F-4,  and  the  relative  amplitude 
response  of  seismometer  6,  subarray  E-2. 

The  relative  phase  responses  and  the  polynomials  which  were 
fit  to  these  responses  are  further  illustrated  in  Figures  III  - 1 0  and  III— 1 1.  The 
maximum  difference  between  measured  relative  phase  and  polynomial - 
estimated  phase  was  1°  for  most  of  the  125  seismometers,  which  is  sufficiently 
accurate  for  application.  In  Figure  III-10,  seismometer  18  of  subarray  E-2 
and  seismometer  6  of  subarray  F-2  are  representative  of  the  seismometers 
for  which  the  best  results  were  obtained;  seismometer  12  of  subarray  C-3  is 
representative  of  a  typical  result.  The  poorest  polynomial  fits  were  obtained 
for  six  seismometers  (6,  12,  17,  19,  23,  and  25)  of  subarray  F-4.  Figure 
III -11  presents  the  measured  phase  responses  and  the  polynomial -e stimated 
phase  responses  for  three  of  these  instruments.  The  difference  in  measured 
relative  phase  responses  and  polynomial -estimated  responses  was  less  than 
4^  for  all  six  seismometers. 

Figures  III  - 10  and  III  - 1 1  indicate  that  certain  of  the  polynomials 
have  large  excursions  between  4.  0  and  6.  0  cps,  that  are  not  justified  by 
the  measured  data.  However,  these  variations  were  not  considered  in  the 
signal  model  design  since  the  signal  model  was  band-limited  from  0.  5  to  4.  0 
cps. 


3.  Results 

Relative  phase  and  amplitude  responses  as  a  function  of  frequency 
are  represented  in  Figures  C-l  through  C-10,  Appendix  C.  Each  figure  con¬ 
tains  six  plots,  each  of  which  presents  the  relative  phase  or  amplitude  responses 
for  the  four  channels  from  one  of  the  six  arms  of  the  subarray.  No  relation 
was  observed  between  the  relative  responses  of  instruments  within  any  one 
particular  arm  of  the  subarray,  nor  between  relative  responses  of  channels 
equidistant  from  the  center  of  the  subarray. 

A  more  visual  illustration  of  the  response  variations  among  the 
25  seismometers  of  subarray  F-3  is  presented  in  Figures  C-ll  through  C-14, 
Appendix  C.  Most  relative  responses  are  well-behaved  for  frequencies  below 
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Figure  III-10.  Comparison  of  Relative  Phase  Responses  and  Polynomial  Fit 
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FREQUENCY  <  C  PS> 


Figure  III -11.  Comparison  of  Relative  Phase  Responses  and  Polynomial  Fit 
for  Subarray  F-4 
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3  cps.  The  exceptions  are  seismometer  16  of  subarray  C-3  (Figure  C-8)  and 
seismometer  6  of  subarray  F -3  (Figure  C-l).  Also,  the  amplitude  response 
for  channel  24  of  subarray  E-2  (Figure  C -4)  has  an  abrupt  deflection  at  1.0 
cps,  and  the  phase  response  for  channel  6  of  subarray  F-3  (Figure  C -6)  is  less 
than  that  for  the  other  instruments. 

The  means,  standard  deviations  and  extremes  of  the  relative 
phase  and  amplitude  responses  for  each  subarray  at  all  calibration  frequencies 
are  illustrated  in  Figures  III-12  through  III-16.  The  phase  response  of  channel 
1  6  of  subarray  C-3  and  the  phase  and  amplitude  responses  of  channel  6  of 
subarray  F -3  were  not  included  in  determining  these  values  (for  reasons 
previously  mentioned).  Of  interest  is  the  fact  that  the  minimum  standard 
deviation  from  the  mean  amplitude  response  occurs  at  1  cps  for  every  subarray. 
This  is  not  surprising  since  the  damped  resonant  frequency  of  each  seismometer 
is  approximately  1  cps. 

All  figures  indicate  substantial  variations  in  both  relative  phase 
and  amplitude  responses  within  each  subarray  for  frequencies  greater  than 
3  cps.  For  frequencies  below  3  cps  the  relative  responses  differ  noticeably 
from  seismometer  to  seismometer  as  Figures  C-l  through  C-14  illustrate. 
However,  the  standard  deviations  from  the  mean  are  less  than  those  at  the 
higher  frequencies* 

4.  Conclusions 

From  this  multifrequency  calibration  analysis,  it  has  been 
determined  that  there  are  strong  variations  in  phase  and  amplitude  responses 
among  the  seismometers  within  a  subarray.  These  variations  are  very 
pronounced  at  frequencies  greater  than  3  cps.  It  is  apparent  that  equalization 
of  these  response  differences  would  result  in  significant  improvement  in  MCF 
performance.  These  variations  also  would  also  seriously  distort  frequency- 
wavenumber  spectra  computed  without  proper  compensation  for  instrument 
response  differences. 
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and  Amplitude  Responses  for  Subarray  F 
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Figure  III-13.  An  Illustration  of  the  Mean,  Standard  Deviations  and  Extremes  of  Relative 
and  Amplitude  Responses  for  Subarray  F-4 
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Figure  III- 15 .  An  Illustration  of  the  Mean,  Standard  Deviations  and  Extremes  of  Relative  Phase 
and  Amplitude  Responses  for  Subarray  E-2 
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C.  DEVELOPMENT  OF  INSTRUMENT  RESPONSE  EQUALIZATION  FILTERS 
FOR  THE  B I  SUBARRAY 

Significant  differences  in  phase  and  amplitude  responses  exist 
among  the  instruments  within  the  LASA.  The  response  variation  from  channel 
to  channel  within  a  particular  subarray  must  be  properly  compensated  for  in 
the  design  of  multichannel  filters  for  that  subarray  in  order  to  achieve  maximum 
S/N  improvement.  Response  variations  within  the  array  also  make  invalid  the 
assumption  of  space -stationarity  which  is  necessary  to  the  computation  of 
f r equency -wavenumber  spectra. 

1.  Purpose 

The  experiment  described  herein  was  designed  to  investigate 
the  development  of  digital  equalization  filters  which  would  compensate  for 
differences  in  phase  and  amplitude  responses.  These  filters  were  designed 
on  the  assumption  that  the  response  differences  could  be  modeled  as  having 
been  produced  by  minimum -phase  filters.  Calibration  data  recorded  by  Texas 
Instruments  Digital  Field  System  (DFS)*  on  13  May  1965  during  the  engineering 
and  seismological  evaluation  of  subarray  Bl**  was  used  to  specify  filter 
response. 


The  filters  were  evaluated  on  the  basis  of  their  ability  to 
equalize  phase  response  to  that  of  an  arbitrarily  chosen  reference  seismometer. 
Whereas  phase  responses  varied  as  much  as  37  at  1  cps,  as  determined  from 
the  calibration,  application  of  the  equalization  filters  reduced  this  variation 
to  7°.  In  general,  considerable  reduction  in  the  phase  response  variation  from 

instrument  to  instrument  was  achieved  at  all  frequencies  for  the  22  channels 
for  which  equalization  filters  were  developed. 


Trade  of  Texas  Instruments  Incorporated. 

Texas  Instruments,  1965,  Large  Aperture  Seismic  Array:  Final 
Specifications  Rpt.  ,  Contract  AF33( 657)  - 1  3899,  August. 
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2.  Filter  Design 

All  channels  were  calibrated  on  a  0.  1-cps  frequency  increment 
from  0.  5  cps  to  1.  5  cps  and  at  2.  5,  3.  0,  4.  0,  5.  0,  6.  0,  and  7.  0  cps.  Two 
channels,  seismometers  10  and  33,  were  not  used  for  this  test  because  of 
large  deviations  from  the  average,  which  were  believed  to  be  due  to  other 
causes.  Calibration  data  were  recorded  digitally  on  a  0.  024- sec /sample 
basis.  All  calibrations  were  Fourier-transformed  and  amplitude  and  phase 
responses  plotted.  Channel  5,  which  corresponded  to  seismometer  44  for 
the  DFS  recording  format,  was  chosen  as  reference  channel  since  its  ampli¬ 
tude  response  seemed  to  be  representative  of  the  average  response.  Table 
III-27  presents  the  recording  format  used  for  this  data. 

The  ratio  of  the  amplitude  response  for  the  i^  channel  to  the 
amplitude  response  of  channel  5  was  formed  at  each  calibration  frequency. 
Then,  a  least-squares  fit  of  a  third-  to  seventh-degree  polynomial,  depending 
upon  which  appeared  to  fit  best  (but  usually  a  fifth-degree)  was  made  to  es¬ 
timate  the  amplitude  ratio  function  for  each  of  the  22  channels  to  be  equalized 
to  channel  5.  Evaluation  of  each  polynomial  on  0.  05-cps  increments  yielded 
the  required  estimate  of  the  amplitude  response  ratio  functions.  Figure  III-  1  7 
shows  in  decibels  the  relative  amplitude  responses  of  the  reference  channel 
and  channels  10  (35)  and  25  (53),  which  are  representative  of  the  remaining 
channel  responses. 


Table  III -27 


TRACE  IDENTIFICATION 

Channel 

Seismometer 

Channel 

Seismometer 

No. 

No. 

No. 

No. 

1 

71 

14 

26 

2 

31 

15 

33 

3 

10 

16 

73 

4 

24 

17 

51 

5 

44 

18 

21 

6 

84 

19 

64 

7 

82 

20 

62 

8 

42 

21 

25 

9 

22 

22 

55 

10 

35 

23 

66 

1 1 

75 

24 

23 

12 

86 

25 

53 

13 

46 
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Figure  III-17*  Relative  Amplitude  Responses  for  13  May  1965 
Calibration  Data 


The  requirements  for  equalization  filters  in  this  case  were 
that  each  be  minimum  phase  and  that  its  amplitude  response  be  the  reciprocal 
of  the  estimated  amplitude  response  ratio  function  for  that  channel.  The 
existing  program  for  designing  minimum- phase  whitening  filter  designs  in 
the  time  domain.  Input  to  the  program  is  an  autocorrelation  function  of  the 
time  series  to  be  whitened.  The  square  of  the  amplitude- re sponse  ratio 
function  was  Fourier-transformed  as  indicated  below  to  obtain  the  proper 
function  for  input  to  this  program: 

416 

}n<niT>  -  I  , 

m=-4  1  6 


A.(mAf) 

A^mAf) 


j2TTmAfmAt  Ar  _  _  oo 

e  Af  n=0, 1 f  •  *  • •  co 


whe  re 

.  (nAT)  “  ’’autocorrelation”  function  corresponding 
to  channel  i 

A.(mAf)  =  amplitude  response  of  channel  i 
Af  =0.  05  cps 
At  =  0.  072  sec 

The  22  autocorrelation  functions  were  input  to  the  filter  design 
program  and  a  29-point  equalization  filter  developed  for  each.  These  filters 
were  applied  to  the  channels  for  which  they  were  designed  for  two  sets  of 
calibrations.  One  calibration  was  that  used  in  estimating  relative  amplitude 
responses,  and  the  other  was  the  same-type  calibration  which  was  performed 
2  weeks  later.  The  degree  of  equalization  achieved  was  very  nearly  the  same, 
indicating  very  little  change  in  instrument  responses  over  the  2-week  period. 
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3.  Filter  Evaluation 


The  relative  phase  responses  for  the  22  channels  studied 
are  presented  in  Figure  III-18(a-d).  This  figure  series  also  shows  relative 
phase  following  application  of  equalization  filters.  Relative  phase  was 
measured  by  Four ier -transforming  the  calibration  data  after  equalization 
filtering.  Only  the  six  discrete  frequenc ie s  —  0.5,  0.7,  1.0,  1.5,  2.5, 
and  3.0  cps  —  were  examined,  although  a  continuous  curve  estimate  of  the 
relative  phase  was  sketched  through  the  measured  values.  Almost  all 
channels  show  considerable  improvement  in  similarity  of  phase  response 
to  that  for  channel  5.  Those  which  show  the  least  improvement  are  always 
the  channels  for  which  the  phase  response  closely  agrees  (within  ±5°)  with 
that  for  channel  5  before  equalization. 

The  relative  phase  following  application  of  equalization  filters 
is  shown  for  both  sets  of  calibration  data  in  Figure  III-19(a-d).  The  close 
agreement  at  all  frequencies  for  all  channels,  except  channel  21,  demon¬ 
strates  the  stability  of  instrument  response  over  the  2-week  interval.  The 
behavior  of  channel  21  is  unexplained  but  may  be  due  to  the  replacing  or 
adjusting  of  some  component  associated  with  this  channel  between  calibrations. 

Figure  III-20  presents  the  range  of  variation  in  relative  phase 
response  among  the  22  channels  for  the  calibration  set  used  in  filter  design 
and  the  same  data  after  equalization.  Variations  of  as  much  as  38°  at  1.  0 
cps  before  filtering  have  been  reduced  to  less  than  4°  by  equalization. 

While  marked  improvement  has  been  achieved  at  all  frequencies,  the  greatest 
percentage  of  improvement  has  been  obtained  at  1.  0  cps.  This  better  filter 
performance  at  1.  0  cps  may  indicate  that  the  assumption  that  the  differences 
in  seismometer-amplifier  responses  may  be  modeled  as  the  effect  of  a 
minimum-phase  filter  is  better  approximated  at  1.0  cps.  Figure  III- 2  1  shows 
the  same  comparisons  for  the  second  set  of  calibrations.  Very  similar  results 
were  obtained,  since  there  was  little  change  in  relative  phase  response  between 
calibrations . 


Ill -27 


science  services  division 
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Figure  III-18a.  Phase  Responses  for  May  13,  1965  Relative  to  Channel  5  Before 
(solid  curve)  and  After  (dashed  curve)  Filtering 
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Figure  III  - 1 8b. 


Phase  Responses  for  May  13,  1965  Relative  to  Channel  5  Before 
(solid  curve)  and  After  (dashed  curve)  Filtering 
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Figure  III-18c.  Phase  Responses  for  May  13,  1965  Relative  to  Channel  5  Before 
(solid  curve)  and  After  (dashed  curve)  Filtering 
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Figure  III  - 1 8d.  Phase  Responses  for  May  13,  1965  Relative  to  Channel  5  Before 
(solid  curve)  and  After  (dashed  curve)  Filtering 
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Figure  III-19a.  Comparison  of  Phase  Responses  for  13  May  (dashed  curve)  and 
27  May  1965  (solid  curve)  Relative  to  Channel  5  After  Filtering 
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Figure  III  - 1  9b.  Comparison  of  Phase  Responses  for  13  May  (dashed  curve)  and 
27  May  1965  (solid  curve)  Relative  to  Channel  5  After  Filtering 
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Figure  III-19c. 


Comparison  of  Phase  Responses  for  13  May  (dashed  curve)  and 
27  May  1965  (solid  curve)  Relative  to  Channel  5  After  Filtering 
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FREQUENCY  (CPS) 

Figure  III-l9d.  Comparison  of  Phase  Responses  for  13  May  (dashed  curve)  and 
27  May  1965  (solid  curve)  Relative  to  Channel  5  After  Filtering 
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Figure  III-20.  Seismometer  Phase  Response  Variations  Relative  to  Channel  5 
for  13  May  1965  Before  (upper  figure)  and  After  Filtering 
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Figure  III -21.  Seismometer  Phase  Response  Variations  Relative  to  Channel  5 
for  27  May  1965  Before  (upper  figure)  and  After  Filtering 
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Table  IU-3 
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SUBARRAY  F-3  MEANS,  VARIANCES  AND  STANDARD  DEVIATIONS  OF 
AMPLITUDE  RATIOS  AND  PHASE  DIFFERENCES 


> 

v 

Q 

"d 

4-> 

CO 


cMojcMcor—  r-'-Or-ocoincoococM^oo 

N  O'  O  H  rn  CT'CO['-'<tfMa'COfMincMCOO 


o 


H  CO 

o  co  in  ^ 


OO-H-HOOOO-HOOOOOOOOOOOOOOOO 


w 

CO 


Oh 


0) 

a 

c 

2 

u 

rt 

> 


NHfo^^inoHOHO(\j^OHin^o^ro^Hoo^^ 
O  ^  LD  — <  C'  in  h  h  N  COfM'^COO^HOfMOOCO^DOOOLO 

T^LDLna^^rvj^Hm^Hcn^oou^^-^-moo^Horo^HfMoo 
o  oo  o  n  ~-,fv]a''-H{\](\)Oco— <Ofvjoooor^co— <m— * 


OO^H^H  OOOOcO 


OOOOOOOOOOt^ 


o  o  o  o 


c 

aJ 

0) 

2 


o 

CO 

00 

OJ 

CO 

<M 

o 

o 

58 

CM 

vO 

33 

o 

o 

o 

o 

o 

00 

00 

o 

r- 

84 

CM 

O 

55 

23 

50 

43 

r- 

22 

37 

o 

CO 

00 

CO 

CO 

00 

00 

i-H 

o 

IT) 

vO 

00 

CM 

o 

CO 

iTi 

LTi 

CO 

h- 

CO 

> 

4) 

P 

xi 

4— > 

in 


OOrOO'HsOHin^DmO^COOOOO^oOHON(>fO\DN 

OOO^HO-HOrOr^OrOOLOi— <  OrOOssD 

oooooooooooooooooooooooo 

oooooooooooooooooooooooo 


w 

Q 

D 

a) 

00 

00 

in 

n 

CO 

oo 

00 

oo 

in 

o 

O' 

00 

in 

in 

o 

Tt< 

CO 

CO 

CO 

o 

in 

o 

00 

f-H 

CO 

i-H 

CM 

vD 

o 

co 

Is- 

r- 

CM 

o 

o 

o 

CO 

00 

CO 

in 

H 

d 

o 

o 

o 

o 

i-H 

o 

CO 

CO 

o 

co 

o 

CO 

I-H 

o 

o 

CM 

Tt< 

CM 

f-H 

o 

I-H 

o 

00 

►H 

aJ 

o 

o 

o 

o 

o 

o 

o 

in 

o 

o 

CO 

o 

o 

o 

CM 

o 

O 

o 

o 

»-H 

o 

CO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

& 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

2 

<; 

> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

d 

0) 

2 


N'tOrOrOrO^COCOfn^[s''^OO^OHOOO'XivDcO(T'sO 

oor-'tmm^oooooo^cocooo^^HsOcooooHH 

OOsOOOO,sOOCTsOOOOOOOCTsOCT''^hOC7vOO 

I- H  i-H  I-H  ^H  ^H  ^H  — H  ^H  ^H  i-H  I— H  i-H  i-H  H  i-H  i-H  i-H  i-H 


a) 

CO 


<Mco^ir»'Or-oooo^H 


Nro^iDvor^oooOHNro^m 

HHHHHHMpHfvJfNjfVjfVjfVjfVj 


IH -41 


science  services  division 


Table  IH-6 


IH -42 


science  services  division 


Table  HI-7 


HI-43 


science  services  division 


SUBARRAY  F-4  MEANS,  VARIANCES  AND  STANDARD  DEVIATIONS  OF 
AMPLITUDE  RATIOS  AND  PHASE  DIFFERENCES 


w 

co 

s 

& 


> 

Q 

d 

4-4 

CO 


0) 

a 

CO 

lO 

co 

xO 

OJ 

00 

CO 

xO 

m 

xO 

o 

o- 

m 

m 

OJ 

00 

OJ 

f-H 

o 

Is- 

a 

xO 

CO 

O) 

i-H 

in 

o 

o 

oo 

00 

CO 

OJ 

r- 

co 

xO 

xO 

xO 

CO 

^■H 

o 

CO 

Is- 

in 

2 

Is- 

o 

co 

OJ 

OJ 

OJ 

iO 

i-H 

OJ 

f-H 

in 

CO 

i-H 

o 

^■H 

CO 

(M 

CO 

o 

OJ 

CO 

o 

o 

o 

o 

o 

o 

o 

o 

^H 

o 

o 

»— H 

o 

o 

o 

o 

o 

o 

O 

o 

o 

«-H 

o 

> 

CO 

o 

o 

o 

o 

o 

o 

o 

o 

in 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

a 

V 

2 


Tf  xO  OLnxor-xOor—sOTtHincoT^cocoxoco  ^  \0 

00  O  CO  iH  <-H  fsj  Os)  H  «-H  o)  i-H  H  CO  OJ  *“H  i-H  Os)  ^  *"H  iH  iH  O  CO  Os) 

HoddddddoNod  o  doddddddddd 


r-r-^ONoa'^N^oohNcoco^m^sO 
in  h  o)  HHO'h-fOsOOooina'O'cocoNcoN 


OJ  CO  Is- 
oo  xO 


mrs-coo^Hcoin-HOojoo-HcocoLn'<t<LncocoLncoooxo 


W 

P 

O 

H 

•— i 
0< 


> 

<D 

Q 

d 

CO 


c 

rt 

2 


OOinoOsO^rOOO^CM^NO^CMnh^^Nh^  cO  h-  (M 

co-h-h^Ooj— <o)cocooj^hOOojOojoojo-hcooj— < 

oooooooooooooooooooooooo 

do  dddddddddddddo  dodo  dodo 


<u 

CO 

OJ 

CO 

xO 

in 

CO 

CO 

xO 

in 

in 

f-H 

CO 

cO 

OJ 

OJ 

xO 

f-H 

OJ 

Is- 

m 

cO 

u 

in 

CO 

OJ 

o- 

CO 

o 

00 

CO 

00 

CO 

CO 

00 

OJ 

m 

CO 

in 

CO 

in 

f-H 

xO 

CO 

Tf< 

c 

OJ 

CO 

o 

o 

in 

CO 

in 

in 

^H 

OJ 

o 

o 

xO 

o 

in 

o 

o 

OJ 

o 

Is- 

1-4 

rt 

»-H 

o 

o 

OJ 

o 

o 

o 

o 

^H 

^H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

»-H 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

rt 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

> 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

HHO'h-0'00'HNhNh-1<N^^ir)OHaNO^POOrO 

com— <oO'-OTfroTfOcoinot^'^cors-or^cocoojrs-'^co 
oo— icoooo-hoooooooococoo-hoooo' 

i-H  i-H  i-H  O  i-H  i-H  — H  i-H  i-H  i-H  i-H  i-H  i-H  i-H  CD  >-H  CD  CD  i-H  i-H  i-H  i-H  i-H  CD 


<D 

7} 


NO^m\Or^coa'OH 


NfO^in^hooa'OHrom^in 

I-H  I-H  I-H  I-H  I-H  I-H  I-H  I-H  Oj  Os)  Os)  Os)  Os)  Os) 


in- 44 


science  services  division 


m-45 


science  services  division 


o 

I 


,-Q 

cd 

H 


m 

W 

O 

£ 

W 

05 

W 


w 

co 

< 

X 

X 

PO 

I 

o 

< 

X 

X 

< 

(0 

X 

in 


O 

NO 

00 

CO 

m 

o 

m 

03 

co 

m 

03 

o 

00 

CO 

o 

*  1 

CO 

03 

m 

o 

i-H 

CO 

o 

0- 

O 

00 

CO 

00 

i-H 

i-H 

O 

03 

03 

i-H 

03 

03 

Dec 

CO 

00 

o 

03 

00 

03 

CO 

vO 

i— H 

i-H 

i-H 

vO 

CO 

03 

03 

o 

i-H 

i 

» 

1 

i 

1 

l 

i 

i— H 

i 

1 

1 

1 

i— H 

1 

i 

03 

vO 

CO 

O' 

00 

^H 

00 

03 

CO 

i— H 

o 

00 

03 

m 

m 

LD 

03 

CO 

lO 

03 

o 

CO 

t—t 

V O 

I—) 

m 

CO 

o 

00 

i— H 

03 

o 

03 

o 

in 

i-H 

co 

i— H 

r-H 

> 

0 

CO 

00 

i—t 

i—t 

o 

03 

00 

CO 

CO 

sO 

03 

i— H 

i— H 

so 

CO 

03 

03 

co 

o 

■ — « 

Z 

1 

1 

fH 

1 

i 

1 

1 

i-H 

i 

1 

1 

1 

1 

i 

i— H 

1 

i 

CO 

Oj 

CO 

00 

CO 

oo 

CO 

o 

m 

CO 

o 

m 

O 

N O 

00 

m 

00 

O 

vO 

CO 

m 

> 

03 

o 

i— H 

CO 

CO 

o 

00 

00 

CO 

CO 

00 

co 

CO 

03 

03 

00 

o 

o 

CO 

i— H 

00 

o 

o 

0 

Z 

CO 

00 

r-H 

o 

03 

co 

CO 

o 

i-H 

i— H 

i-H 

vO 

CO 

CO 

03 

03 

CO 

vO 

i-H 

1 

1 

i 

i 

1 

i-H 

1 

^H 

i 

1 

1 

1 

^H 

1 

i 

f— H 

OJ 

CO 

CO 

03 

in 

NO 

CO 

nO 

03 

CO 

03 

vO 

o 

CO 

CO 

Is- 

CO 

m 

00 

00 

o 

03 

> 

ot< 

03 

03 

03 

o 

CO 

03 

in 

00 

o 

i-H 

03 

co 

03 

o 

CO 

o 

CO 

i-H 

03 

0 

£ 

CO 

00 

_ 

i-H 

o 

03 

00 

co 

CO 

sO 

03 

03 

1— H 

sO 

CO 

CO 

03 

03 

CO 

-o 

^H 

i 

1 

i— H 

1 

i 

1 

i-H 

1 

i-H 

i 

1 

i 

i— H 

1 

i 

o 

OJ 

xO 

m 

Is- 

O' 

oo 

m 

CO 

o- 

m 

03 

CO 

nO 

CO 

— 1 

03 

lT) 

03 

o 

CO 

o 

CO 

co 

00 

O 

o 

03 

03 

CO 

03 

o 

i-H 

CO 

o 

i-H 

0 

CO 

00 

r- 

r— 1 

o 

03 

00 

CO 

CO 

sO 

03 

03 

•-H 

'O 

CO 

CO 

03 

03 

CO 

NO 

i— H 

£ 

1 

1 

r-H 

1 

i 

1 

1 

r-H 

1 

r-H 

1 

i 

I 

1 

1 

r-H 

i 

vO 

r— 1 

CO 

00 

CO 

m 

vO 

i— H 

CO 

03 

CO 

i-H 

O- 

CO 

Is- 

00 

m 

CO 

i— H 

O 

m 

m 

> 

vO 

o 

nO 

03 

i— H 

i— H 

m 

00 

i— H 

o- 

m 

CO 

i-H 

00 

i-H 

vO 

03 

i— H 

0 

00 

00 

o 

r- 

r-H 

o 

CO 

00 

co 

CO 

o- 

i— H 

^H 

o 

m 

00 

03 

03 

sO 

i— H 

1 

1 

(H 

1 

i 

i 

i— H 

1 

' 

i 

1 

1 

i 

t 

1 

i 

LD 

i— H 

o 

CO 

m 

00 

03 

co 

i-H 

03 

03 

Is- 

m 

CO 

o 

CO 

00 

03 

■O 

r- 

> 

CO 

00 

00 

03 

nO 

CO 

o 

in 

CO 

i-H 

0- 

m 

s-O 

o- 

o 

m 

o 

03 

njD 

m 

03 

co 

0 

►  y 

00 

00 

o 

Is- 

r-H 

o 

CO 

00 

CO 

CO 

Is- 

•^o 

i— H 

i-H 

o 

vO 

CO 

03 

r- 

03 

sO 

i-H 

Z* 

1 

1 

1 

' 

1 

i— H 

1 

i 

1 

1 

i— H 

1 

i 

CO 

00 

03 

03 

CO 

03 

r-H 

CO 

CO 

00 

00 

i-H 

in 

nO 

m 

CO 

o 

03 

03 

n£> 

co 

oo 

sO 

co 

> 

CO 

oo 

CO 

CO 

00 

CO 

CO 

CO 

00 

03 

o- 

00 

co 

i-H 

CO 

o 

CO 

i— H 

i-H 

o 

00 

0 

o 

o 

03 

CO 

o 

Z 

m 

o 

m 

i— H 

CO 

oo 

CO 

CO 

CO 

oo 

03 

r-H 

O 

co 

o 

i-H 

o 

£ 

i— H 

1 

1 

r-H 

1 

1 

1 

1 

i 

03 

i 

1 

i 

i 

i 

i 

i— H 

•— H 

1 

i 

i 

1 

03 

CO 

oo 

V.O 

CO 

o- 

03 

in 

o 

o 

03 

r- 

m 

03 

03 

CO 

r- 

co 

vO 

> 

i— H 

oo 

i— H 

LD 

^H 

CO 

N0 

oo 

•— H 

m 

CO 

00 

co 

03 

o 

•— H 

i— H 

vO 

r-H 

^H 

00 

0 

o 

03 

CO 

O 

i-H 

m 

O 

03 

r-H 

CO 

oo 

CO 

CO 

03 

oo 

03 

i— H 

o 

co 

O 

i— H 

o 

Z 

i— H 

1 

1 

r-H 

1 

1 

i 

1 

i 

03 

i 

i 

i 

i 

i 

i 

i-H 

r-H 

1 

i 

1 

CO 

03 

CO 

m 

vO 

00 

CO 

O 

»-H 

03 

CO 

m 

vO 

Is- 

oo 

co 

o 

03 

co 

m 

r-H 

i— H 

i-H 

i-H 

i-H 

^H 

—1 

i— H 

i-H 

i— H 

03 

03 

03 

03 

03 

03 

*03 

CO 

III-46 


science  services  division 


Table  III-ll 


SUBARRAY  C-3  MEANS,  VARIANCES  AND  STANDARD 
DEVIATIONS  FOR  PHASE  DIFFERENCES 

Seis 

Mean 

Variance 

Std.  Dev. 

2 

-9. 

17 

0.  0584 

0.  24 

3 

-8. 

44 

0.  1082 

0.  32 

4 

-11. 

07 

0.  0561 

0.  23 

5 

-7. 

35 

0.  0294 

0.  17 

6 

1. 

39 

0.  0425 

0.  21 

7 

0. 

08 

0.  0326 

0.  18 

8 

-3. 

04 

0.  0199 

0.  14 

9 

-18. 

18 

0.  0343 

0.  18 

10 

9. 

56 

21. 7735 

4.  66 

11 

3. 

54 

0.  1145 

0.  33 

12 

-7. 

30 

0. 0288 

0.  16 

13 

-6. 

77 

0.  0169 

0.  13 

14 

-1. 

85 

0. 0519 

0.  22 

15 

-1. 

87 

0.  0310 

0.  17 

16 

-1. 

06 

0.  0401 

0.  20 

17 

-6. 

21 

0.  0234 

0.  15 

18 

14. 

06 

0. 0372 

0.  19 

19 

-9- 

13 

0.  0237 

0.  15 

20 

2. 

10 

0. 0156 

0.  12 

21 

-7. 

51 

0.  0421 

0.  20 

22 

2. 

26 

0. 0375 

0.  19 

23 

3. 

97 

0.  0267 

0.  16 

24 

6. 

16 

0.  0187 

0.  13 

25 

1. 

17 

0.  0080 

0.  08 
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Table  III -23 


COMPARISON  OF  1  DEC  1965  AND  MID -MARCH  1966  CALIBRATION 

Amplitude 

Ph? 

'  se 

Percent  Change 

No.  of  Channels 

Diff. 

No.  of  Channels 

<5 

61 

<2 

37 

5-10 

18 

2-5 

39 

10-15 

28 

5-10 

35 

15-20 

10 

>  10 

9 

>  20 

3 

Table  III -24 

CALIBRATION  FREQUENCIES 


F-2 

F-3 

C-3 

F-4 

E-2 

0.  202 

0.  202 

0.  200 

0.  202 

0.  202 

0.  398 

— 

0.  398 

0.  400 

0.  398 

0.  600 

— 

0.  598 

0.  598 

0.  598 

0.  802 

0.  802 

0.  802 

0.  800 

0.  802 

1. 000 

1. 006 

1. 002 

1. 000 

0.  980 
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1. 206 
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1.424 

1. 418 

1. 430 

1. 420 
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2.  012 

2.  006 

2.  012 

3.  040 

3.  030 

3.  020 

3.  024 

3.  014 

4.  000 

4.  000 

4.  000 

4.  000 

3.  984 

_ 

6.  000 

6.  006 

6.  004 

6.  000 
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Table  III-25b 


POLYNOMIAL  -  ESTIMATED  RELATIVE  RESPONSES 


FREQUENCY 

C-3  (212) 
RELATIVE 

PHASE 

RESPONSE 

C-3  (212) 

RELATIVE 

AMPLITUDE 

RESPONSE 

F -4  (217) 

RELATIVE 

PHASE 

RESPONSE 

F-4  (217) 
RELATIVE 
AMPLITUDE 
RESPONSE 

0.  0 

10.  038 

1. 3081 

-0. 431 

0. 9267 

0.  25 

-1. 555 

1.  2059 

-1. 499 

0. 9254 

0.  50 

-7. 099 

1.  1309 

-1. 798 

0. 9173 

0.  75 

-8. 607 

1. 0775 

-1. 534 

0. 9034 

1. 00 

-7. 673 

1. 0409 

-0. 876 

0. 8852 

1.  25 

-5. 521 

1. 0169 

0.  038 

0. 8637 

1.  50 

-3. 038 

1. 0024 

1.  101 

0. 8400 

1.75 

-0.  817 

0. 9946 

2.  230 

0. 8151 

2.  00 

0.  796 

0.  9915 

3.  366 

0.  7900 

2.  25 

1.  669 

0. 9915 

4.  470 

0.  7654 

2.  50 

1.  839 

0. 9935 

5.  517 

0.  7419 

2.  75 

1. 467 

0. 9968 

6.  498 

0. 7200 

3.  00 

0.  804 

1.  0009 

7.  411 

0. 7002 

3.  25 

0.  145 

1.  0054 

8.  261 

0. 6826 

3.  50 

-0. 211 

1. 0102 

9.  059 

0. 6673 

3.  75 

-0. 003 

1.  0153 

9.  812 

0. 6542 

4.  00 

0.  952 

1. 0204 

10. 528 

0. 6430 

4.  25 

2.  719 

1.  0254 

11.  204 

0. 6331 

4.  50 

5.  200 

1. 0299 

11. 831 

0. 6240 

4.  75 

8.  093 

1.  0333 

12.  385 

0.  6145 

5.  00 

10. 858 

1.  0345 

12. 828 

0. 6036 

5.  25 

12. 664 

1.  0322 

13.  099 

0. 5900 

5.  50 

12. 358 

1.  0247 

13.  118 

0. 5720 

5.  75 

8.  419 

1. 0094 

12. 777 

0. 5477 

6.  00 

-1. 082 

0. 9835 

11.  940 

0. 5151 
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COMPARISON  OF  MEASURED  RELATIVE  RESPONSES 
AND  POLYNOMIAL  -  ESTIMATED  RELATIVE  RESPONSES 


U 

o 

2 

H 

tC 

w 

Dm 

Dm 


U  W 

£  2 
H  O 
<  CD 
J  w 

M  W 
Od  Od 


DJ 

W  w 
>  2 

H  ^ 

< 

J 
W 
cc 


w 

u 

2 

w 

OS 

u 

Dm 

Dm 


w 

> 

H 

1  u 
-J  w 

W  2 

£8 

SS 

< 

5  U 
a  w 
H  < 

w  a 

W  CD 


w 

> 

<  w 

3  g 

u  o 
01 2 
Q  w 
M  U 
OtJ  * 

|  s 

W  2 

2  Eu 


>< 

o 

2 

w 

D 

2 

cd 

Dm 


(» 

— «  r- 


oo  vO  rsj  qo 


oo 

o 


-h  rO 
nO  O  — 
O  — '  rvj 


O' 

fM 


aor-f'ooor- 
o  — .  O'  Tf  o 

vO  m  Tf  lD 

o  o  o  o  — < 


N  N  ro 


O  O 
in  cn 
r-  — ■ 


rvi  \0  oo  oo  no 
O  o  O  —  — i 


o  —  — <  -h 


Tf  fsj  rn  rvj 

^  <  m  rg 

o  o  o  o 


'T 

o 


co  co  m  m 


00  *—  fM  CO  co  co 


r-.  r- 

ro  . — •  r~-  0^ 
O'  T»<  o 


vO  (\1  Tj* 

o  —  — 

'TOO 


III-62 


science  services  division 


Table  III -2  6b 


POLYNOMIAL  -  ESTIMATED  RELATIVE  RESPONSES 

FREQUENCY 

F-3 

F-3 

E-2  (26) 

E-2  (26) 

RELATIVE 

RELATIVE 

RELATIVE 

RELATIVE 

PHASE 

AMPLITUDE 

PHASE 

AMPLITUDE 

RESPONSE 

RESPONSE 

RESPONSE 

RESPONSE 

0.  00 

0.  298 

1.  1105 

-8. 344 

0.  7753 

0.  25 

-4.  597 

1. 0827 

6.  553 

0.  7112 

0.  50 

-8. 264 

1. 0642 

15.  017 

0.  7645 

0.  75 

-11.  113 

1.  0542 

18. 617 

0. 8822 

1. 00 

-13.  466 

1. 0517 

18. 693 

1.  0241 

1.  25 

-15.  566 

1.  0560 

16.  367 

1.  1613 

1.  50 

-17.  587 

1. 0661 

12. 568 

1.  2749 

1.  75 

-19.  645 

1.  0812 

8.  039 

1. 3546 

2.  00 

-21. 804 

1.  1004 

3.  358 

1.  3974 

2.  25 

-24. 084 

1.  1233 

-1. 045 

1. 4061 

2.  50 

-26.  475 

1.  1491 

-4.  875 

1. 3880 

2.75 

-28.  944 

1.  1778 

-7. 955 

1. 3537 

3.  00 

-31. 441 

1. 2091 

-10. 208 

1.  3152 

3.  25 

-33.  912 

1.  2432 

-11. 644 

1.  2852 

3.  50 

-36. 309 

1. 2806 

-12. 340 

1.  2752 

3.  75 

-38.  593 

1. 3221 

-12.  427 

1. 2946 

4.  00 

-40.  750 

1. 3688 

-12. 074 

1. 3486 

4.  25 

-42.  797 

1. 4225 

-11. 469 

1.  4379 

4.  50 

-44.  789 

1. 4849 

-10. 806 

1. 5561 

4.  75 

-46. 835 

1. 5587 

-10. 268 

1. 6894 

5.  00 

-49.  098 

1 . 6469 

-10. 008 

1.  8146 

5.  25 

-51.  811 

1.  7529 

-10. 139 

1.  8979 

5.  50 

-55.  285 

1. 8810 

-10.710 

1. 8935 

5.  75 

-59.  914 

2.  0359 

-11. 698 

1. 7423 

6.  00 

-66. 190 

2. 2229 

-12. 986 

1.  3706 
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SECTION  IV 

DETERMINATION  OF  OPTIMUM  MCF  DESIGN  PARAMETERS 


Presented  in  this  section  are  results  of  the  evaluation  and 
comparison  of  13  multichannel  filters  (MCF)  for  the  purpose  of  determining 
an  optimum  combination  of  design  parameters.  The  series  of  filters  was 
developed  in  such  a  way  that  pairs  of  filters  could  be  compared,  with  only 
a  single  design  parameter  differing  between  the  two  filters. 

The  evaluation  was  based  on  four  types  of  measurements, 
with  more  or  fewer  measurements  for  certain  filters:  impulse  response, 
frequency  response,  frequency- wavenumber  (f-k)  response,  and  ambient 
seismic-noise  attenuation.  Four  of  the  filters  were  applied  to  recorded 
P-wave  signals  and  the  degree  of  signal  distortion  was  estimated. 

A  brief  description  of  each  filter  is  presented  first  and  a 
list  of  the  comparisons  are  presented  in  a  table.  Also,  a  general  des¬ 
cription  of  the  four  basic  measurements  is  presented,  leaving  exceptions 
to  be  mentioned  within  the  appropriate  section. 

A.  DESCRIPTION  OF  FILTERS 

Each  filter  used  in  determining  the  optimum  MCF  design 
parameters,  their  characteristics  and  comparisons  are  listed  in  Table  IV-1. 
All  filters  are  1.0  sec  long  unless  otherwise  noted.  The  filters  were  de¬ 
signed  for  the  LASA  standard  subarray  shown  in  Figure  IV-1.  The  F-3 
subarray  from  which  the  noise  samples  were  taken  is  oriented  so  that  the 
subarray  arm  containing  the  number  1  seismometer  of  each  ring  is  6°  west 
of  north. 


The  25 -channel  noise  samples  were  obtained  by  retaining 
each  of  the  25  seismometer  outputs  as  a  separate  channel.  The  8 -ring 
noise  samples  were  obtained  by  summing  the  noise  output  of  all  the  seis¬ 
mometers  on  a  given  ring  for  each  of  the  eight  rings  shown  in  Figure  IV-1. 
The  5 -ring  noise  samples  were  obtained  by  summing  rings  3  and  4,  rings 
5  and  6  and  rings  7  and  8  of  the  8 -ring  noise  sample.  (In  discussing  ring- 
model  filters,  the  terms  channel  and  ring  will  be  used  interchangeably.) 

The  filters  designed  using  this  ring-forming  technique, 
being  linear  filters,  are  such  that  the  set  of  filter  weights  listed  for  a 
particular  ring  may  be  applied  either  to  each  seismometer  output  on  the 
corresponding  ring  or  to  the  ring  sum. 
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SUMMARY  OF  MULTICHANNEL  FILTER  CHARACTERISTICS 
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IV -2 


0° 


■< - 7.0  Km - ► 


NOTE: 

Seismometer  No.  10  Is  300'  Deep 
Seismometers  21-65  Are  200'  Deep 


Figure  IV -  1 .  LASA  Standard  Subarray 
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Several  noise  samples  were  used  in  designing  the  filters. 
Noise  sample  B  is  a  330-sec  noise  sample  recorded  at  the  F-3  subarray 
on  24  November  1965  between  03:18  and  03:24  MST  (Figure  IV-2).  Noise 
sample  D  used  in  designing  most  of  the  short  filters  is  the  last  150  sec 
of  sample  B.  Noise  sample  A  used  only  in  evaluating  the  filters  is  the 
150  sec  immediately  preceding  sample  D.  Noise  sample  C  is  a  150-sec 
sample  recorded  at  F-3  about  25  hours  earlier  (23  November  1965  between 
02:19  and  02:22  MST). 


1 

B- 

- ► 

*1 

|  ^  U  "  r 

1  1 - 

\**— 150  SEC  150  SEC — ► 

-< - 330  SEC - ► 

r —  1  50  SEC— ^ 

K - 330  SEC- 


Noise  Sample  Recorded 
24  November  1  965  Between 
02 :1  8  and  02 :24  MST 


Noise  Sample  Recorded 
23  November  1965  Between 
01:1  9  and  01:22  MST 


\< - 150  SEC - H 

Noise  Sample  Recorded 
31  March  1966  Between 
13:58  and  14:01  MST 


Figure  IV-2.  Noise  Samples  Used  in  Designing  and  Evaluating 
Multichannel  Filters 
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1.  Design  Parameters  for  MCF-3 


MCF-3  is  a  25-channel  maximum-likelihood  filter,  i.  e.  , 
an  MCF  constrained  to  have  an  impulse  response  of  exactly  a  unit  impulse 
which  thus  permits  no  distortion  of  broadband  infinite  velocity  signals.  The 
effective  signal  model  (no  signal  model  as  such  is  required)  for  the  maxi¬ 
mum-likelihood  filter  design  is  an  infinite  velocity  model  with  an  infinite 
S/N  ratio. 


The  noise  used  in  designing  this  filter  was  noise  sample  D. 
The  mean-square -error  as  a  function  of  filter  length  generated  by  the  stan¬ 
dard  MCF  design  program  is  not  meaningful  in  this  case. 

MCF-3  was  used  for  comparison  with  MCF  -  5  and  MCF-6, 
both  ring-model  filters  for  which  the  signal  model  was  a  disk  in  wavenumber 
s  pace . 

2.  Design  Parameters  for  MCF-4 

MCF-4  is  an  8-ring  MCF  utilizing  an  infinite  velocity  signal 
model.  The  method  of  summing  the  channels  into  eight  rings  was  described 
in  previous  paragraphs.  The  signal  correlations  for  simulating  an  infinite 
velocity  set  are  all  identical.  The  noise  sample  D  autocorrelation  for 
channel  1  was  used  to  generate  this  set  and  scaled  to  give  a  S/N  ratio  of 
4  on  channel  1. 


This  filter  was  used  in  the  comparison  between  an  infinite 
velocity  signal  model  and  a  disk  signal  model  (MCF-5)  and  in  the  compari¬ 
son  between  the  use  of  whitened  noise  (MCF-8)  and  the  use  of  nonwhitened 
nois  e  (MCF  -4) . 

3.  Design  Parameters  for  MCF-5 

MCF-5  is  an  8-ring  disk  signal  model  MCF.  The  signal 
model  corresponds  to  constant  power  density  over  a  disk  in  wavenumber 
space  (k-space)  with  a  radius  corresponding  to  a  horizontal  component 
velocity  of  12  km/sec.  The  signal  spectrum  was  shaped  to  the  average 
spectrum  of  the  5-ring  noise  sample  D.  This  spectrum  differs  very  little 
from  the  average  spectrum  of  the  25-channel  noise  sample.  The  noise 
sample  used  in  designing  MCF-5  was  noise  sample  D,  and  the  S/N  ratio 
was  4. 


This  filter  was  used  for  comparison  with  the  maximum- 
likelihood  filter  (MCF-3)  and  was  used  for  comparison  with  an  MCF 
designed  with  an  infinite  velocity  signal  model  (MCF-4). 
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4.  Design  Parameters  for  MCF-6 

MCF-6  is  a  5-ring  disk  signal  model  MCF.  The  disk  signal 
model  and  the  noise  model  (measured  correlations  from  noise  sample  D) 
are  the  same  as  those  used  in  the  design  of  MCF-5  except  that  the  eight 
rings  were  summed  into  five  rings,  as  described  in  previous  paragraphs. 

Again,  the  signal  spectrum  was  shaped  to  give  a  constant  S/N  ratio  of  4. 

This  filter  was  used  in  several  comparisons:  between  the 
disk  signal  ring-model  filter  design  and  the  maximum -likelihood  filter 
design  (MCF-3);  between  using  measured -noise  statistics  for  the  noise 
model  (MCF-6)  and  using  a  theoretical  noise  model  (MCF-7);  between  an 
MCF  designed  with  a  constant  S/N  ratio  (MCF-6)  which  does  not  frequency- 
filter  and  an  MCF  designed  with  a  S/N  ratio  which  varies  with  frequency 
(MCF- 11),  thus  permitting  frequency-filtering  as  well  as  velocity -filte ring . 

5.  Design  Parameters  for  MCF-7 

MCF-7  is  a  5-ring  filter  designed  from  a  theoretical  noise  model. 
The  noise  model  corresponds  to  constant  power  density  within  an  annulus  in  k-space 
with  an  inne  r  radius  corresponding  to  a  velocity  of  8  km  /  sec  and  an  oute  r  radius 
corresponding  to  a  velocity  of  1  km /sec.  The  spectrum  of  the  signal  model  used 
for  the  design  of  MCF-5  and  MCF-6  was  used  also  for  the  signal  and  nois e  spectral 
models  for  MCF-7  .  The  S/N  ratio  was  4  at  all  frequencies . 

MCF-7  was  used  in  the  comparison  between  using  a  theoretical 
noise  model  and  using  measured -nois e  statistics  (MCF-6)  for  filter  design. 

6.  Design  Parameters  for  MCF-8 

MCF-8  is  an  8-ring  filter  designed  from  a  noise  model  of 
whitened  measured  correlations  from  noise  sample  D.  To  whiten  the 
noise  spectrum,  a  5-point  whitening  filter  was  designed  from  the  noise 
sample  D  autocorrelation  of  channel  1.  The  8-channel  correlations  were 
convolved  with  the  autocorrelation  of  the  whitening  filter.  This  resulted  in 
a  correlation  set  exactly  equivalent  to  the  set  which  would  have  been  obtained 
by  applying  the  whitening  filter  to  each  channel  and  then  computing  the  cor¬ 
relation  functions.  The  infinite  velocity  signal  model  consisted  of  a  set  of 
correlations,  each  identical  to  the  whitened-noise  autocorrelation  of  channel  1. 

In  the  whitening  process,  the  correlations  were  necessarily 
shortened.  Consequently,  MCF-8  is  only  0.8  sec  long. 

This  filter  was  used  in  the  comparison  between  using  whitened 
measured-noise  statistics  and  using  nonwhitened  measured-noise  statistics 
(MCF-6)  for  the  noise  model  in  MCF  design. 
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7.  Design  Parameters  for  MCF-9 


MCF-9  is  a  5-ring  filter  designed  from  an  infinite  velocity 
signal  model  and  a  measured  nighttime  noise  model.  The  noise  model 
was  obtained  by  calculating  the  5-ring  correlations  of  seven  separate  noise 
samples  recorded  at  the  F-3  subarray  between  01:00  and  05:35  MST  and 
between  5  November  1965  and  1  December  1965.  A  150-sec  portion  of  each 
recorded  sample  was  used  to  compute  the  correlations.  The  seven  corre¬ 
lation  sets  were  summed  to  obtain  the  average  nighttime  noise  correlations 
used  to  design  the  filter.  Each  of  the  signal  correlations  for  the  infinite 
velocity  model  were  identical  to  the  channel  1  noise  autocorrelation.  The 
S/N  ratio  was  4. 

MCF-9  was  one  of  the  filters  used  to  determine  the  noise 
reduction  gains  to  be  attained  by  applying  multichannel  filters  designed 
from  daytime  noise  statistics  (MCF-10)  during  the  day  and  multichannel 
filters  designed  from  nighttime  noise  statistics  (MCF-9)  during  the  night. 

8.  Design  Parameters  for  MCF-10 

MCF-10  is  a  5-ring  filter  designed  from  an  infinite  velocity 
signal  model  and  a  measured  daytime  noise  model.  The  four  daytime  noise 
samples  used  were  recorded  at  the  F-3  subarray  between  11:00  and  17:00 
MST  and  between  24  March  1966  and  5  April  1966.  A  5-channel  correlation 
set  was  computed  for  each  noise  sample,  and  the  four  sets  were  then  sum¬ 
med  to  produce  a  single  5-channel  noise  correlation  set.  The  noise  auto¬ 
correlation  for  channel  1  was  used  to  form  an  infinite  velocity  signal  cor¬ 
relation  set  with  S/N  ratio  equal  to  4. 

This  filter  was  used  in  the  comparison  between  the  effects 
of  using  daytime  or  nighttime  (MCF-9)  noise  statistics. 

9.  Design  Parameters  for  MCF-11  and  MCF-12 

MCF-11  and  MCF-12  were  designed  for  a  5-channel  ring- 
model  subarray  configuration.  The  noise  model  was  identical  for  each, 
consisting  of  a  correlation  set  computed  over  a  150-sec  interval  (noise 
sample  D)  recorded  at  subarray  F-3.  Both  filters  were  designed  with 
a  disk  signal  model  in  order  to  pass  the  k-  space  region  from  12  km/sec 
to  infinite  velocity.  However,  the  signal  model  for  MCF-12  incorporated 
variations  in  instrument  responses.  The  signal  model  for  each  filter  was 
band-limited  from  1.  0  to  2.  25  cps.  The  total  signal  energy  to  total  noise 
power  ratio  was  set  to  2  to  insure  a  S/N  ratio  greater  than  1  everywhere 
within  the  disk  in  k-space. 

MCF-11  and  MCF-12  were  directly  compared  to  evaluate 
the  decrease  in  signal  distortion  or  attenuation  to  be  achieved  by  correcting 
the  signal  model  to  include  instrument  response  variations.  MCF-11  was 
compared  also  with  MCF-6  to  determine  the  effectiveness  of  frequency¬ 
filtering  combined  with  velocity-filtering. 
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10.  Design  Parameters  for  MCF-13,  MCF-14  and  MCF-15 


MCF-13,  MCF-14  and  MCF-15  were  designed  for  a  5-channel 
ring  model  subarray  configuration.  The  noise  model  consisted  of  a  corre¬ 
lation  set  computed  over  a  330-sec  interval  (noise  sample  B)  of  ambient 
noise  recorded  at  the  F-3  subarray.  The  signal  model  consisted  of  a  point 
in  k-space  corresponding  to  an  infinite  velocity  signal.  A  constant  signal- 
to-organized-noise  ratio  of  4  was  used.  The  three  filters  differed  only  in 
length:  MCF-13  was  1-sec  long,  MCF-14  was  2  sec  long  and  MCF-15  was 
5.65-sec  long.  The  signal  model  was  not  corrected  for  variations  in  instru¬ 
ment  responses. 

MCF-13,  MCF-14  and  MCF-15  were  directly  compared  to 
examine  the  relationship  between  noise  rejection  capabilities  and  filter 
length. 

B.  CALCULATION  OF  FILTER  RESPONSES 

1.  Impulse  Response 

The  impulse  response  of  each  filter  was  calculated  on  a 
0.05-sec  increment  over  the  length  of  the  filter.  All  filters  were  designed 
at  this  sample  rate  from  data  recorded  at  the  same  rate. 

2.  Frequency  Response 

The  frequency  response  of  all  filters  was  calculated  from 
0.  00  cps  to  10.  00  cps  in  0.  25-cps  increments  except  for  the  5.  6-sec  filter 
(MCF-15)  and  the  1-sec  (MCF-13)  and  2-sec  (MCF-14)  filters  compared 
with  it.  These  three  responses  were  calculated  on  a  0.  05-cps  increment. 

3.  F requency - Wa venumbe r  Response 

The  f-k  response  of  each  filter  was  calculated  in  0.  25-cps 
increments  from  0.  25  cps  to  1.50  cps,  and  up  to  2.  00  cps  in  some  cases. 
However,  the  wavenumber  responses  are  presented  only  for  frequencies  of 
0.  25,  0.  50,  1.  00,  and  1.  50  cps.  The  responses  at  these  frequencies  are 
representative  in  character  and  correspond  to  the  signal  and  noise  areas 
of  primary  interest.  Velocity  rings  are  overlaid  in  red  on  each  plot.  In 
all  cases  the  velocity  at  the  edge  of  each  plot  corresponds  to  1 .  0  km/sec. 

The  contour  levels  used  vary  from  plot  to  plot  because  of 
great  differences  in  complexity  of  the  contour  plots.  In  general,  each  plot 
is  contoured  in  6-db  increments  to  a  level  which  seemed  appropriate  for 
that  plot.  Contours  around  areas  of  very  low  response  are  not  shown  if 
they  are  fairly  narrow  and  do  not  correspond  to  some  known  noise  peak. 

Also,  in  general,  the  contour  of  -3  db  is  shown  only  around  the  peak 
response  at  the  origin. 
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The  f-k  responses  for  the  ring-model  multichannel  filters 
have  hexagonal  or  triangular  symmetry.  However,  due  to  the  coarse 
wavenumber  increment  used,  differences  appear  in  the  contour  levels  of 
symmetrically  corresponding  features. 

The  f-k  response  at  1.50  cps  shows  an  area  slightly  larger 
than  a  unit  cell.  The  kx  and  ky  shown  have  a  maximum  value  of  1.50  cycles/ 
km  at  the  edge.  The  unit  cell  has  a  maximum  kx  of  1.  154  cycles /km  and  a 
maximum  k^r  of  1.  333  cycles/km. 

The  frequency-wavenumber  spectrum  of  noise  sample  D 
is  shown  in  Figure  II- 18  (Section  II)  for  comparison  with  the  f-k  responses 
of  the  filters. 

4.  Noise  Attenuation 

Most  of  the  multichannel  filters  were  applied  to  two  noise 
samples.  As  mentioned  in  paragraph  B,  noise  samples  A  and  D  were 
recorded  at  the  F-3  subarray  between  03:18  and  03:24  MST  on  24  November 
1965  (Figure  IV  -2) .  Noise  sample  C  was  recorded  between  02:19  and  02:22 
MST  the  previous  day  (23  November  1965)  and  noise  sample  E  was  recorded 
between  13:58  and  14:01  MST  on  31  March  1966. 

The  noise  attenuation  of  each  filter  was  calculated  as  the 
ratio  of  the  MCF  output  power  density  spectrum  to  the  input  noise  power 
density  spectrum.  Three  methods  of  representing  the  input  noise  were 
chosen. 


•  Summing  the  input  noise  traces, 
dividing  by  the  number  of  seismometers 
used  to  obtain  the  noise,  and  calculating 
the  power  density  spectrum  of  this  trace. 

This  power  density  spectrum  is  actually 
the  spectrum  of  the  noise  filtered  by  a 
straight  sum  process.  The  ratio  of  an 
MCF  output  power  density  spectrum  to 
the  straight-sum  power  density  spectrum 
will  be  referred  to  as  an  MCF/ST  SUM  ratio. 

•  Calculating  the  power  density  spectrum  of 
the  noise  on  the  center  channel  input.  The 
ratio  of  the  MCF  output  power  density  spec¬ 
trum  to  this  center  channel  input  power  density 
spectrum  will  be  referred  to  as  an  MCF/CC 
ratio. 

•  Calculating  the  average  of  the  power  density  spectra 

of  the  25  noise  traces  originally  used  to  obtain  the  ring- 
summed  noise  samples.  The  ratio  of  the  MCF  output 
power  density  spectrum  to  this  average  spectrum  will 
be  referred  to  as  an  MCF/ AV G  ratio. 
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When  it  is  necessary  to  refer  to  the  ratios  calculated  by  any  two  or  all  three 
of  these  methods,  they  may  be  termed  noise  attenuation  plots. 

After  calculating  the  MCF/AVG  noise  ratio  for  a  few  filters, 
it  was  decided  that  no  extra  information  was  gained  from  this  technique. 

The  ratios  are  shown  for  the  four  filters  for  which  they  were  computed. 

It  should  be  noted  that  all  the  filters  were  designed  to  estimate 
the  signal  at  the  center  channel. 

5.  Other  Responses 

In  some  cases  the  random  noise  response  is  shown  as  a 
function  of  frequency.  The  response  represents  the  MCF  output  power 
relative  to  input  power  when  the  input  consists  of  unit  power  random-noise 
on  all  channels . 

Also,  for  some  filters  the  mean-square -error  vs  frequency 
is  shown.  The  mean -square -error  is,  in  general,  a  very  reasonable  quantity 
by  which  to  rate  the  overall  performance  of  a  multichannel  processing  system. 
The  objective  for  using  such  a  system  is  to  combine  several  channel  outputs 
(each  containing  both  signal  and  noise)  in  a  linear  fashion  so  that  the  resultant 
single  output  is  the  best  possible  estimate  of  the  expected  signal.  The  lin¬ 
ear  combination  of  these  channel  outputs  is  accomplished  by  sets  of  filter 
weights  derived  from  an  iterative  solution  to  the  Wiener  optimum  multi¬ 
channel  theory.  These  weights  are  constrained  to  yield  the  minimum  aver¬ 
age  power  difference  (mean-square-error)  between  the  expected  output 
signal  and  its  estimate,  assuming  given  signal  and  noise  models  and  given 
length  for  each  set  of  weights.  A  measure  of  the  quality  of  this  estimation 
is  the  average  power  of  the  difference,  or  error,  between  the  expected  signal 
and  its  estimate.  This  error  may  be  thought  of  as  the  sum  of  the  error  in 
estimating  the  expected  signal  caused  by  the  presence  of  noise  and  that 
caused  by  distortion  in  the  estimated  signal.  Therefore,  if  signal  distortion 
is  negligible,  the  reduction  in  the  mean-square-error  is  a  direct  measure 
of  the  reduction  in  total  noise  power  achieved. 

C.  COMPARISON  OF  FILTERS 

Comparisons  between  the  filters  described  in  the  previous 
paragraphs  will  now  be  discussed.  In  most  cases  the  only  difference  between 
the  two  filters  being  compared  is  a  single  design  parameter.  The  filter 
design  parameters  and  a  list  of  comparisons  are  summarized  in  Table  IV  -  1 . 

1.  Comparison  of  Maximum-Likelihood  (MCF -3),  8-Ring  (MCF-5) 
and  5-Ring  (MCF-6)  Models 

The  maximum-likelihood  filter  (MCF-3)  was  compared  with 
both  an  8  -  ring  filter  (MCF-5)  and  a  5-ring  filter  (MCF-6)  which  were  de¬ 
signed  from  the  same  measured  noise  (noise  sample  D)  used  to  design  MCF-3. 
MCF-5  and  MCF-6  were  designed  from  the  same  signal  model  which  was  a 
disk  in  wavenumber  space  with  a  radius  co r r esponding  to  12  km/sec. 
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As  mentioned  before,  MCF-3  was  constrained  to  have  an 
impulse  response  of  exactly  a  unit  impulse.  An  indication  of  the  constraint 
accuracy  is  given  by  the  figures  presenting  impulse  response  and  frequency 
response  measurement.  The  impulse  response  (Figure  IV-6)  has  a  zero 
time  value  of  1.  0000000  and  all  other  points  are  less  than  4.82x10"®  in 
magnitude.  The  phase  response  (Figure  IV-6)  is  0.00°  ±  0.005°  at  every 
calculated  frequency  (0.00  cps  to  10.00  cps  in  0.25-cps  increments),  and 
the  amplitude  response  varies  between  0.99999972  and  1.0000001. 

The  mean-square-error  as  a  function  of  filter  length  generated 
by  the  standard  MCF  design  program  is  not  meaningful  in  this  case. 

The  frequency-wavenumber  (f-k)  response  of  MCF-3  is  shown 
in  Figure  IV-16.  Examination  of  the  f-k  spectra  of  noise  sample  D  (Figure  II- 
18,  Section  II)  from  which  this  filter  was  designed  shows  a  source  of  high- 
velocity  noise  to  the  northeast.  This  peak  accounts  for  the  rapid  falloff 
of  the  filter  response  toward  the  northeast. 

At  0.  25  cps  the  main  response  lobe  centered  at  the  origin 
in  k-space  is  very  broad,  being  down  less  than  3  db  at  3  km/sec.  At  0.  50  cps 
the  main  response  lobe  is  much  narrower  than  at  0.25  cps,  but  the  response 
has  no  deep  holes  except  at  a  few  spots  outside  the  2-km/sec  ring.  At  1.  00 
cps  the  main  peak  is  still  narrower  and  the  nearby  reject  zone  is  much  deeper 
than  at  0.  50  cps. 

Several  sharp  reject  areas  which  appear  in  the  f-k  response 
deserve  comment.  Outside  of  the  main  noise  peak,  the  noise  power  is  so 
low  that  the  f-k  spectra  plots  (Figure  11-18  ,  Section  II)  of  noise  sample  D 
are  dominated  by  the  spectral  window  which  appears  on  the  plots  as  peaks 
of  about  -12  db.  Therefore,  it  is  possible  that  there  are  other  weak, 
though  coherent,  noise  sources  which  do  not  appear  on  the  f-k  spectra  plots 
and  that  MCF-3  has  designed  to  reject  these  noise  components.  However, 
most  of  the  holes  shown  on  the  MCF-3  response  remain  at  constant  k  as  the 
frequency  increases,  indicating  a  possible  array  geometry  effect.  Since  the 
transform  ^of  a  discrete  function  is  the  sum  of  a  finite  and  usually  small  num¬ 
ber  of  periodic  functions  which  are  harmonics  of  one  function,  the  transform 
may  be  expected  to  have  an  oscillating  character  which  is  basically  a  function 
of  the  sample  interval  of  the  original  variable.  The  severity  of  the  oscillation, 
the  peak-to-trough  amplitude  variation,  generally  decreases  as  the  number  of 
samples  increases.  For  the  f-k  response,  the  number  of  samples  is  small 
so  that  the  oscillatory  character  might  be  expected  to  be  noticeable  for  the 
transform  of  some  functions,  i.e.,  for  the  f-k  response  of  some  filters  at 
certain  frequencies.  Furthermore,  the  original  samples  are  unevenly  located 
in  x-y  space,  making  almost  impossible  calculating  by  hand  the  variation 
to  be  expected  in  the  f-k  response  due  to  this  sampling  problem.  It  should 
be  noted  that  this  oscillatory  response  is  not  erroneous,  but  is  the  actual 
transform  of  the  original  discrete  function  in  x-y  space. 
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The  wavenumber  spectrum  for  1.50  cps  shows  the  same 
trends  noted  at  1.  00  cps,  viz.  ,  sharper  main  lobe,  and  a  few  deep  holes 
at  the  same  wavenumbers  as  before. 

The  noise  attenuation  plots  for  MCF-3  are  shown  in  Figure 
IV-24  as  MCF/ST  SUM  ratio  and  in  Figure  IV-32  as  MCF/CC  ratio.  At 
frequencies  below  0.5  cps,  MCF-3  achieves  more  noise  attenuation  than  the 
straight  sum  process,  with  a  maximum  difference  between  the  two  processes 
of  10.  5  db  at  0.  16  cps  for  noise  sample  D  and  10.  0  at  0.  16  cps  for  noise  sam¬ 
ple  A.  At  high  frequencies  the  rejection  of  MCF-3  is  about  10  db  less  than 
that  of  a  straight  sum. 

MCF-5  is  the  constant  S/N  ratio  ring-model  filter  designed 
from  a  disk  signal  model  and  the  same  noise  sample  used  in  the  design 
of  MCF-3.  The  impulse  and  frequency  responses  of  MCF-5  are  shown  in 
Figure  IV-8,  and  the  mean-square  -  error  reduction  is  shown  in  Figure  IV- 
14.  The  impulse  response  has  a  peak  of  0.96  amplitude  at  zero  time, 
but  also  has  a  doublet  of  0.  1  amplitude  on  the  right  end.  This  doublet 
results  from  the  periodic  component  of  increasing  amplitude  in  the  fre¬ 
quency  response. 

The  f-k  response  is  shown  in  Figure  IV-18.  At  0.50  cps 
the  response  drops  sharply  at  velocities  just  below  12  km/sec,  having 
a  deep  reject  zone  around  4  km/sec.  This  sharp  drop  is  necessary 
because  of  the  high-velocity  noise  source.  Since  the  noise  source  is 
slightly  broader  in  k  -space  than  is  the  signal  model,  the  MCF  is  able  to 
reduce  substantially  the  total  noise  power  by  a  sharp  response  peak  at 
high  velocity.  The  response  behaves  in  a  similar  manner  at  0.  25  cps  and 
1.  00  cps.  At  1.  50  cps  the  main  noise  source  is  narrower  than  the  signal 
model,  and  the  filter  is  unable  to  reject  much  of  the  main  noise  source. 

MCF- 6  is  a  5 -ring  MCF  with  constant  S/N  ratio.  It  was  designed 
from  the  same  disk  signal  model  and  measured  noise  as  that  of  MCF-5. 

The  impulse  and  frequency  responses  are  shown  in  Figure  IV-9  and  the 
reduction  in  mean  -  square  -  error  is  shown  in  Figure  IV- 14.  There  is  no 
significant  difference  between  these  responses  and  those  of  MCF-5. 

The  noise  attenuation  plots  for  both  MCF-5  (Figures  IV-26, 

IV -34  and  IV-41)  and  MCF-6  (Figures  IV-27,  IV-35  and  IV-42)  have  several 
characteristics  in  common  with  the  plots  for  most  of  the  other  filters.  Some 
of  the  general  characteristics  will  be  described  now. 

An  examination  of  the  MCF/CC  ratios  indicates  that  at  some 
frequencies,  usually  the  high  frequencies,  the  MCF  noise  power  output  is 
greater  than  the  center  channel  noise  power  input.  However,  a  comparison 
of  the  MCF/CC  ratios  and  the  MCF/AVG  ratios  indicates  that  above  2.5  cps 
the  noise  power  on  the  center  channel  is  6  or  7  db  below  the  average  noise 
power  in  sample  D  and  is  8  to  10  db  below  the  average  noise  power  in  sam¬ 
ple  A. 
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This  difference  between  the  center  channel  noise  power  and  the 
average  input  noise  power  results  in  multichannel  filters  whose  noise  power 
output  at  high  frequencies  is  greater  than  the  center  channel  noise  power  in¬ 
put  but  is  less  than  the  average  noise  power  input.  The  difference  between 
the  spectrum  of  the  center  seismometer  (at  a  500-ft  depth)  and  the  maxi¬ 
mum  and  minimum  of  the  spectra  from  the  six  seismometers  on  ring  2 
(at  a  200-ft  depth)  is  shown  in  Figure  IV-3.  The  difference  between  the 
spectrum  of  the  center  channel  input  and  that  of  the  output  of  the  straight  sum 
process  is  shown  in  Figure  IV-4. 

Figures  IV-3  and  IV-4  indicate  that  the  bulk  of  the  noise  is 
below  1.  0  cps.  Above  1.0  cps  the  power  density  is  more  than  20  db  below 
the  peak  power.  Because  so  much  of  the  noise  power  is  concentrated  at 
these  low  frequencies,  filters  designed  in  the  time  domain  may  be  expected 
to  concentrate  on  this  frequency  band  and  neglect  the  response  at  the  higher 
frequencies.  In  general,  the  MCF/AVG  ratios  show  that  at  low  frequencies 
the  filters  have  high  attenuation.  The  attenuation  decreases  with  increasing 
frequency  until  1.  5  cps  to  2.  0  cps;  at  the  higher  frequencies  the  noise  be¬ 
comes  more  random  and  the  noise  attenuation  ratio  decreases  toward  the 
straight  sum  limit  of  -14  db.  The  plots  actually  remain  3  or  4  db  above  this 
limit;  this  difference  seems  to  be  caused  by  a  change  in  the  S/N  ratio  as 
described  below. 

The  disk  signal  model  filters  were  designed  by  shaping  the 
signal  model  to  the  spectra  average  of  noise  sample  D.  The  ratio  of  the 
total  signal  power  to  total  noise  power  on  channel  1  was  set  to  4.  The 
great  bulk  of  noise  power  is  below  1.  5  cps;  at  these  low  frequencies  the 
channel  1  noise  power  spectrum  is  virtually  identical  to  the  average  noise 
power  spectrum.  However,  at  high  frequencies  the  channel  1  noise  power 
spectrum  is  about  8  to  10  db  below  the  average  power  spectrum.  Setting 
the  S/N  ratio  to  4  adds  another  6  db  to  the  signal  spectrum,  giving  a  total 
of  14  to  16  db,  or  a  S/N  ratio  of  about  32.  This  variation  in  S/N  ratio 
allows  the  filter  some  frequency-filtering.  Since  the  S/N  ratio  increases 
with  increasing  frequency,  the  effect  will  be  to  allow  the  MCF  to  pass  more 
high-frequency  information  than  would  normally  be  passed  by  a  straight 
sum  process. 


Another  characteristic  in  some  of  the  plots  of  relative  noise 
attenuation  is  the  periodicity  at  frequencies  above  1.0  or  2.  0  cps.  This 
periodicity  may  be  attributed  to  spectral  window  effects  when  generating 
the  signal  spectrum  from  the  noise  autocorrelation  for  which  the  low- 
frequency  peak  dominates  the  spectrum.  Thus,  at  high  frequencies,  in 
addition  to  the  gradual  increase  in  S/N  ratio,  there  is  a  periodic  variation 
which  produces  more  frequency-filtering  that  results  in  periodicity  as 
observed  in  the  noise  attenuation  plots. 
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Figure  IV-3.  Power  Density  Spectra  for  the  Center  Seismometer  and  the 
Maximum  and  Minimum  of  the  Six  Seismometers  on  Ring  2 
for  a  Nighttime  Noise  Sample 
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Figure  IV-4.  Straight  Sum  and  Center  Channel  Power  Density  Spectra 
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The  MCF/ST  SUM  ratios  show  re  suits  similar  to  those  of  the 
MCF/AVG  ratios,  i.  e.  ,  for  most  of  the  filters  the  ratios  are  about  +4  to 
+  6  db  from  2.  5  cps  to  5.  0  cps,  indicating  that  at  these  frequencies  the  straight 
sum  process  has  4  to  6  db  greater  noise  rejection  than  the  multichannel  filters. 

Many  of  the  attenuation  curves  slope  sharply  downward  at 
5.  0  cps.  In  every  case  the  curves  were  calculated  to  higher  frequencies  and 
show  only  a  sharp  dip  at  5.  0  cps  with  the  general  level  of  the  curve  above 
5.  0  cps  returning  to  the  level  of  the  curve  below  5.  0  cps.  Examination  of 
the  center  channel  spectra  and  of  the  straight-sum  trace  spectra  for  samples 
D  and  A  reveals  large  peaks  at  5.  0  cps.  This  noise  component  apparently 
was  coherent  across  the  array  and  some  rejection  was  possible  on  the  basis 
of  wavenumber. 

Similarly,  the  large  dip  in  the  noise  sample  D  spectrum 
between  1.5  cps  and  2.0  cps  results  in  less  attenuation  than  at  adjacent 
frequencies  as  evidenced  by  most  of  the  noise  attenuation  plots. 

As  previously  mentioned,  these  characteristics  of  the  noise 
attenuation  plots  are  common  to  most  of  the  filters. 

The  results  of  analyzing  the  f-k  responses  of  MCF-3 
(Figure  IV  -16)  and  MCF-5  (Figure  IV- 18)  indicate  that  at  the  two  lower 
frequencies  (0.25  cps  and  0.50  cps)  the  response  to  high  velocities  falls  off 
much  more  rapidly  with  decreasing  velocity  and  to  a  lower  value  in  the  re¬ 
jection  zone  for  MCF-5  than  for  MCF-3.  However,  at  low  velocities  (below 
2  km/sec)  MCF-3  remains  lower  than  MCF-5  which  comes  back  up  to  -6  db 
at  0.  25  cps  and  to  +6  db  at  0.  50  cps.  At  1.  00  cps,  MCF-3  continues  to  fall 
off  rapidly,  dropping  to  -6  db  over  about  90  percent  of  the  azimuth  between 
6  and  7  km  /sec;  in  some  directions  MCF  -3  drops  to  -18  dbat  the  same  veloc¬ 
ities  MCF-5  only  drops  to  -6  db,  and  over  less  than  half  the  azimuth  at 
6  km/sec;  at  3  km/sec  it  is  back  to  0  db  and,  except  for  a  -6  db  rejection 
ring  around  1.  6  km/ sec,  the  response  is  high  (above  0  db)  over  most  of  the 
region  at  velocities  below  3  km/sec.  At  1.50  cps,  integration  of  the  two 
responses  over  the  area  inside  the  2-km/sec  circle  would  appear  to  show 
much  less  power  in  the  MCF-3  response  than  in  the  MCF-5  response.  The 
MCF-3  response  has  a  sharp  peak  at  infinite  velocity  whereas  MCF-5  is 
flat,  within  ±3  db  at  velocities  above  approximately  4  km/sec. 

The  MCF-6  response  (Figure  IV  -19)  has  a  sharper,  more 
isotropic  peak  at  high  velocity  and  stays  lower  at  low  velocity  than  does  the 
response  of  either  MCF-3  or  MCF -5  at  0.25  cps.  At  0.50  cps,  the  response 
of  MCF-6  falls  off  more  rapidly  and  drops  to  a  lower  value  than  the  MCF-3 
response;  however,  it  is  about  the  same  as  the  MCF-5  response  in  this 
respect.  MCF-6  does  not  have  peaks  at  low  velocity  which  rise  as  high  as 
those  of  MCF-5,  but  MCF-6  seems  to  have  a  larger  area  in  k-space  above 
0  db.  Again,  MCF-6  is  more  isotropic  than  MCF-5  or  MCF-3.  At  1.  00  cps, 
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the  response  of  MCF-6  has  a  more  uniform  rejection  ring  around  4  km/sec 
than  does  that  of  MCF-5.  Furthermore,  the  MCF-6  response  does  not  come 
back  up  as  much  at  2  km/ sec  as  does  that  of  MCF-5.  At  1. 6  km/ sec  MCF-5 
has  a  rejection  ring  which  MCF-6  does  not  have,  though  MCF-6  appears  to 
be  trying  to  reject  around  this  velocity.  At  this  frequency,  MCF-3  falls  off 
much  more  sharply  from  high  velocity  to  a  much  lower  value  in  most  directions 
than  does  MCF-6;  but  to  the  south  and  northeast,  MCF-3  has  very  little  rejec¬ 
tion  at  velocities  above  2  km/sec.  Below  2  km/sec  the  MCF-3  response  is 
above  0  db  over  a  large  area  of  k-space.  At  1.  50  cps,  MCF-6  attempts  to 
form  a  rejection  ring  around  5  km/sec  whereas  MCF-5  does  very  little  at 
this  velocity.  Both  filter  responses  rise  between  4  and  3  km/sec;  however, 
below  2  km  /  sec,  MCF-6  seems  to  hold  its  response  down  over  a  much  larger 
area  in  k-space  than  does  MCF-5.  MCF-3  has  a  good  rejection  zone  between 
4  and  8  km/sec  in  most  directions,  though  again  it  is  far  from  isotropic 
(having  little  rejection  in  the  north  at  any  velocity).  MCF-3  also  tends  to 
fall  off  rapidly  in  the  signal  area,  being  down  more  than  6  db  in  some  places 
at  velocities  greater  than  12  km/sec. 

In  general,  the  f-k  response  of  MCF-6  seems  to  be  as  good 
as,  or  better  than,  MCF-5  at  the  four  frequencies  examined.  The  response 
at  high  velocity  usually  falls  off  more  rapidly  and  the  response  at  low  velocity 
is  generally  lower  for  MCF-6  than  for  MCF-5.  Furthermore,  MCF-6  is 
considerably  more  isotropic  than  MCF-5. 

MCF-3  and  MCF-6  are  not  as  easily  compared.  MCF-6 
generally  has  sharper  falloff  from  the  signal  area  and  a  deeper  rejection 
zone  adjacent  to  the  signal  area  at  the  two  lowest  frequencies.  Over  all 
k-space,  MCF-6  obviously  has  greater  rejection  at  0.  25  cps  and  appears 
to  have  greater  overall  rejection  at  0.  50  cps  also.  But  at  the  two  highest 
frequencies,  the  situation  changes  slightly.  MCF-3  has  generally  sharper 
falloff  from  high  velocity.  The  total  rejection  at  velocities  below  2  km/sec 
is  hard  to  determine,  but  MCF-6  appears  to  do  better  than  MCF-3.  The 
sharp  falloff  of  MCF-3  at  high  frequencies  is  offset  by  its  lack  of  isotropy. 

The  isotropic  response  at  high  velocities  should  compensate  for  the  slightly 
shallower  rejection  zone  of  MCF-6. 

Examination  of  the  noise  attenuation  plots  (Figures  IV-24, 

IV-26,  IV-27,  IV-32,  IV- 34,  and  IV-35)  allows  further  comparison  among 
the  filters.  Compared  to  the  straight  sum  process  for  sample  D,  MCF-6 
does  better  than  MCF-5  at  most  frequencies  below  1.  5  cps.  The  exceptions 
are  at  0.  16  cps  where  MCF-5  is  about  1  db  better  than  MCF-6  and  at  0.  79  cps 
where  MCF-5  is  again  about  1  db  better  than  MCF-6.  Around  1.  5  cps,  MCF- 
6  has  a  much  smaller  peak,  the  difference  being  about  5  db.  Above  1.5  cps, 
the  two  filters  perform  similarly  at  frequencies  up  to  5.  0  cps.  The  compari¬ 
son  of  the  two  filters  for  sample  A  is  about  the  same  as  that  described  for 
sample  D. 
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Again  comparing  MCF/ST  SUM  ratios  for  sample  D,  MCF-3 
has  about  4.  5  db  more  attenuation  at  0.  16  cps  than  does  MCF-6  and  about 
1 . 2  db  more  at  0.32  cps.  At  0.48  cps  the  two  filters  perform  about  the 
same,  but  above  0.5  cps  MCF-6  has  more  attenuation  than  MCF-3,  except 
at  0.79  cps.  The  improvement  of  MCF-6  over  MCF-3  is  about  5  db  from 
2.  0  cps  to  5.  0  cps.  Again,  for  sample  A  similar  results  are  obtained. 

This  difference  in  noise  attenuation  at  high  frequencies  com¬ 
pared  to  a  straight  sum  process  between  ring -model  and  maximum -likelihood 
filters  is  easily  explained.  As  mentioned  above,  most  of  the  filters  have 
from  4  to  6  db  less  rejection  than  does  the  straight  sum  process  at  high  fre¬ 
quencies.  However,  those  filters  are  ring-model  filters  for  which  a  few 
traces  were  straight  -  summed  before  input.  The  resulting  straight-sum 
noise  rejection  cannot  be  directly  calculated  because  the  number  of  traces 
summed  into  a  ring  varies  from  ring  to  ring.  But  considering  an  average 
of  about  3  traces  per  ring  for  the  8  -  ring  filters,  one  would  expect  about  5  db 
more  random  noise  rejection  from  an  8-ring  filter  than  from  a  25-channel 
filter.  Comparing  the  MCF/ST  SUM  ratio  for  MCF-3  with  that  for  MCF-4, 
which  is  an  8-ring  filter  designed  from  an  infinite  velocity  signal  model, 
a  difference  of  about  5  db  is  observed  at  high  frequencies. 

Comparing  MCF/ST  SUM  ratios  of  MCF-6  for  samples  A  and 
D  shows  about  the  same  attenuation  for  both  samples  at  0.  16  cps  (0.  2  db 
more  attenuation  for  sample  A),  with  slightly  less  attenuation  for  sample 
A  at  other  frequencies.  The  differences  between  the  two  samples  are  about 
the  same  as  those  for  MCF-3,  except  that  at  0.  16  cps  MCF-3  has  0.  5  db 
less  attenuation  for  sample  A  than  for  sample  D. 

A  comparison  of  the  MCF/CC  ratios  for  MCF-5  and  MCF-6 
shows  the  same  results  as  the  MCF/ST  SUM  ratios.  Likewise,  comparison 
of  the  MCF/CC  ratios  for  MCF-3  and  MCF-6  shows  results  similar  to  the 
MCF/ST  SUM  ratios. 

Thus,  in  a  very  narrow,  very  low  frequency  band  (frequencies 
below  0.4  cps)  the  MCF-3  has  from  1  to  4.  5  db  greater  attenuation  than 
MCF-6.  At  higher  frequencies  MCF-6  has  greater  attenuation,  the  average 
difference  being  from  3  to  5  db  between  2.  0  and  5.  0  cps. 

Other  techniques  can  greatly  improve  the  attenuation  of  the 
5-ring  MCF  at  frequencies  below  0.4  cps,  as  will  be  shown  later. 

In  addition  to  the  response  measurements,  an  estimation  of 
signal  distortion  was  obtained  for  MCF-3  and  MCF-6  by  applying  each  filter 
to  teleseismic  events.  Figure  IV-5  shows  the  center  seismometer,  the 
straight  sum  output  and  the  MCF-3,  -6,  -11,  and  -12  output  traces  for  each 
event  used.  (Some  filters  were  not  applied  to  every  event.  )  The  distortion 
of  the  signal  seems  to  be  negligible  for  both  filters.  Another  noteworthy 
feature  is  the  superiority  of  MCF-6  over  MCF-3  with  regard  to  noise  rejection. 


IV  -  1 8 


science  services  division 


This  comparison  indicates  that  the  5 -ring  MCF  designed 
from  a  disk  signal  model  does  better  than  the  8-ring  disk  signal  model 
MCF  and  performs  comparably  to  a  maximum-likelihood  filter.  Further¬ 
more,  the  ring-model  filter  has  the  advantage  of  an  isotropic  passband  in 
the  area  of  k-space  corresponding  to  velocities  above  12  km/sec,  a  feature 
which  the  maximum-likelihood  filter  does  not  have.  These  results  led  to  the 
choice  of  the  5-ring,  disk  signal  model  MCF  over  the  maximum-likelihood 
design. 

2.  Comparison  of  Infinite  Velocity  (MCF-4)  and  Disk  Signal  Models  (MCF-5) 

MCF-4,  an  8 -ring  filter  designed  from  measured  noise  and  an 
infinite  velocity  signal  model,  was  compared  with  MCF-5,  an  8-ring  filter 
designed  from  the  same  measured  noise  as  for  MCF-4  and  from  a  disk  signal 
model  with  the  same  power  spectrum  as  the  average  of  the  5-channel  noise 
autopower  spectra. 

The  impulse  and  frequency  responses  of  MCF-4  are  shown 
in  Figure  IV-7.  The  peak  at  zero  time  has  an  amplitude  of  0.  96  and  the 
side  lobes  are  less  than  0.  03  in  magnitude.  The  mean-square  -  er  ror 
reduction  is  shown  in  Figure  IV- 15. 

The  f-k  response  of  MCF-4  is  shown  in  Figure  IV- 17.  At 
0.25  cps  the  filter  is  unable  to  cut  very  sharply,  but  does  have  a  rejection 
ring  at  a  low  wavenumber  value  which  corresponds  to  the  high-velocity 
noise  source.  At  0.50  cps  the  rejection  zone  appears  at  higher  velocity, 
probably  due  to  the  increasing  resolution.  Except  for  a  few  deep  rejection 
regions,  the^  response  is  quite  high  at  velocities  below  2.  5  km/ sec.  However, 
the  noise  f-k  spectrum  (Figure  11-18,  Section  II)  shows  very  little  noise  at 
these  velocities.  The  response  at  1.00  cps  is  similar  to  that  at  lower 
frequencies.  At  1.50  cps  the  filter  response  is  down  even  at  infinite 
velocity.  At  this  frequency  the  center  channel  noise  autopower  density 
spectrum  from  which  the  signal  model  was  formed  has  dropped  1  to  2  db 
below  the  average  noise  power  density  spectrum,  causing  a  decrease  in 
the  S/N  ratio.  Furthermore,  the  high-velocity  noise  source  is  concentrated 
in  a  very  small  area  in  k-space  and  the  resolution  is  insufficient  to  enable 
the  filter  to  reject  the  noise  without  some  rejection  at  infinite  velocity. 

MCF-5  is  identical  to  MCF-4  except  that  MCF-5  was  designed 
from  a  disk  signal  model  rather  than  from  an  infinite  velocity  signal  model. 

The  impulse  and  frequency  responses  are  shown  in  Figure  IV-8  and  were 
explained  in  the  previous  paragraphs. 
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The  first  comparison  is  that  of  the  f-k  responses  of  the  two 
filters  (Figures  IV-17  and  IV-18).  At  0.25  cps  the  MCF-5  response  falls 
more  sharply  and  goes  to  a  lower  value  in  some  directions  than  does  that 
of  MCF-4.  Furthermore,  at  velocities  below  about  1.  5  km/sec  the  MCF-4 
response  comes  up  quite  high,  with  peaks  of  412  db.  The  MCF-5  response 
has  a  few  small  -f-6-db  peaks  at  low  velocity;  however,  over  all  the  area  of 
k-space  shown,  corresponding  to  velocities  below  2  km/sec,  the  MCF-5 
response  is  considerably  lower  than  that  of  MCF -4.  At  0.  50  cps  MCF -5 
still  has  a  sharper  falloff  and  deeper  rejection  ring  between  3  and  4  km/sec. 
MCF-5  stays  lower  (between  1 . 7  and  3 .  0  km/ sec)  than  MCF -4,  but  MCF -4 
is  still  slightly  more  isotropic  than  MCF-5.  At  1.  00  cps  MCF-4  has  a 
deeper,  more  uniform  rejection  zone  of  6  to  12  db  around  2  km/sec  and 
extending  to  about  1.6  km/sec  than  MCF-5.  The  rejection  zone  of  MCF-5  is 
from  about  1.7  to  1.4  km/sec.  At  1.50  cps  a  major  problem  arises  in 
MCF-4:  the  response  is  down  more  than  6  db  over  the  entire  signal  area 
in  k-space.  This  signal  rejection  was  explained  above. 

Both  filters  have  little  rejection  at  velocities  below  8  km/sec, 
but  the  response  of  MCF-5  seems  to  be  slightly  lower  at  these  velocities. 

MCF-4  behaves  well  at  all  frequencies  but  the  highest,  al¬ 
though  it  does  not  have  as  sharp  a  response  at  high  velocities  or  as  much 
rejection  at  low  velocities  as  does  MCF-5  at  the  two  low  frequencies.  At 
1.  00  cps  MCF-4  is  about  the  same  as  MCF-5,  but  at  1.  50  cps  MCF-4 
attenuates  the  signal.  This  problem  might  be  reduced  by  the  use  of  a  signal 
model  identical  in  form  to  the  average  autopower  spectrum  rather  than  to 
an  individual  channel  so  that  there  would  be  less  variation  in  S/N  ratio.  The 
chief  advantage  of  MCF-4  over  MCF-5  is  its  greater  isotropy;  the  chief 
disadvantage,  besides  the  signal  attenuation  at  1.50  cps,  is  the  more  gradual 
falloff  of  its  response  with  increasing  velocity. 

The  noise  attenuation  plots  show  about  the  same  results  in 
comparing  the  two  filters  with  each  other,  whether  by  MCF /ST  SUM  ratios 
(Figures  IV-25  and  IV-26)  or  by  MCF/CC  ratios  (Figures  IV-33  and  IV-34) 
and  whether  using  sample  A  or  sample  D. 

At  the  two  lowest  frequencies  (0.  16  and  0.  32  cps)  the  ratios 
are  the  same  for  both  filters.  As  frequency  increases,  the  attenuation  of 
MCF-5  becomes  less  than  that  of  MCF-4  until  at  about  1.  5  cps  where  it  is 
about  1.5  to  2.  0  db  less  than  that  of  MCF-4.  From  1.5  to  5.  0  cps  the 
ratios  are  similar,  with  each  filter  alternately  having  the  greater  attenuation. 
Thus,  on  the  basis  of  noise  attenuation,  MCF-4  seems  to  do  better  than 
MCF-5. 
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The  comparison  of  these  two  filters  indicates  that  MCF-5 
maintains  a  much  better  signal  passband  at  all  frequencies  than  does  MCF-4. 
The  slight  improvement  in  noise  attenuation  of  MCF-4  over  MCF-5  in  a  small 
frequency  band  (1.5  to  2.0  cps)  is  not  sufficient  to  offset  MCF-4's  disadvan¬ 
tage  of  6  db  or  more  attenuation  in  the  signal  area  of  k-space.  Thus,  the 
disk  signal  model  results  in  a  better  filter  than  an  infinite  velocity  signal 
model. 

3.  Comparison  of  a  Theoretical  Noise  Model  (MCF-7)  and  Measured- 
Noise  Statistics  (MCF-6) 

A  5-ring  filter  (MCF-7)  designed  from  a  theoretical  annular 
noise  model  in  an  area  of  k-space  corresponding  to  an  annulus  between  the 
velocities  of  1  and  8  km/sec  was  compared  with  a  5-ring  filter  (MCF-6) 
designed  from  measured  noise.  The  same  disk  signal  model  was  used  for 
both  filte  r s . 


The  impulse  and  frequency  responses  of  MCF-7  are  shown 
in  Figure  IV-10  and  the  mean-square-error  reduction  is  shown  in  Figure  IV- 14. 
Since  both  the  signal  and  noise  are  perfectly  isotropic,  the  impulse  response 
is  perfectly  symmetrical  and  the  phase  angle  of  the  frequency  response  is 
0.  0  ±  0.  05°  everywhere.  The  location  of  the  two  doublets  on  the  impulse 

response  corresponds  to  the  half-cycle -per-unit-frequency  periodicity  in 
the  frequency  response;  their  character  as  doublets  results  from  the  in¬ 
creasing  amplitude  of  the  envelope  of  the  frequency  response. 

The  f-k  response  shown  in  Figure  IV-20  deserves  comment. 

The  two  major  differences  between  this  response  and  the  others  previously 
discussed  are  that  the  response  is  not  as  sharply  peaked  at  high  velocities 
nor  does  it  have  any  high  spots  at  low  velocities.  Both  characteristics  are 
due  to  the  difference  in  noise  model.  In  the  earlier  cases,  the  actual  mea¬ 
sured  noise  was  similar  in  shape  and  location  in  k-space  to  the  signal  model. 

In  order  to  reject  noise,  the  filter  had  to  have  a  response  which  would  fall  off 
very  rapidly  at  velocities  just  below  12  km/sec.  Furthermore,  since  there 
was  very  little  noise  at  low  velocities,  the  response  was  allowed  to  rise  at 
low  velocities,  if  necessary,  to  improve  the  rapid  falloff  at  the  higher  veloc¬ 
ities.  In  MCF-7,  there  is  complete  separation  between  signal  and  noise  and 
much  less  need  for  rapid  falloff  from  the  high-velocity  response.  Since  the 
noise  covers  the  low-velocity  region,  the  response  must  stay  down  at  these 
velocities . 
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The  noise  attenuation  plots  for  MCF-7  are  shown  in  Figures 
IV-28,  IV  -  36  and  IV-43.  The  shape  of  these  plots  is  different  from  the  ones 
previously  examined.  The  MCF/ST  SUM  ratios  show  that  MCF-7  is,  at 
most,  only  about  2  or  3  db  better  than  a  straight  sum  process.  Since  the 
MCF  was  designed  from  models  which  have  all  signal  and  no  noise  at  high 
velocity,  it  will  make  no  effort  to  attenuate  actual  noise  at  high  velocity. 

But  the  two  noise  samples  to  which  the  filter  was  applied  contain  mostly 
high-velocity  noise  at  low  frequency.  Thus,  the  filter  will  not  attenuate 
much  of  the  noise  in  these  samples. 

The  responses  of  MCF-6  were  explained  in  a  preceding 
paragraph.  Comparing  the  f-k  responses  of  the  two  filters  shows  that  the 
MCF-7  response  falls  off  with  decreasing  velocity  much  more  slowly  than 
that  of  MCF-6.  The  broad  response  lobe  at  high  velocity  would  pass  a  large 
amount  of  noise.  The  performance  of  MCF-7  at  low  velocities  is  markedly 
superior  to  that  of  MCF-6,  but  previous  analysis  has  pointed  out  the  lack  of 
major  low-velocity  noise  sources. 

The  noise  attenuation  plots  (Figures  IV-27 ,  IV-28,  IV-35,  IV-36, 
IV-42,  and  IV-43)  show  that  at  the  very  low  frequencies  MCF-6  has  about 
4  to  5  db  more  attenuation  than  does  MCF-7. 

— 

In  summary,  MCF-7  lacks  the  sharp  falloff  in  f-k  response 
from  high  velocities  needed  for  the  actual  noise  field,  and  has  less  noise 
attenuation  than  does  MCF-6  at  low  frequencies  where  there  is  the  bulk  of 
the  noise.  These  factors  indicate  that  measured -nois e  statistics  should  be 
used  in  the  design  of  the  filters. 

4.  Comparison  Between  the  Use  of  Nonwhitened  Noise  (MCF-4)  and  the 
Use  of  Whitened  Noise  (MCF-8) 

An  8-ring  filter  (MCF-4)  designed  from  an  infinite  velocity 
signal  model  and  measured  noise  was  compared  with  an  8-ring  filter  (MCF-8) 
designed  from  an  infinite  velocity  signal  model  and  whitened  noise  obtained 
by  whitening  the  noise  correlations  used  in  designing  MCF-4.  The  infinite 
velocity  signal  model  used  for  MCF-8  was  obtained  from  the  whitened  noise 
autocorrelation  for  channel  1. 

The  impulse  and  frequency  responses  of  MCF-8  are  shown 
in  Figure  IV  -  1  1  and  the  mean  -  square  -  e  r  ro  r  reduction  is  shown  in  Figure 
IV- 15.  The  amplitude  response  has  a  variation  of  less  than  approximately 
0.  25  db  between  0.  0  and  5.  0  cps.  This  small  variation  indicates  that  the 
use  of  whitened  noise  reduces  the  S/N  ratio  variations  which  allow  frequency¬ 
filtering  . 
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The  f-k  response  in  Figure  IV-21  shows  a  peak  at  infinite 
velocity  surrounded  by  a  rejection  ring  of  12  to  24  db  between  2  and  3  km/sec 
at  0.  25  cps.  The  peak  narrows  and  the  ring  moves  in  at  0.  50  cps  and  again 
at  1.  00  cps,  though  the  ring  is  breaking  up  and  goes  down  only  to  -18  db. 
Finally  at  1.  50  cps  the  first  ring  is  only  down  12  db,  but  the  falloff  from  high 
velocity  is  very  sharp. 

At  low  velocities  below  4  km/sec  the  response  is,  in  general, 
lower  than  for  the  nonwhitened  noise  filters.  At  0.  50  cps  this  response  is  up 
around  -3  to  +3  db  over  much  of  the  area  between  1  and  2  km/sec  where  there 
is  little  noise.  However,  at  other  frequencies  the  response  is  as  high  as  0  db 
only  twice;  both  of  the  peaks  are  very  narrow  with  rapid  falloff  and  occur  at 
1.  00  cps  at  a  velocity  less  than  1  km/sec. 

The  noise  attenuation  plots  are  shown  in  Figures  IV-29  and 
IV -37.  The  chief  difference  between  these  ratios  and  those  of  the  nonwhitened 
noise  filters  is  that  these  ratios  are  relatively  flat  over  the  entire  frequency 
range.  The  MCF/ST  SUM  ratio  for  sample  D  shows  about  4.5-db  improvement 
for  the  MCF  at  the  two  lowest  frequencies  measured  (0.  16  and  0.32  cps)  and 
shows  4.  5  and  3. 7  db  for  the  two  lowest  frequencies  on  sample  A.  The 
MCF/CC  ratios  show  that  the  MCF  deli  vers  6.  2-  and  7.  0-db  improvement 
over  the  center  channel  on  sample  D  at  the  two  lowest  frequencies  while  it 
gives  6.5  db  at  each  frequency  on  sample  A. 

The  impulse  and  frequency  responses  (Figure  IV -7)  and  the 
f-k  response  (Figure  IV-  17)  of  MCF-4  were  described  in  paragraph  C.  2. 

An  examination  of  the  f-k  responses  (Figures  IV -17  and  IV-21) 
of  the  two  filters  shows  that  the  response  of  MCF-8  falls  off  from  high  veloc¬ 
ities  more  rapidly  than  that  of  MCF-4.  Also,  MCF-8  has  a  noticeably  deeper 
rejection  zone  around  the  high-velocity  peak.  This  rejection  is  good  even  at 
1.50  cps  where  MCF-4  has  no  rejection  between  2  and  8  km  /  sec.  MCF-8 
has  a  response  peak  at  infinite  velocity  rather  than  a  rejection  area  as  has 
MCF-4. 

The  f-k  response  of  MCF-8  seems  to  be  much  better  than 
that  of  MCF-4,  maintaining  a  sharp  response  peak  at  high  velocities  and  at 
all  frequencies  considered.  It  should  be  noted  that  with  this  infinite  velocity 
signal  model  obtained  from  a  correlation  set  identical  in  form  to  that  of  the 
center  channel  noise  autocorrelation,  the  filters  were  unable  to  maintain  a 
high  response  over  the  entire  k-space  area  within  the  12-km/sec  circle 
at  1.50  cps . 
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The  noise  attenuation  plots  are  shown  in  Figures  IV-25,  IV-29, 
IV-33,  and  IV-37.  A  comparison  of  the  MCF/ST  SUM  ratios  for  the  two  filters 
shows  that  MCF-8  has  more  attenuation  than  MCF-4  at  all  frequencies  except 
the  first  (0.  16  cps).  At  0.  16  cps  the  attenuation  of  MCF-4  is  3  db  better  for 
sample  A  and  2 .  5  db  better  for  sample  D.  F rom  1.0  to  5.0  cps  MCF-8  has 
an  average  of  about  7-db  improvement  over  MCF-4  for  sample  D  and  about 
5  db  for  sample  A.  Similar  results  are  obtained  from  the  MCF/CC  ratio  com¬ 
parison  except  that  at  0.32  cps  for  sample  D,  MCF-4  is  still  better  than  MCF-8, 
but  only  by  about  0.25  db.  Also,  at  5.0  cps  MCF-4  has  a  sharp  dip  in  the  ratio 
which  goes  about  0.  75  db  below  that  of  MCF-8  for  sample  D  and  about  2.  25  db 
below  for  sample  A.  This  dip  was  explained  in  the  earlier  paragraph  on  noise 
attenuation  (paragraph  B.4). 

These  comparisons  indicate  that  substantial  (3  to  10  db)  improve¬ 
ment  in  noise  attenuation  was  obtained  over  almost  the  entire  frequency  range 
of  interest  by  whitening  the  noise  correlations  before  designing  the  filter. 
Furthermore,  significant  improvement  in  the  f-k  response  resulted  from  the 
whitening  process.  These  results  indicate  that  the  noise  used  in  the  filter 
design  should  be  whitened. 

5.  Comparison  of  Nighttime  (MCF-9)  and  Daytime  (MCF-10)  Noise 
Statistics 

This  study  was  made  to  determine  if  noise  reduction  could  be 
significantly  improved  by  applying  a  filter  designed  from  nighttime  noise 
samples  to  nighttime  noise  and  applying  a  filter  designed  from  daytime  noise 
samples  to  daytime  noise.  Two  5-ring  filters  were  designed  from  infinite 
velocity  signal  models,  but  from  different  noise  samples.  The  nighttime 
filter  (MCF-9)  was  designed  from  an  average  of  seven  noise  samples  re¬ 
corded  between  01:00  and  06:00  MST  while  the  daytime  filter  (MCF-10)  was 
designed  from  an  average  of  four  noise  samples  recorded  between  11:00 
and  17:00  MST. 


The  impulse  and  frequency  responses  of  MCF-9  shown  in 
Figure  IV- 12  indicate  that  this  technique  of  averaging  correlations  over 
several  noise  samples  produces  a  filter  less  prone  to  fr equency-filte ring 
than  does  the  single  sample  design.  MCF-9  has  an  amplitude  response 
flat  within  40  and  about  -0.5  db.  The  mean-square -error  reduction  is 
shown  in  Figure  IV-15. 

The  f-k  response  of  MCF-9  shown  in  Figure  IV-22  follows 
the  pattern  already  described  for  filters  designed  from  a  single  nonwhitened 
noise  sample  at  0.  25,  0.  50  and  1. 00  cps.  A  slight  difference  is  noted  at 
1.50  cps.  The  previously  mentioned  filters  had  trouble  separating  signal 
and  noise  at  high  velocities  because  the  high-velocity  noise  peak  was  so 
narrow  that  it  almost  matched  the  signal  model.  In  order  to  reject  the  noise, 
the  filter  either  had  to  reject  extremely  high  velocities,  resulting  in  a  re¬ 
sponse  which  was  down  3  to  6  db  at  infinite  velocity  as  in  the  case  of  the 
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infinite  velocity  filter  (MCF-4),  or  had  to  compromise  and  pass  signal, 
resulting  in  less  than  3-db  rejection  at  lower  velocities  (8  to  1Z  km/sec). 
The  process  of  averaging  the  correlations  produced  a  model  with  a  more 
uniform  spectrum  over  the  frequency  range  of  interest.  Thus,  the  high- 
velocity  noise  peak  in  k-space  does  not  seem  to  narrow  down  drastically 
at  1.50  cps  as  it  did  with  the  single  sample.  The  resulting  filter  response 
remains  above  -3  db  over  most  of  the  area  within  the  lZ-km/sec  velocity 
ring,  but  drops  as  low  as  -18  db  at  8  km/sec. 

The  noise  attenuation  plots  are  shown  in  Figures  IV-30 
and  IV-38.  MCF-9  and  MCF-10  (applied  to  a  different  pair  of  noise 
samples  than  were  the  other  filters)  were  designed  from  an  average  of 
noise  samples  and  were  applied  to  one  sample  of  each  type.  Sample  A 
was  the  sample  used  for  the  other  filters  and,  being  recorded  between 
03:18  and  03:21  MST,  was  chosen  as  the  nighttime  noise  sample.  Though 
not  used  in  the  design  of  MCF-9,  sample  A  is  immediately  adjacent  to 
one  of  the  seven  noise  samples  used.  In  the  design  of  MCF-10,  four  day¬ 
time  noise  samples  of  150-sec  duration  were  used.  One  of  these  samples 
was  chosen  to  represent  the  daytime  noise  in  the  measurement  of  noise 
attenuation.  This  daytime  sample  (sample  E)  is  a  150-sec  section  of  a 
noise  sample  recorded  between  13:58  and  14:01  MST. 

In  summary,  the  noise  attenuation  plots  show  the  results  of 
applying  a  filter  (MCF-9)  designed  from  nighttime  noise  to  a  nighttime 
noise  sample  (sample  A)  and  to  a  daytime  noise  sample  (sample  E). 

A  comparison  of  the  MCF/ST  SUM  ratios  for  the  two  samples 
shows  that  at  0.  16  cps  the  nighttime  filter  improvement  in  noise  rejection 
over  a  straight  sum  process  is  about  0.  5  db  better  for  the  nighttime  sample 
than  for  the  daytime  sample.  At  0.  32  and  0.  48  cps  the  behavior  for  both 
samples  is  the  same.  Above  0.48  cps  the  MCF/ST  SUM  ratio  for  the  night¬ 
time  sample  becomes  increasingly  better  over  a  small  band  of  frequencies 
than  does  that  for  the  daytime  sample  with  a  maximum  improvement  of 
about  5  db  from  1.0  to  1.5  cps.  Above  1.5  cps  the  two  ratios  approach 
the  same  value  at  about  1. 9  cps  and  closely  follow  each  other  up  to  5.  0  cps. 

Figure  IV- 13  shows  the  impulse  and  frequency  responses 
of  MCF-10.  This  filter  has  an  amplitude  response  flat  within  40  and  about 
-0.  5  db  from  0.  0  to  5.  0  cps.  The  mean-square-error  reduction  is  shown 
in  F igur  e  IV  -  1 5  . 

The  f-k  response  shown  in  Figure  IV -23  is  very  similar  to 
that  of  MCF-9  except  that  at  1.50  cps  MCF-10  has  a  rejection  ring  around 
7  km/sec  which  is  12  to  18  db  below  that  for  MCF-9.  Other  features  of 
the  f-k  response  follow  the  patterns  already  discussed  for  the  other  filters. 
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Noise  attenuation  plots  are  shown  in  Figures  IV -31  and  IV-39. 
At  0,  16  and  0.  32  cps  the  daytime  filter  MCF/ST  SUM  ratio  is  about  1  db 
better  for  the  daytime  sample  than  that  for  the  nighttime  sample.  From 
about  0.  5  to  1.  0  cps  the  daytime  filter  works  better  on  the  nighttime  sample, 
the  ratio  showing  about  3-db  improvement  over  the  daytime  sample  at  about 
0.  75  cps.  From  1. 2  to  5.  0  cps  the  two  ratios  follow  each  other  in  general. 

Comparing  the  f-k  responses  of  the  two  filters  (Figures 
IV-22  and  IV-23),  they  are  similar  at  0.25  cps,  except  that  the  response 
of  MCF-9  remains  high  over  a  wider  area  of  k-space  at  high  velocity,  but 
drops  more  sharply  to  a  lower  value  than  does  that  of  MCF-10.  The  response 
of  MCF-9  is  generally  lower  at  the  low  velocities  than  is  the  response  of 
MCF-10,  At  0.  50  cps  MCF-9  keeps  its  peak  at  high  velocity  sharper  than 
does  MCF-10,  but  the  MCF-9  response  does  not  fall  as  sharply  or  fall  to 
such  low  values  in  its  first  rejection  ring  as  does  that  of  MCF-10.  However, 
at  velocities  below  3  km/sec  MCF-9  has  a  generally  lower  response  than 
does  MCF-10.  At  1.  00  cps  the  response  of  MCF-9  is  higher  over  most  of 
k-space  than  is  the  response  of  MCF-10.  Between  2  and  3  km/sec  MCF-10 
seems  to  be  trying  to  reject  some  noise  whereas  MCF-9  does  nothing  notice¬ 
able  in  this  frequency.  At  1.  50  cps  MCF-9  exhibits  a  sharp  falloff  below 
12  km/sec  and  a  deep  rejection  ring  between  6  and  8  km/sec.  In  achieving 
this  deep  rejection,  the  response  was  allowed  to  come  up  quite  high  around 
3  km/sec.  MCF-10,  not  having  the  deep  rejection  ring,  does  not  rise  quite 
as  high  as  3  km/sec. 

These  responses  indicate  that  the  noise  does  vary  from  day 
to  night  or  from  month  to  month:  the  nighttime  noise  samples  were  recorded 
in  November  1965  while  the  daytime  samples  were  recorded  four  months 
later  in  late  March  and  early  April  1966. 

The  noise  attenuation  plots  are  shown  in  Figures  IV-30, 

IV  -  3  1 ,  IV  -38  ,  and  IV-39.  A  comparison  of  the  MCF/ST  SUM  ratios  shows 
that,  in  general,  the  daytime  filter  (MCF-10)  performs  better  on  the  daytime 
noise  (sample  E)  than  does  the  nighttime  filter;  however,  the  difference  is 
only  1  db  at  0.  16  cps.  At  0.  32  cps  both  filters  are  the  same,  but  the  night¬ 
time  filter  is  about  1  db  better  at  0.48  cps.  At  higher  frequencies,  the 
daytime  filter  is  slightly  better  than  the  nighttime  filter,  with  about  5-  or 
6-db  improvement  from  1  to  2  cps.  Considering  the  nighttime  noise 
(sample  A),  the  nighttime  filter  is  almost  1  db  better  than  the  daytime 
filter  at  0.  16  cps,  about  0.  5  db  better  at  0.  32  cps,  about  0.  5  db  better  at 
0.48  cps,  and  generally  worse  at  higher  frequencies.  The  MCF/CC  ratios 
agree  that  the  nighttime  filter  does  better  than  the  daytime  filter  in  filtering 
daytime  noise  at  a  couple  of  frequencies,  being  a  fraction  less  than  1  db 
better  at  0.  16  cps  and  0.48  cps.  At  the  other  frequencies  the  daytime 
filter  does  better. 
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At  the  low  frequencies  below  0.  5  cps,  both  filters  do  com¬ 
parable  work  on  both  samples.  At  the  frequencies  above  1.  0  cps,  the 
daytime  filter  seems  to  do  better  on  both  samples.  This  suggests  that  the 
daytime  filter,  unlike  the  nighttime  filter,  was  designed  from  noise  samples 
containing  slightly  more  power  at  the  high  frequencies. 

On  the  whole,  the  difference  in  performance  between  the 
two  filters  does  not  seem  sufficient  to  warrant  separate  filters  for  daytime 
and  nighttime  use.  In  fact,  since  the  daytime  samples  were  recorded  about 
four  months  after  the  nighttime  samples,  it  is  further  concluded  that  a  filter 
designed  from  an  average  of  several  noise  samples  may  be  expected  to 
perform  consistently  well  over  a  period  of  several  months.  These  results 
indicate  that  the  filters  should  be  designed  from  an  average  of  several 
noise  samples. 

6.  Comparison  of  Constant  (MCF -  6)  and  Variable  (MCF-11) 
Signal-To-Noise  Ratios 

Multichannel  filters  described  in  previous  paragraphs  have 
utilized  signal  models  whose  power  spectra  have  been  shaped  to  approxi¬ 
mate  the  spectrum  of  the  noise  model.  This  shaping  procedure  is  normally 
followed  to  insure  that  the  filter  has  a  good  transient  response.  It  also 
insures  that  MCF  rejection  of  a  particular  frequency  is  on  a  velocity  basis 
rather  than  due  to  spectral  differences  between  signal  and  noise.  A  constant 
multiplier  is  usually  applied  to  the  signal  model  to  achieve  a  desired  S/N 
ratio.  The  following  study  was  initiated  to  determine  the  extent  to  which 
detection  capabilities  could  be  enhanced  by  allowing  the  MCF  to  reject  on 
both  a  frequency  and  velocity  basis. 

MCF  -  6  and  MCF-11  were  used  to  perform  this  evaluation. 
They  were  identical  in  design,  except  for  the  shape  of  the  signal  model 
power  spectrum.  Both  are  1.  0  sec  long  and  are  designed  for  the  5-channel 
ring  model  case.  The  noise  model  in  each  case  was  from  a  150-sec  section 
of  data  recorded  on  24  December  1965  at  approximately  03:00  MST  (noise 
sample  D) .  No  equalization  to  incorporate  instrument  response  variations 
was  performed.  Both  filters  used  a  disk  signal  model  designed  to  pass  the 
k-space  region  from  12  km/sec  to  infinite  velocity.  This  type  of  signal 
model  introduces  less  P-wave  distortion  in  the  medium- velocity  range 
(12  to  20  km/sec)  than  does  infinite  velocity  signal  models  which  generally 
are  6  to  12  db  down  at  these  velocities.  MCF-6  was  described  previously 
in  this  section,  being  used  in  a  comparison  between  a  measured-noise  model 
and  a  theoretically  derived  noise  model.  Also,  figures  pertaining  to  MCF-6 
were  presented  and  will  be  referred  to  during  this  discussion. 
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MCF-6  used  a  constant  signal -to-organized-noise  ratio  of  4. 
The  signal  model  spectrum  was  shaped  to  resemble  the  noise  model  spectrum 
at  all  frequencies.  Therefore,  this  MCF  would  be  expected  to  achieve  noise 
rejection  primarily  on  a  velocity  basis.  This  is  verified  by  the  wavenumber 
responses  previously  pre  sented  in  Figure  IV-19.  Arrivals  with  apparent 
horizontal  velocities  above  1Z  km/sec  are  essentially  unattenuated  at  all  fre¬ 
quencies.  Noise  rejection  occurs  primarily  in  the  velocity  range  from  4  to 
8  km/ sec  for  all  frequencies.  An  additional  reject  occurs  between  1  and 
Z  km/  sec. 


MCF- 11  used  a  theoretical  shaping  function  for  the  signal 
model  passing  the  frequency  band  from  1.  0  to  Z.  Z5  cps.  The  theoretical 
input  function  is  shown  on  Figure  IV-44  overlaid  with  the  amplitude  response 
of  the  weighted  stack  of  the  filters.  An  additional  S/N  factor  of  Z  was  used 
to  insure  that  signal  power  would  exceed  noise  power  at  all  points  in  the 
k-space  passband.  The  filter  is  allowed  to  reject  noise  on  a  frequency  basis 
at  all  frequencies  outside  the  defined  passband.  Inside  this  band,  rejection 
is  a  combination  of  velocity-  and  frequency- filtering  which  varies  as  a  function 
of  the  relative  power  in  the  signal  and  noise  spectra  at  any  frequency. 

Figure  IV-44  is  a  three-part  figure  showing  the  impulse 
response  of  MCF- 11  and  the  amplitude  and  phase  responses  obtained  by 
taking  the  Fourier  transform  of  the  impulse  response.  As  the  parts  of 
this  figure  indicate,  the  MCF  response  resembles  a  bandpass  filter  in  many 
respects.  Amplitude  response  peaks  in  the  higher  frequency  reject  regions 
are  apparently  due  to  low  spectral  power  in  the  noise  spectrum  and  to 
spectral  window  effects.  The  spectral  window  is  created  when  the  inverse 
transform  of  the  frequency-domain-shaped  spectrum  is  taken.  The  exact 
inverse  transform,  in  this  case,  requires  a  ±Z-sec  time  domain  function. 
However,  this  function  was  truncated  to  a  il-sec  function  for  the  signal  model. 
Also,  the  reject  regions  of  the  signal  model  were  restrained  to  a  -30  db  level. 

Figure  IV-45  shows  the  mean-square- error  reduction 
achieved  in  the  design  of  MCF-11.  Total  noise  power  was  reduced  to 
less  than  1  percent  of  the  input  mean-square  value.  The  mean-square-error 
curve  for  MCF-6  (Figure  IV- 14)  indicated  total  noise  power  reduced  to  5.  5 
percent  of  the  original.  These  figures  assume  no  signal  distortion.  Obviously, 
total  noise  power  has  been  reduced  to  an  additional  6  db  with  MCF-11. 

Wavenumber  responses  for  MCF-11  (Figure  IV-46)  indicate 
attenuation  greater  than  1Z  db  at  all  points  in  k-space  for  frequencies  located 
outside  the  signal  model  passband.  For  example  at  0.  Z5  cps,  infinite  velocity 
energy  is  down  48  db  from  unity  and  falls  off  further  to  a  low  of  60  db  between 
8  and  1Z  km/ sec.  All  points  in  the  k-plane  are  at  least  Z4  db  down  from  unity. 
Comparatively,  MCF-6  (Figure  IV-19)  passes  energy  at  velocities  greater 
than  1Z  km/sec  with  less  than  3-db  rejection  and  is  generally  18  db  down 
between  Z  and  4  km/sec. 
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At  0.5  cps  MCF-11  is  down  at  least  9  db  over  the  entire 
k-plane,  including  the  high-velocity  region,  and  ranges  to  48  db  down 
between  1  and  2  km/sec.  MCF-6  is  less  than  3  db  down  from  infinity  to 
8  km/sec  with  rejection  regions  18  to  30  db  down  between  2  and  4  km/sec. 

In  the  passband  of  the  signal  model  at  1.00  cps,  MCF-11 
is  within  3  db  of  unity  to  8  km/sec  and  down  12  to  30  db  between  2  and  4km/ 
sec.  MCF-6  is  similar,  although  the  region  between  2  and  4  km/sec  is 
less  than  12  db  down. 

— ¥ 

At  1.5  cps  MCF-6  rejects  very  little  in  the  k-plane  at  veloc¬ 
ities  greater  than  2  km/sec.  MCF-11  still  passes  the  disk  from  12  km/sec 
to  infinity  with  less  than  3-db  attenuation,  but  generally  rejects  the  region 
between  2  and  8  km/sec  by  amounts  varying  from  12  to  24  db. 

Summarizing  the  wavenumber  response  comparisons,  it 
appears  that  allowing  MCF-11  to  filter  the  spectral  regions  of  high-noise 
energy  on  a  frequency  basis  has  enabled  that  filter  to  concentrate  its  capa¬ 
bility  to  velocity-f ilte r  on  the  low-velocity  noise  energy  contained  in  the 
passband,  thereby  achieving  improved  overall  nois e- rej ection  capability. 

The  passband  is  designed  to  encompass  the  frequency  range  expected  for 
P-phase  body-wave  energy,  thus  signal  energy  should  be  passed  without 
significant  attenuation  from  12  km/sec  to  infinity. 

Both  filters  were  applied  to  noise  samples  A  and  D. 

Figure  IV-47  illustrates  a  section  of  noise  from  each  sample,  showing  the 
center  seismometer  (ring  1),  the  straight  sum  output  and  the  filtered  outputs 
from  MCF-6,  -11  and  -12.  MCF- 12  will  be  discussed  later.  The  MCF-6 
output  shows  a  6-  to  9-db  reduction  in  the  dominant  5-sec  microseism  com¬ 
pared  with  ring  1  and  the  straight  sum.  MCF-11  has  removed  all  visible 
evidence  of  this  low-frequency  energy. 

Power  spectra  were  computed  for  ring  1,  the  straight  sum 
and  the  MCF-6  and  MCF-11  outputs  over  the  total  noise  interval  for  both 
noise  samples.  The  relative  noise  rejection  as  a  function  of  frequency  was 
computed  by  forming  the  ratio  of  power  density  at  each  frequency  between 
each  MCF  output  trace  and  the  center  seismometer  (ring  1).  Also  computed 
was  the  ratio  between  each  output  MCF  trace  and  the  straight  sum  trace. 

The  ratios  for  MCF-6  are  shown  in  Figures  IV-35  and  IV-27  and  those  for 
MCF-11  are  in  Figure  IV-48. 

The  ratios  allow  comparison  of  MCF  performance  over  noise 
not  included  in  the  MCF  design.  Superior  noise  rejection  is  evident  at  all 
frequencies  for  MCF-11  when  compared  with  MCF-6.  The  dominant  5-sec 
microseism  is  especially  attenuated,  with  MCF-11  yielding  20  db  greater 
rejection  than  MCF-6  over  this  region.  Performance  is  consistent  in  and 
out  of  the  design  interval,  indicating  a  relative  insensitivity  to  minor  noise 
field  variations  with  the  ring  models.  The  increased  noise  rejection  achieved 
by  MCF-11  (  as  described  by  these  ratios)  is  consistent  with  the  mean-square- 
error  and  frequency-wavenumber  plot  indications. 
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Signal  distortion  was  evaluated  by  apply  MCF-6  and  MCF-11 
to  3  teleseismic  events.  Figure  IV-5  shows  the  center  seismometer,  the 
straight  sum  output  and  the  MCF-3,  -6,  -11  and  -1Z  output  traces  for  each 
event.  Of  particular  interest  was  relative  noise  rejection,  signal  distortion 
(both  amplitude  and  phase)  and  possible  precursor  effects  introduced  by  the 
frequency-filtering  capabilities  of  MCF-11.  Noise  rejection  is  markedly 
superior  for  MCF-11  compared  with  MCF-6,  and  signal  amplitude  and  phase 
distortion  seem  negligible  for  both.  A  slight  precursor  is  discernible  for 
MCF-11;  however,  it  does  not  appear  to  be  of  major  concern.  For  purposes 
of  event  detection,  a  slight  precursor  is  of  relatively  small  consequence  if  a 
large  gain  in  S/N  ratio  can  be  achieved. 

In  conclusion,  it  appears  that  event  detection  capabilities  can 
be  increased  substantially  by  the  use  of  an  MCF  incorporating  a  combination 
of  velocity-  and  f r equency-filtering  characteristics  as  opposed  to  velocity¬ 
filtering  alone. 

7.  Comparison  Between  Equalized  and  Unequalized  Signal  Models 

MCF-  1  1  was  designed  using  an  idealized  signal  model,  i.  e.  , 
the  assumption  was  made  that  amplitude  and  phase  responses  were  identical 
for  each  channel.  The  noise  model,  having  been  constructed  from  measured 
noise,  contains  departures  from  this  idealized  case  due  to  existing  variations 
in  the  instrument  responses.  Depending  upon  their  extent,  such  variations 
were  utilized  by  MCF-11  to  achieve  noise  rejection.  However,  signals 
recorded  by  these  instruments  were  affected  by  the  same  variations  and 
were  also  distorted  to  some  extent  by  the  MCF. 

MCF- 12  was  designed  in  a  fashion  similar  to  MCF-11;  both 
used  band-limited  disk  signal  models  to  pass  the  k-space  area  from  12  km/sec 
to  infinity,  both  used  measured  noise  from  subarray  F-3  (noise  sample  D) 
and  both  had  S/N  ratios  of  2.  However,  the  signal  model  for  MCF- 12  was 
perturbed  to  incorporate  the  F-3  instrument  response  variations  described 
in  Section  III.  Therefore,  signals  also  containing  these  same  variations 
should  be  undistorted  by  the  MCF. 

The  design  of  MCF- 12  resulted  in  a  total  noise  power  de¬ 
crease  to  1.8  percent  of  the  input  mean-square  value  (Figure  IV-45)  as 
compared  with  1.  0  percent  obtained  with  MCF-11.  The  difference  is  a 
result  of  prohibiting  MCF- 12  from  using  the  noise  variations  due  to  instru¬ 
ment  responses  as  a  basis  for  rejection.  However,  actual  S/N  improve¬ 
ment  with  MCF  -  12  should  be  comparable  to  that  for  MCF  -  1  1  due  to  the 
expected  decrease  in  signal  distortion. 
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The  evaluation  of  the  impulse  response  for  MCF-12  required 
a  departure  from  previously  used  methods.  The  impulse  response  was  ob¬ 
tained  by  convolving  the  filter  for  each  channel  with  an  impulse  that  had  been 
perturbed  to  incorporate  the  instrument  response  variations  for  the  channel 
and  then  summing  the  convolved  outputs.  Figure  IV-49  shows  the  impulse, 
amplitude  and  phase  response  of  MCF-12.  The  signal  model  perturbations 
have  caused  considerable  character  change  from  the  MCF-11  impulse  response 
(Figure  IV-44).  The  bandpass  character  remains  evident  in  the  amplitude 
response,  although  superimposed  on  a  curve  describing  the  instrument 
response.  Some  spectral  window  effects  are  also  present  due  to  the  signal 
model  truncation.  (The  signal  model  was  shaped  in  the  frequency  domain 
and  the  exact  inverse  transform  required  a  ±2 -sec  time  domain  representation). 

The  signal  distortion  (2  to  3  db)  in  the  passband  apparently  is 
caused  by  several  factors:  the  truncation  of  the  time  domain  signal  model  to 
±1  sec,  a  relatively  low  S/N  ratio  of  2:1  and  the  rather  narrow  signal  band¬ 
width  with  high  roll-off  rates  on  both  sides.  The  low  overall  S/N  ratio 
causes  difficulty  when  high-velocity  noise  is  concentrated  in  a  small  section 
of  the  signal  passband,  creating  a  situation  where  the  S/N  ratio  is  actually 
less  than  1  and  the  MCF  is  caused  to  reject.  The  existence  of  such  noise  is 
indicated  by  the  behavior  of  the  MCF-3  f-k  response  (Figure  IV-16)  which 
suggests  a  high-velocity  noise  peak  in  the  northeast  quadrant. 

Frequency -wavenumber  plots  for  MCF-12  are  presented  in 
Figure  IV-50.  The  departure  from  ring-model  symmetry  is  a  result  of  the 
perturbation  of  the  signal  model  to  include  the  instrument  response  variations. 
Considerably  greater  noise  rejection  is  evident  with  MCF-11  (Figure  IV-46) 
at  0.  25  cps  (40  db  more  rejection  at  8  km/sec,  and  18  db  more  at  infinite 
velocity).  Both  filters  are  comparable  at  0.  5  cps  except  at  infinite  velocity 
where  MCF-12  is  6  db  down  from  MCF-11.  At  1.  0  cps  both  filters  exhibit 
some  signal  distortion  across  the  12-km/sec  disk,  with  MCF-12  down  3  db 
lower  than  MCF-11.  The  multichannel  filters1  responses  are  comparable 
through  the  remainder  of  the  k-plane.  At  1.  5  cps,  responses  for  both 
filters  are  very  similar,  indicative  of  the  relatively  less  coherent  noise 
field  at  this  frequency. 

In  summary,  the  f-k  comparisons  indicate  a  reduced  capability 
in  noise  rejection  accompanied  by  increased  signal  distortion  in  the  passband 
for  MCF-12.  The  former  effect  was  expected,  since  the  design  departure 
from  an  idealized  signal  model  eliminated  some  basis  for  noise  rejection. 
However,  the  latter  effect  is  undesirable.  Measures  taken  to  minimize 
such  distortion  in  designing  the  final  filters  include  a  broader  signal  passband 
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with  lower  roll-off  rates,  longer  (±  1.8-sec)  filters  (thus  minimizing  the 
effects  of  truncation),  and  whitening  of  the  noise  model  to  stabilize  S/N 
ratios  through  the  passband. 

MCF-1Z  was  applied  to  noise  samples  A  and  D  (Figure  IV-47). 
Its  noise  rejection  appears  to  be  roughly  comparable  with  that  for  MCF-11. 
Spectral  ratios  were  computed  for  the  MCF-1Z  output  with  both  the  center 
seismometer  (ring  1)  and  the  straight  sum  output  (average)  and  are  presented 
in  Figure  IV-51.  Again,  its  noise  rejection  is  comparable  with  that  of  MCF-11 
(Figure  IV-48). 


Figure  IV-5  illustrated  the  relative  signal  distortions  obtained 
when  MCF-11  and  -1Z  were  applied  to  three  teleseismic  events.  Rejection 
of  the  noise  preceding  each  event  is  comparable  for  the  two  filters.  However, 
the  MCF-1Z  filtered  output  indicates  signal  distortion  to  the  extent  of  3  to  6 
db  of  loss  in  amplitude  as  compared  to  the  MCF-11  output.  This  is  directly 
contrary  to  the  initial  hypothesis,  viz.  ,  that  signal  distortion  would  be  de¬ 
creased  by  incorporating  the  instrument  response  variations  in  the  signal 
model.  Part  of  the  observed  distortion  occurs  due  to  the  filter  impulse  and 
frequency  responses  just  described.  An  additional  and  perhaps  equally  im¬ 
portant  cause  is  the  evidence  presented  in  Section  III  concerning  the  observed 
changes  in  the  relative  instrument  response  variations  from  November  1965 
(when  the  1.0-cps  calibrations  were  first  taken)  to  March  1966  (when  the 
multifrequency  calibrations  were  conducted).  Especially  with  F-3,  large 
changes  occurred  in  the  responses  of  instruments  during  this  period,  making 
the  responses  measured  in  March  nonrepresentative  for  other  time  periods. 
For  this  reason,  it  was  decided  to  exclude  the  response  variations  from  the 
final  filter  design. 

8.  Effect  of  Filter  Length  on  Noise  Rejection  Properties  (MCF-13, 

-14,  -15) 

This  study  was  conducted  to  determine  the  variation  in  noise 
rejection  properties  when  the  time  span  of  the  optimum  MCF  is  increased. 
LASA  subarray  F-3,  reduced  to  a  5-channel  ring  model,  was  used.  A 
measured -nois e  model  was  constructed  from  a  330-sec  sample  of  noise 
(sample  B)  recorded  on  Z4  December  1965  at  approximately  03:00  MST. 

The  signal  model  was  a  point  in  k-space  corresponding  to  an  infinite  velocity 
signal.  A  constant  signal -to -organized -noise  ratio  of  4  was  used.  Three 
filters  were  developed:  MCF-13  (1-sec  long),  MCF-14  (2-sec  long)  and 
MCF-15  (5.  65-sec  long). 
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The  three  filters  were  compared  in  terms  of  their  responses  and 
noise  rejection  capabilities.  Figure  IV-52  illustrates  the  relative  reduction  in 
the  mean-square-error  obtained  with  the  design  of  each  filter.  Total  mean- 
square-error  reduction  is  comparable  in  all  cases,  with  MCF-13  providing 
94.  9  percent  reduction  and  MCF-15  yielding  97.  5  percent  reduction.  Since 
the  initial  mean -square -error  corresponds  to  the  total  noise  power,  assuming 
no  signal  distortion,  this  represents  a  reduction  in  total  noise  power  to  5.  1 
percent  and  2,  5  percent  of  the  initial  value,  respectively. 

The  filter  transient  responses  are  presented  in  Figure  IV-53. 
The  figure  has  three  parts:  the  impulse  responses  for  each  filter,  the 
amplitude  responses  and  the  phase  responses.  The  impulse  responses  were 
formed  from  a  weighted  sum  of  the  five  individual  filters  in  each  MCF.  The 
amplitude  and  phase  responses  correspond  to  the  Fourier  transform  of  the 
impulse  responses.  MCF-15  yields  the  best  approximation  to  an  ideal 
impulse,  with  side -lobe  peak  amplitude  lower  than  1  percent  of  the  main- 
lobe  peak  amplitude. 

The  amplitude  and  phase  responses  are  correspondingly  flat 
(amplitude  within  1  db  of  unity  from  0.  0  to  8.  0  cps,  within  2  db  from  8.  0  to 
10.0  cps  and  less  than  3°  phase  shift  at  all  frequencies).  However,  MCF-14 
and  MCF-13  also  have  impulsive  responses  with  side-lobe  peak  amplitudes 
less  than  10  percent  of  main-lobe  peak  amplitudes,  and  have  correspondingly 
flat  amplitude  and  phase  responses  between  0.0  and  5.0  cps.  The  shorter 
filters  maintain  less  control  between  5.0  and  10.  0  cps,  but  maximum  de¬ 
viation  is  still  only  4  db  in  amplitude  (MCF -14)  and  20°  in  phase  (MCF -13). 

An  interesting  fact  is  that  MCF-13  maintains  a  flatter  amplitude  response 
and  allows  the  phase  to  vary  while  MCF-14  suppresses  phase  shift  at  the 
expense  of  amplitude.  Very  little  power  was  present  in  the  design  spectrum 
at  these  higher  frequencies.  Consequently,  the  multichannel  filters  con¬ 
centrated  on  increasing  S/N  ratio  at  the  lower  frequencies,  at  the  expense 
of  controlling  the  higher  frequencies. 

Figure  IV-54  presents  the  relative  power  spectrum  of  the 
center  seismometer  in  the  ring  model  for  the  design -inter val  noise  sample. 
This  figure  also  indicates  that  spectral  power  in  the  region  above  5.  0  cps 
ranges  from  30  db  to  more  than  50  db  below  peak  power.  Instrument 
response  characteristics  have  not  been  removed  from  these  curves. 

The  random-noise  response  is  presented  as  a  frequency  function 
for  each  MCF  in  Figure  IV-55.  Also  indicated  is  the  reduction  obtainable  for 
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random  noise  by  the  straight  sum  process.  Each  MCF  shows  a  rapidly 
increasing  response  with  decreasing  frequency  below  1.0  cps  where,  in 
previous  reports  ‘  ,  the  ratio  of  random  noise  to  coherent  noise  has  been 
shown  to  be  small. 

Figures  IV-56,  IV-57  and  IV-58  present  the  wavenumber 
response  evaluation  for  MCF-13,  -14  and  -15,  respectively,  at  frequencies 
0.Z5,  0.50,  1.0,  and  1.5  cps .  Regions  of  more  than  1 8  -db  r  ejection  are 
shaded.  In  general,  contour  levels  shown  are  in  6 -db  increments,  although 
the  -6  db  contour  generally  was  omitted  from  the  more  complex  plots  for 
clarity.  The  ±3-db  contour  levels  are  included  only  at  the  centers  of  each 
plot. 


In  general,  these  response  curves  show  fairly  sharp  peaks 
around  the  origin  in  k-space  which  corresponds  to  the  infinite  velocity  signal 
model.  Since  P-phase  arrivals  with  epicentral  angles  between  80°  and  100° 
will  exhibit  apparent  horizontal  velocities  greater  than  20  km/sec,  P-phase 
signal  attenuation  is  limited  in  all  cases  to  less  than  6  db  and  usually  is  less 
than  3  db. 


At  0.  25  and  0.  50  cps,  noise  rejection  is  primarily  in  the 
region  from  1  to  4  km/sec  for  all  three  filters,  although  additional  re¬ 
jection  at  4  km/sec  is  obtained  with  MCF-15  at  0.  5  cps.  The  region  from 
2  to  4  km  / sec  generally  is  r ejected  by  MCF-15  and  MCF-14  at  1.0  cps,  but 
the  shorter  MCF-13  has  a  rejection  ring  in  the  1.  5  km/sec  region.  This 
type  of  tradeoff,  or  compromise,  will  occur  when  different  filter  lengths 
are  compared,  since  each  filter  is  constrained  by  design  to  be  the  optimum 
filter  for  a  given  length.  In  this  case,  optimum  is  defined  to  mean  the  set 
of  filter  points  that  achieve  the  maximum  reduction  in  mean -square -e rror. 
Therefore,  filters  of  different  lengths  may  satisfy  this  restriction  by  re¬ 
jecting  different  regions  in  k-space,  according  to  their  capabilities.  Varia¬ 
tions  in  the  reject  regions,  however,  generally  will  occur  only  at  frequencies 
where  noise  power  is  relatively  small  over  the  whole  k-plane  and  where  such 
tradeoffs  are  reasonably  small -order.  Figure  IV -52  indicated  that  all  three 
filters  achieved  approximately  equal  reductions  in  total  noise  power.  The 
k-space  region  between  4  and  7  km/sec  is  attenuated  by  all  three  filters  at 
1 .  5  cps . 


Comparison  of  the  three  filters  in  terms  of  wavenumber 
response  leads  to  the  conclusion  that  MCF-15  achieves  the  best  rejection 
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at  both  low  and  high  frequencies.  However,  both  MCF-14  and  MCF-13 
have  comparable  responses,  and  it  is  probable  that  the  increase  in  total 
noise  rejection  is  not  sufficient  to  warrant  the  increased  cost  of  design 
and  application.  An  overall  improvement  in  noise  rejection  is  apparently 
achieved  with  increased  filter  length,  although  this  may  not  necessarily  be 
true  for  specific  frequencies  or  areas  in  k-space. 

The  three  filters  were  applied  to  noise  samples  B  (part  of 
the  design  interval)  and  C  (recorded  25  hours  earlier).  Figure  IV-59  shows 
the  center  seismometer  (ring  1)  of  the  subarray,  the  straight  sum  trace 
and  each  of  the  three  MCF  output  traces  for  the  two  noise  samples.  Little 
difference  is  apparent  in  MCF  performance  for  the  two  noise  samples, 
indicating  that  the  filters  are  not  highly  tuned  to  the  noise  contained  in  the 
design  interval.  This  relative  insensitivity  to  minor  variations  in  the  noise 
field  is  a  particularly  advantageous  feature  of  a  ring  model  which  makes  no 
assumptions  as  to  noise  directionality.  In  both  cases  MCF-15  generally 
has  rejected  the  low-frequency  noise  better  than  MCF-14  and  MCF-13,  but 
all  three  filters  have  achieved  from  6  to  12  db  greater  rejection  of  the  pre¬ 
dominate  low-frequency  noise  energy  than  has  the  straight  sum.  Cancellation 
of  the  less  -coherent  high-frequency  energy  is  roughly  equivalent  for  MCF-15 
and  the  straight  sum,  with  somewhat  less  rejection  achieved  by  MCF-14  and 
MCF-13  as  might  be  expected  after  examination  of  the  random  noise  responses 
presented  earlier.  It  should  be  noted  that  the  center  seismometer  is  buried 
about  500  ft  deep  and  the  other  seismometers  are  buried  200  ft  deep,  and  thus 
the  center  seismometer  does  not  necessarily  give  a  representative  indication 
of  the  high-frequency  energy  content  across  the  subarray  (Figure  IV-3). 

Power  density  spectra  using  a  150-sec  data  gate  were  computed 
for  each  of  the  traces  illustrated  in  Figure  IV-59.  The  relative  power  density 
spectra  for  the  input  noise  sample  (ring  1),  the  straight  sum  output  and  the 
three  MCF  outputs  are  shown  in  Figure  IV-54.  The  relative  noise  rejection 
as  a  frequency  function  was  computed  by  forming  the  ratio  of  power  density 
at  each  frequency  between  each  output  MCF  trace  and  the  center  seismometer. 
These  curves  are  shown  in  Figure  IV-60  for  noise  samples  B  and  C.  Also 
computed  was  the  ratio  between  each  output  MCF  trace  and  the  straight  sum 
output.  These  ratios  are  presented  in  Figure  IV-61  for  noise  samples  B 
and  C. 


Figure  IV-60  indicates  comparable  noise  rejection  performance 
with  the  two  noise  samples  and  about  1 -db  measurable  deterioration  outside 
the  design  area  over  the  frequency  range  from  0.  0  to  5.  0  cps.  The  behavior 
of  each  MCF  is  markedly  similar,  with  differences  in  degree  of  rejection 
only. 
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Figure  IV-61  follows  the  same  trend  of  similarity  in  behavior 
and  about  1  -db  deterioration  in  performance  outside  the  design  interval. 

Also,  the  relative  performance  of  the  three  filters  and  the  straight  sum 
trace  are  compared.  Superior  performance  is  evident  for  the  filters  in  the 
region  of  greatest  noise  power  (below  1.  0  cps).  Above  1.  0  cps  MCF-15 
and  the  straight  sum  are  roughly  equivalent  while  MCF-14  and  MCF-13 
show  3-  to  5-db  less  rejection.  These  results  compare  favorably  with  the 
results  of  applying  the  filters  (Figure  IV-59)  and  the  random  noise  responses 
(Figure  IV-55).  However,  Figure  IV-54  indicates  that  the  actual  noise  re¬ 
jection  at  these  higher  frequencies  is  relatively  insignificant  when  compared 
with  the  rejection  at  frequencies  where  maximum  power  is  concentrated. 

In  summary,  the  effect  of  increasing  the  length  of  the  MCF 
is  apparently  an  increase  in  overall  noise  rejection  and  specifically  an 
improvement  in  the  higher  frequency  energy  rejection.  This  is  probably 
because  the  filter  design  constrains  each  MCF  to  reject  first  the  energy 
which  more  easily  minimizes  the  mean -square -error  (i.  e.  ,  the  large 
quantities  of  low-velocity  low-frequency  energy)  and  later  achieves  the 
relatively  minor  gains  possible  at  the  higher  frequencies.  The  gains  are 
minor,  however,  especially  for  the  nearly  6-to-l  increase  in  filter  length 
(in  one  instance),  and  generally  would  not  be  economically  justifiable.  The 
remarkable  stability  in  the  noise  rejection  as  a  frequency  function  indicates 
that  differences  in  MCF  performance  as  a  function  of  filter  length  are 
primarily  a  matter  of  degree  rather  than  emphasis. 
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Figure  IV-5.  Signal  Applications 
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Figure  IV-6. 
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Figure  IV-7.  Impulse  and  Frequency  Response  of  MCF-4 
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Figure  IV-8.  Impulse  and  Frequency  Response  of  MCF-5 


Figure  IV -9.  Impulse  and  Frequency  Response  of  MCF-6 
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Figure  IV  -  10. 


Impulse  and  Frequency  Response  of  MCF-7 
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Figure  IV -11.  Impulse  and  Frequency  Response  of  MCF-8 
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Figure  IV-12. 


Impulse  and  Frequency  Response  of  MCF-9 
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Figure  IV-13.  Impulse  and  Frequency  Response  of  MCF-10 
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Figure  IV- 14.  Percent  Mean-Square-Error 
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Figure  IV-15.  Percent  Mean-Square-Error 
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Figure  IV-16.  The  f-k  Responses  of  MCF-3 
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F  igure  IV -  17. 
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Figure  IV-18.  The  f-k  Responses  of  MCF-5 
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Figure  IV-19.  The  f-k  Responses  of  MCF-6 
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Figure  IV-20.  The  f-k  Responses  of  MCF-7 
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Figure  IV-21.  The  f-k  Responses  of  MCF-8 
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Figure  IV-22.  The  f-k  Responses  of  MCF-9 
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Figure  IV-23.  The  f-k  Responses  of  MCF-10 
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Figure  IV-24.  MCF-3  Output/Straight  Sum 
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Figure  IV-25.  MCF-4  Output/Straight  Sum 
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Figure  IV-26.  MCF-5  Output/Straight  Sum 
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Figure  IV-27.  MCF-6  Output /Straight  Sum 
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NOISE  SAMPLE  D 


Figure  IV-28.  MCF-7  Output/Straight  Sum 


Figure  IV-29.  MCF-8  Output /Straight  Sum 
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Figure  IV-30.  MCF-9  Output/Straight  Sum 
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Figure  IV-31.  MCF-10  Output/Straight  Sum 
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Figure  IV-32.  MCF-3  Output/CC 
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Figure  IV-33.  MCF-4  Output/CC 
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Figure  IV-34.  MCF-5  Output/CC 
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Figure  IV-35.  MCF-6  Output/CC 
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Figure  IV-36.  MCF-7  Output/CC 


Figure  IV -37 . 


MCF-8  Output/CC 
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NOISE  SAMPLE  E 


Figure  IV-38.  MCF-9  Output/CC 


NOISE  SAMPLE  E 


Figure  IV-39.  MCF-10  Output/CC 
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Figure  IV-40.  MCF-4  Output/Avg 


Figure  IV-41.  MCF-5  Output/Avg 
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Figure  IV-43.  MCF- 
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FREQUENCY  (CPS) 

Figure  IV-44.  Impulse,  -Amplitude  and  Phase  Responses  of  MCF-11  with 
Theoretical  Frequency-Domain  Shaping  Function 
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Figure  IV -46.  F requency- Wavenumber  Responses  of  MCF -  1 1;  f=  0.  25,  0.50, 
1 . 00,  and  1,50  cps 
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CENTER  SEISMOMETER  (RING  1) 
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Figure  IV-47.  Noise  Samples  A  and  D  Showing  Center  Seismometer,  Straight  Sum  Output 
and  MCF-6,  -11  and  -12  Outputs 
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Figure  IV-48.  Noise  Attenuation  as  a  Ratio  of  Power  Density  Spectra;  MCF-11  Output 
to  Center  Seismometer  and  to  Straight  Sum  Output  (Noise  Samples  A  and 
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F igure  IV - 50.  F requency-  Wavenumber  Responses  for  MCF -12;  f-  0.25,  0.50, 
1.0,  and  1.50 
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Figure  IV-51.  Noise  Attenuation  as  a  Ratio  of  Noise  Power  Density  Spectrum  of  MCF  Output 
to  Noise  Power  Density  Spectrum  of  Center  Channel  Input  (MCF/CC  Ratio) 
for  MCF  - 12 
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Figure  IV-52.  Percentage  Mean-Square-Error  as  a  Function  of  Filter  Length  for  the 
Design  of  MCF-13,  -14  and  -15 
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Figure  IV-53.  Impulse  Amplitude  and  Phase  Responses  of 
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Figure  IV-54.  Relative  Power  Density  Spectra:  Center  Seismometer 
Straight  Sum  Output,  Overlaid  with  MCF-13,  -14  and 
-15  Outputs 
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F igure  IV-5  5.  Random  Noise  Response  of  MCF  -  1  3,  -  14  and 
-15  as  a  Frequency  Function 
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Figure  IV-56.  Frequency -Wavenumber  Responses  for  MCF-13; 
f=0.25,  0.50,  1.0,  and  1.50 
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Figure  IV-57.  Frequency- Wavenumber  Responses  for  MCF-14;  f  =  0.25, 
0.  50,  1.  0,  and  1.  50 
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Figure  IV-58.  Frequency -Wavenumber  Responses  for  MCF-15;  f  =  0.  25, 
0  50,  1. 0,  and  1.50 
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Figure  IV-59.  Noise  Samples  B  and  C  Showing  Center  Seismometer,  Straight  Sum  Output, 
MCF-13,  -14  and  -15  Output 
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Figure  IV-60.  Noise  Attenuation  as  a  Ratio  of  Power  Density  Spectra, 
and  -15  Outputs  to  Center  Seismometer,  Noise  Sample 
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Figure  IV-61.  Noise  Attenuation  as  a  Ratio  of  Power  Density  Spectra, 
Outputs  to  Straight  Sum  Output,  Noise  Samples  B  and  C 


SECTION  V 


DESIGN  AND  EVALUATION  OF  MULTICHANNEL  FILTERS 
AND  THE  DETECTION  PROCESSING  PROCEDURE 

On  the  basis  of  the  optimum  design  techniques  determined 
by  the  study  described  in  the  previous  section,  one  MCF  was  developed 
for  each  of  five  LASA  subarrays.  This  section  presents  the  design  and 
evaluation  of  these  five  filters  (MCF-16  through  MCF-20)  and  a  brief 
discussion  of  their  capabilities  as  detection  processors.  Another  set  of 
filters  (MCF-21  through  MCF-25)  designed  with  a  different  signal  model 
are  also  briefly  described  in  this  section.  The  subarrays  for  which  the 
filters  were  designed  are  C-3,  E-2,  F-2,  F-3,  and  F-4  (Figure  V-l). 

A.  DESIGN  OF  THE  MULTICHANNEL  FILTERS 

The  type  of  filter  chosen  was  the  5-ring  filter  designed  from 
a  disk  signal  model  and  from  measured-noise  statistics.  As  was  shown 
in  Section  IV,  the  5 -ring  filter  performed  as  well  as  or  better  than  the 
8 -ring  or  the  25 -channel  maximum -likelihood  filter.  The  disk  signal 
model  was  shown  to  offer  some  improvement  over  the  infinite  velocity 
signal  model  in  maintaining  a  high,  almost  flat  response  to  signal  over 
the  velocity  range  from  12  km/sec  to  infinity.  The  superiority  of  the 
use  of  measur  ed -nois e  statistics  as  opposed  to  a  theoretical  noise  model 
was  shown  in  its  sharp  response  to  signal  in  frequency -wavenumber  space 
and  in  suppressing  the  noise  at  low  frequencies.  As  was  mentioned  before, 
the  averaging  of  the  correlation  sets  from  several  noise  samples  recorded 
on  different  days  eliminates  many  of  the  eccentricities  of  a  single  noise 
sample  and  produces  a  filter  which  operates  consistently  over  a  long  period 
of  time. 

In  order  to  improve  the  low-frequency  noise  attenuation, 
frequency -filtering  was  used  to  cut  off  all  information  below  a  certain 
frequency.  The  high-frequency  attenuation  was  improved  by  whitening 
the  noise  correlations  to  reduce  the  effect  of  the  large  low-frequency  noise 
peak  in  the  design  of  the  filters. 

The  study  of  the  effect  of  filter  length  indicated  that  the 
gain  of  a  5.  85-sec  length  over  a  2-sec  length  was  insufficient  to  justify 
its  use.  The  length  of  1.  8  sec  resulted  from  the  shortening  produced  in 
whitening  the  correlations  calculated  for  a  2-sec  filter  and  seemed  to  be 
a  good  compromise  between  the  1 -sec  and  2 -sec  filters  which  did  not 
differ  greatly  in  quality. 
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The  large  variations  in  seismometer  response  within  only  a 
few  months  led  to  the  decision  to  not  include  instrument  response  variations 
in  the  signal  model. 


The  filters  were  designed  from  an  average  of  several  noise 
correlations  sets  which  were  obtained  from  noise  recorded  simultaneously 
at  each  subarray;  however,  in  some  cases  the  record  of  a  particular  sample 
at  a  particular  subarray  was  considered  to  be  a  poor  speciman  and  was 
discarded.  Table  V-l  gives  the  noise  samples  that  were  used  in  the  design 
of  the  filters.  A  150-sec  gate  was  used  in  each  case,  the  data  being 
selected  generally  at  the  beginning  of  each  330 -sec  record. 

The  signal  model  corresponds  to  constant  power  density 
over  a  disk  in  wavenumber  space  (k-space)  with  a  radius  corresponding 
to  a  horizontal  velocity  component  of  12  km/sec.  The  power  density 
spectrum  was  shaped  by  the  autopower  density  spectrum  of  the  average  of 
the  noise  autocorrelations  of  a  single  subarray  to  provide  a  separate  signal 
model  for  each  subarray. 

Each  signal  model  was  then  reshaped  to  provide  the  specified 
frequency  filtering.  This  frequency  shaping  function,  used  for  all  five  sub¬ 
arrays  ,  was  flat  from  1.0  to  3.0  cps,  down  10  db  at  0.75  and  3.25  cps , 
down  30  at  0.  50  and  3.  50  cps,  down  50  db  at  3.  75  cps  and  equal  to  0.  0  at 
all  other  frequencies  from  0.00  to  10.00  cps.  The  combination  of  the  two 
shapings  produced  a  signal  and  noise  pair  with  constant  signal -to -nois  e 
ratio  between  1.00  and  3.00  cps  and  a  signal-to-noise  ratio  which  decreased 
rapidly  to  zero  outside  this  band.  The  ratio  of  total  signal  power  to  total 
noise  power  was  then  set  equal  to  2. 

The  second  set  of  filters  (MCF-21  through  MCF-25)  was 
designed  identically  as  the  first  set  except  that  the  signal  model  was  an 
infinite  velocity  model  rather  than  a  disk  model. 

B.  MULTICHANNEL  FILTER  RESPONSES 

The  responses  presented  for  each  filter  include  the  percent 
mean-square  error  as  a  function  of  filter  length  for  design  of  the  MCF, 
the  impulse  response  and  the  frequency  response  and  the  frequency -wave - 
number  response.  The  frequency  response  is  presented  in  two  parts:  the 
amplitude  response  and  the  phase  response. 

The  percent  mean -square -er ror  as  a  function  of  filter  length 
(Figures  V-2  through  V-6)  shows  a  total  noise  power  reduction  of  96  percent 
to  98  percent  for  each  filter.  The  impulse  responses  of  all  six  filters 
(Figures  V-2  through  V-6)  are  almost  identical.  The  amplitude  responses 
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Table  V-l 


NOISE  SAMPLES  USED  IN  DESIGN  OF  MCF-I6  THROUGH  MCF-20 


Subarray 

MCF  No. 

Date 

Time  (MST) 

11-05-65 

03:19-03:25 

1  1  -09-65 

02:39-02:45 

1 1-16-65 

06:29-06:35 

11-20-65 

03:40-03:46 

1 1  -23-65 

02:19-02:24 

11-23-65 

05:49-05:55 

1 1  -24-65 

03:18-03:24 

12-01-65 

10:10-10:15 

03-24-66 

12:47-12:55 

03-29-66 

11:08-11:16 

03-31-66 

13:58-14:06 

04-05-66 

16:10-16:18 

C3 

E2 

F2 

F  3 

F  4 

16 

17 

18 

19 

20 

X 


X 


X 


X  -  Indicates  that  a  sample  from  a  particular 
subarray  was  not  used. 
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of  all  the  filters  are  similar:  showing  a  good  passband;  almost  flat  from 
1.  0  to  3.  0  cps,  as  desired;  and  a  drop  to  about  -20  to  -30  db  above  4.  0 
cps.  The  phase  responses  are  flat  from  0.  5  to  3.  5  cps  and,  in  general, 
oscillate  between  0°  ±  5°  and  180°  ±  5°  outside  this  frequency  band. 

The  second  set  of  filters  (MCF-21  through  MCF-25)  was 
designed  to  attempt  a  slight  increase  in  noise  rejection  at  the  expense  of 
low  velocity  signal  response.  As  in  the  case  of  the  disk  signal  model 
filters,  the  infinite  velocity  filters  have  nearly  identical  responses,  so 
the  responses  are  presented  for  only  one  filter,  MCF-21,  designed  for 
the  F-3  subarray.  The  mean -square -er ror  (Figure  V-7)  has  been  reduced 
to  1.  1  percent  compared  with  2.  8  percent  for  MCF-16  (the  disk  signal 
model  filter  designed  for  the  F-3  subarray).  The  impulse  responses  of 
these  two  filters  are  similar  as  are  the  frequency  responses.  One 
difference  in  the  amplitude  responses  is  that  at  high  frequencies  (above 
4.  0  cps)  the  response  of  MCF-21  is  down  about  30  db  compared  to 
approximately  22  db  for  the  MCF-16  response. 

The  frequency -wavenumber  responses  of  the  five  filters 
are  shown  in  Figures  V-8  through  V-12.  As  expected  from  examination 
of  the  other  responses,  there  is  little  difference  among  the  f-k  responses. 
The  amplitude  responses  dropped  very  low  below  1.  0  cps,  but  the  lowest 
point  in  that  band  occurred  at  0.  30  cps.  This  is  why  the  amplitude 
responses  indicate  a  greater  rejection  than  the  f-k  responses  show  at 
infinite  velocity  at  0.25  cps.  The  difference  of  5  to  10  db  is  easily 
accounted  for  by  the  extreme  slope  of  the  amplitude  responses  at  these 
frequencies.  At  1.  00  cps,  the  f-k  responses  are  up  to  0  db  at  infinite 
velocity  and  are  down  less  than  1.  0  to  1.  5  db  at  12  km/sec  at  this  fre¬ 
quency  and  at  1.  50  cps.  At  velocities  less  than  12  km/sec  the  responses 
fall  sharply  to  levels  below  -12  db.  The  exact  level  and  velocity  at  which 
greatest  rejection  occurs  varies  from  filter  to  filter  since  each  filter 
was  designed  for  a  different  subarray.  At  0.25  cps  a  rejection  ring 
occurs  between  2  and  3  km/sec  with  about  12  to  15  db  more  rejection 
than  at  infinite  velocity.  At  0.  50  cps  the  response  has  risen  several 
decibels  while  maintaining  the  same  general  shape.  At  1.  00  cps  most 
of  the  filters  are  rejecting  noise  most  strongly  around  2.  5  km/sec  and 
1.  5  km/sec,  with  typical  rejections  of  18  db  or  more.  At  1.  50  cps  these 
rejection  rings  are  still  present  but  are  not  as  deep  as  at  1.  00  cps. 

The  f-k  response  of  MCF-21  is  shown  in  Figure  V-13.  This 
response  is  basically  the  same  as  that  of  MCF-16  except  for  two  differences. 
MCF-21  tends  to  have  3  to  6  db  greater  rejection  at  low  velocities  than 
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MCF-16.  At  high  velocities  MCF-21  does  not  have  the  flat  signal  response 
which  MCF-16  has  at  velocities  above  12  km/ sec.  At  1.00  cps  the  response 
of  MCF-21  is  down  about  4  db  at  12  km/sec.  Thus  for  detection  of  signals 
at  high  velocities  above  20  km/sec,  the  infinite  velocity  filter  would  be  as 
good  as  or  better  than  the  disk  signal  model  filter. 

The  general  characteristics  which  were  enhanced  by  proper 
selection  of  design  parameters  are  noticed  in  the  f-k  responses: 

•  High  rejection  over  all  f-k  space  at  low  frequencies 

•  A  nearly  flat  pass  region  in  k-space  corresponding 
to  high  velocity  signals  at  frequencies  from  1. 00  to 
3. 00  cps 

•  Sharp  falloff  of  the  response  at  velocities  below  12  km/sec 

•  Good  rejection  in  all  directions  at  low  velocities 

•  A  generally  isotropic  response 

C.  NOISE  ATTENUATION  AND  SIGNAL  DISTORTION 

The  plots  of  noise  attenuation  are  shown  in  Figures  V-14 
through  V-16.  These  plots  show  the  ratio  of  the  MCF  output  power  density 
spectrum  to  the  straight  summation  power  density  spectrum  for  the  noise 
immediately  preceding  one  of  the  signals  used  in  the  evaluation  of  the 
application  of  the  filters  to  signals.  The  three  signals  are  designated  by 
the  location  of  their  source  as  Algeria,  Kuriles  and  Marianas. 

The  noise  used  in  the  plots  of  noise  attenuation  was  recorded 
at  the  subarray  for  which  the  MCF  applied  to  that  noise  was  designed.  The 
subarray  is  identified  in  parentheses  on  the  figure. 

The  ratios  show  that  the  filters  have  greater  attenuation  of 
noise  than  the  straight  sum  process  at  all  frequencies  below  0.  9  or  1.0  cps. 
The  maximum  improvement  of  the  filters,  compared  to  the  straight  sum, 
occurs  between  0.  25  and  0.  30  cps  and  is  about  25  db  in  most  cases;  going 
as  high  as  30  db  in  the  case  of  MCF-20  applied  to  noise  recorded  at  F -2 
preceding  the  signal  from  Algeria.  The  improvement  goes  as  low  as  2ldb 
in  the  case  of  MCF-20  applied  to  noise  recorded  at  F -2  preceding  the  signal 
from  the  Marianas.  From  1.0  to  about  4.0  cps,  the  MCFs  have  about  3  to 
lOdb  less  noise  rejection  than  the  straight  sum  process.  This  difference 
is  probably  due  to  the  high  proportion  of  random  noise  in  this  band,  little  of 
which  the  MCF  can  reject  when  the  signal  model  is  a  disk  in  k-space .  At  higher 
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frequencies  the  filters  noise  rejection  approaches  that  of  the  straight 
sum.  At  these  frequencies  the  MCF  has  allowed  the  signal  response 
to  drop  to  a  low  value  because  of  the  signal  response  shaping  function 
used  in  the  original  design  of  the  MCFs.  It  should  be  remembered 
that  for  the  25 - seismomete r  LASA  subarrays,  the  straight  sum  process 
has  a  theoretical  maximum  noise  rejection  of  14  db. 

The  MCF  outputs  for  samples  containing  signals  are  shown 
in  Figures  V-17,  V-18  and  V-19.  The  signals  are  from  Algeria,  the 
Kurile  Islands  and  the  Mariana  Islands,  respectively.  The  figures  show 
no  significant  signal  distortion.  A  slight  precursor  can  be  seen  on  some 
of  the  figures,  but  it  is  very  small  and  doesn't  present  any  problem  in 
signal  detection. 

D.  DETECTION  PROCESSING 

The  detection  processing  procedure  consisted  of  a  running 
autocovariance  estimate  over  the  MCF  output.  For  a  trace  G  with  elements 
g. ,  A.  was  computed  for  every  value  of  i  according  to  the  equation 

1  k  =  N’1  9 

Ai  =  "n~~  53  gi  -k 

k=0 

The  gate  used  in  these  cases  was  N  =  50  or  a  time  of  2.  5  sec.  The  outputs 
are  shown  in  Figures  V-20,  V-21  and  V-22.  The  first  signal  arrival  is 
clearly  shown  as  are  some  of  the  later  phase  arrivals.  The  peak  occuring 
11  sec  after  the  P-wave  peak  for  the  Kurile  Islands  event  corresponds  to  a 
pP  arrival  for  the  reported  33  km  depth.  The  peak  occuring  40  sec  after 
the  P-wave  peak  for  the  Mariana  Islands  is  unidentified.  The  reported 
depth  for  this  event  was  also  33  km.  Strong  evidence  of  teleseismic  energy 
was  obtained  in  each  case.  The  indication  at  subarray  F3  of  the  Mariana 
Islands  event  is  especially  impressive  considering  the  poor  signal -to -nois e 
ratio  of  this  event  at  F3 . 
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Figure  V-l.  LASA  Subarrays 
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Figure  V-2.  Mean -Square -Error  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF -16 
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PERCENT  MEAN  SQUARE  ERROR 


Figure  V-3.  Mean -Square -Err or  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF -17 
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PERCENT  MEAN  SQUARE  ERROR 


Figure  V-4.  Mean -Square -Error  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF-18 


V-10 


science  services  division 


PERCENT  MEAN  SQUARE  ERROR 


Figure  V -5.  Mean-Square-Error  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF-19 
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PERCENT  MEAN  SQUARE  ERROR 


Figure  V-6.  Mean -Square -Error  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF-20 
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Figure  V-7.  Mean -Square -Error  Vs  Filter  Length,  Impulse 
Response  and  Frequency  Response  of  MCF-21 
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Figure  V-8.  Frequency -Wavenumber  Response  of  MCF  -  1  6 ;  f  =  0.  25,  0.  50, 
1.00  and  1 .  50  cps 
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Figure  V-9.  Frequency-Wavenumber  Response  of  MCF-17;  f  =  0.  25,  0.  50, 
1. 00  and  1.  50  cps 


V  - 1  5 


science  services  division 


Figure  V-10.  Frequency -Wavenumber  Response  of  MCF-18;  f  =  0.25,  0.50, 
1.00  and  1 .  50  cps 
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Figure  V-ll.  Frequency -Wavenumber  Response  of  MCF - 1  9;  f  =  0.  25,  0.50, 
1 .  00  and  1 .  50  cps 
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Figure  V-1Z.  Frequency-Wavenumber  Response  of  MCF-20;  f  =  0.  25,  0.50, 
1.00  and  1.50  cps 
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Figure  V-13.  Frequency -Wavenumber  Response  of  MCF-21;  f  =  0.  25,  0.  50,  1. 00 
and  1 .  50  cps 
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Figure  V-14.  Ratio  of  MCF  Output  Power  Density  Spectrum  to  Straight  Sum  Power 
Density  Spectrum  for  Noise  Immediately  Preceding  the  Algerian 
Signal.  Subarray  at  Which  Noise  Was  Recorded  is  Indicated 
in  Parentheses 
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Figure  V-15.  Ratio  of  MCF  Output  Power  Density  Spectrum  to  Straight  Sum  Power 
Density  Spectrum  for  Noise  Immediately  Preceding  the  Kuriles 
Signal.  Subarray  at  Which  Noise  Was  Recorded  is  Indicated 
in  Parentheses 
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Figure  V-16.  Ratio  of  MCF  Output  Power  Density  Spectrum  to  Straight  Sum  Power 
Density  Spectrum  for  Noise  Immediately  Preceding  the  Marianas 
Signal.  Subarray  at  Which  Noise  Was  Recorded  is  Indicated 
in  Parentheses 
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Figure  V-17a.  MCF  Outputs  and  Reference  Traces  for  Signal  from 

Algeria.  The  outputs  are  Shown  for  Each  Filter  Applied 
to  the  Signal  From  Each  Subarray.  The  Numbers  in 
Parentheses  Indicate  the  Subarray  for  Which  the  Filter 
was  Designed 
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Figure  V-17b.  MCF  Outputs  and  Reference  Traces  for  Signal  from  Algeria.  The  Outputs  are  Shown  for 
Each  Filter  Applied  to  the  Signal  from  Each  Subarray.  The  Numbers  in  Parentheses 
Indicate  the  Subarray  for  Which  the  Filter  was  Designed 
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Figure  V-18.  MCF  Outputs  and  Reference  Traces  for  Signal  from  the  Kurile 
Islands.  The  Results  are  Shown  Only  for  the  MCF  Applied  to 
Data  for  the  Subarray  for  Which  it  was  Designed 
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-19.  MCF  Outputs  and  Reference  Traces  for  Signal  from  the  Mariana 
Islands.  The  Results  are  Shown  Only  for  the  MCF  Applied  to 
Data  for  the  Subarray  for  Which  it  was  Designed 
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Figure. V-20.  Outputs  of  Detection  Processor  for  Each  Subarray  for  Signal 
from  Algeria 
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APPENDIX  A 


AVERAGED  NOISE  CORRELATION  SETS 
FOR  THE  FIVE  SUBARRAYS 


The  noise  models  used  in  the  design  of  the  filters  presented  in 
Section  V  were  formed  by  averaging  several  ring  model  correlation  sets. 

The  recording  date  and  time  of  the  noise  sample  for  each  cor  relation  set 
is  listed  in  Table  V-l.  The  list  consists  of  seven  night  noise  and  five  day 
noise  samples.  Therefore,  to  equally  represent  day  and  night  noise,  the 
day  noise  samples  are  multiplied  by  a  factor  of  7/5. 

Two  night  noise  samples  (5  Nov  1965  and  l6  Nov  1965)  at  subarray 
F-3  and  one  night  noise  sample  (23  Nov  1965)  at  subarray  E-2  were  found 
to  be  unusable  due  to  severe  non-seismic  noise  bursts.  These  samples 
are  excluded  fiom  the  average. 

Figures  A-la  through  A-5b  present  the  ave raged-noise  correlation 
sets  for  the  five  subarrays.  The  correlations  are  in  the  order  ( ring)  1* ( ring)l, 
1*2,  .  •  •  1*5,  2*2,  ...2*5,  3*3,  ...3*5,  4*4,  4*5,  5*5,  for  each  subarray. 

The  crosscorrelations  are  shown  from  -2.0  sec  to  +2.0  sec,  and  the  auto¬ 
correlations  between  0.  0  sec  and  2.  0  sec.  Each  lag  value  represents  0.  05 
sec  in  time. 

Whitening  filters  were  developed  for  each  subarray  using  the  av¬ 
eraged  autocorrelation  of  the  center  seismometer  (ringl).  These  filters 
are  presented  in  Table  A-l.  Each  filter  was  autocor related  and  subsequently 
convolved  with  the  averaged  correlations  of  that  set.  The  resulting  sets  of 
whitened  correlations  were  truncated  by  the  convolution  process  to  between 
-1.85  sec  and  +1.85  sec  in  length,  and  were  used  in  the  MCF  design. 

The  inclusion  in  this  report  of  the  listings  of  the  unwhitened  cor¬ 
relation  sets  is  based  on  the  belief  that  they  would  be  of  considerable  use 
to  those  interested  in  further  study  along  these  lines.  Much  expensive  and 
relatively  non- technical  effort  was  required  to  assemble,  demultiplex, 
ring-sum,  compute  and  average  correlations  with  such  a  considerable 
quantity  of  data.  The  availability  of  this  data  may  expedite  further  signal 
detection  studies. 
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Figure  A -la.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  F-3 
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Figure  A-lb.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  F -3 
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Ring  Model  of  Subarray  F  - 
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Figure  A-3a.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  F-4 
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Figure  A-3b.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  F -4 
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Figure  A-4a.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  for  Subarray  C-3 
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Figure  A-4b.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  C-3 
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Figure  A -5a.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  for  Subarray  E-2 
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Figure  A-5b.  Averaged  (Day  plus  Night)  Noise  Correlation  Set  for  5-Channel 
Ring  Model  of  Subarray  E-2 


Table  A-l 

WHITENING  FILTERS  FOR  THE  FIVE  SUBARRAYS 


Subarray  C -3 

Subarray  F -2 

1 

1. 0000000E  00 

1 

1. 0000000E  00 

2 

-1. 0757737E  00 

2 

-1.  201 1  120E  00 

3 

-6.  4689597E-02 

3 

3.  5605594E-01 

4 

9.  6869609E-02 

4 

-4.  5989817E-02 

5 

7.  2587480E -02 

5 

-7.  7634451E  -02 

Subarray  F -4 

Subarray  E -2 

1 

1. 0000000E  00 

1 

1 . 0000000E  00 

2 

-1. 1068382E  00 

2 

-1.  3222608E  00 

3 

-3. 3041 908E  -02 

3 

5.  51081  98E  -0 1 

4 

9.  9640668E-02 

4 

-8.  9560338E-02 

5 

7.  2090762E-02 

5 

-9.  5849597E-02 

Subarray  F -3 

1 

1 . 0000000E  00 

2 

-1.  227 641 8E  00 

3 

6.  4279133E-01 

4 

-3.  52050C6E -01 

5 

-2.  6365489E-02 
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APPENDIX  B 

MULTICHANNEL  FILTER  WEIGHTS 

This  report  has  dealt  with  the  design  and  evaluation  of  23  multi¬ 
channel  filters  of  various  types  and  descriptions.  Of  these,  13  were  de¬ 
scribed  and  compared  in  various  ways  in  Section  IV.  They  are  briefly 
outlined  in  Table  IV-1,  presented  earlier.  The  remaining  10  filters 
were  presented  in  Section  V  as  ,,finalM  filters  since  they  incorporated  the 
more  desirable  features  revealed  by  the  Section  IV  comparisons.  This 
group  consists  of  five  disk  model  filters  and  five  inf inite -velocity  filters, 
one  for  each  of  the  five  subarrays. 

In  anticipation  of  further  study  and  to  complete  the  MCF  presen¬ 
tation,  Appendix  B  presents  a  complete  listing  of  all  MCF  filter  weights. 
These  are  contained  in  Figures  B-l  through  B-23. 

MCF -  3  is  a  25-channel  filter,  MCF -4,  -5  and  -8  are  8-channel 
ring  models,  with  all  remaining  filters  corresponding  to  a  5-channel  ring 
model  configuration.  In  each  instance,  the  channel  1  designation  refers  to 
the  center  seismometer.  For  MCF -  3,  the  filters  are  presented  in  the 
order  corresponding  to  seismometers  10,  21,  31,  51,  71,  22,  42,  62,  82, 
23,  33,  53,  73,  24,  44,  64,  84,  25,  35,  55,  75,  26,  46,  66,  and  86. 
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Figure  B-la.  Filter  Weight  Listings  for  MCF-3;  25-Channel 
Maximum -Likelihood  MCF 
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Figure  B-lb.  Filter  Weight  Listings  for  MCF -3;  25-Channel 
Maximum-Likelihood  MCF 
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Figure  B-lc.  Filter  Weights  for  MCF  -3;  25-Channel 
Maximum- Likelihood  MCF 
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Figure  B-2.  Filter  Weight  Listings  for  MCF-4;  8-Ring  IV  Model  MCF 
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Figure  B-3.  Filter  Weight  Listings  for  MCF-5;  8-Ring  Disk  Model  MCF 
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Figure  B-5.  Filter  Weight  Listings  for  MCF-7;  5-Ring  Disk  Model  MCF 
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Figure  B-6.  Filter  Weight  Listings  for  MCF -8;  8-Ring  IV  Model  MCF 
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Figure  B-8.  Filter  Weight  Listings  for  MCF-10;  5-Ring  IV  Model  MCF 
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Figure  B-10.  Filter  Weight  Listings  for  MCF-12;  5-Ring  Disk  Model  MCF 
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Figure  B-ll.  Filter  Weight  Listings  for  MCF -13;  5-Ring  IV  Model  MCF 
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Figure  B-12.  Filter  Weight  Listings  for  MCF-14;  5-Ring  IV  Model  MCF 
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Figure  B-  13. 
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Figure  B-17.  Filter  Weight  Listings  for  MCF-16;  Final  5-Channel 
Disk  Model;  Subarray  F-3 
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Figure  B-14.  Filter  Weight  Listings  for  MCF-18;  Final  5-Channel  Disk 
Model;  Subarray  C-3 
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Figure  B-16.  Filter  Weight  Listings  for  MCF-20;  Final  5-Channel 
Disk  Model;  Subarray  F -2 
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Figure  B-22.  Filter  Weight  Listings  for  MCF-21;  Final  5-Channel 
Infinite  Velocity  Model;  Subarray  F-3 
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Figure  B-23.  Filter  Weight  Listings  for  MCF-22:  Final  5-Channel 
Infinite  Velocity  Model;  Subarray  F -4 
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Figure  B-19.  Filter  Weight  Listings  for  MCF-23:  Final  5-Channel 
Infinite  Velocity  Model;  Subarray  C-3 
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Figure  B-20.  Filter  Weight  Listings  for  MCF-24;  Final  5-Channel 
Infinite  Velocity  Model;  Subarray  E-2 
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Figure  B-21.  Filter  Weight  Listings  for  MCF-25;  Final  5-Channel 
Infinite  Velocity  Model;  Subarray  F-2 
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Figure  C-l.  Relative  Amplitude  Responses  of  Subarray  F- 
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Figure  C-2.  Relative  Amplitude  Responses  of  Subarray  F- 
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Figure  C-3.  Relative  Amplitude  Responses  of  Subarray  C 
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Figure  C-4.  Relative  Amplitude  Responses  of  Subarray  E-2 
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Figure  C-5.  Relative  Amplitude  Responses  of  Subarray  F 
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Figure  C-10.  Relative  Phase  Responses  of  Subarray  F 
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Figure  C-ll.  Relative  Phase  and  Amplitude  Responses  for  Subarray  F-3  at 
Frequencies  0.202  and  0.802  cps 
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Figure  C-12.  Relative  Phase  and  Amplitude  Responses  for  Subarray  F-3  at 
Frequencies  1.  198  and  1.418  cps 
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igure  C-13.  Relative  Phase  and  Amplitude  Responses  for  Subarray  F— 3  at 
Frequencies  2.  016  and  3.  03  cps 
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